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GENERAL  CIRCULATION  EXPERIMENTS 
WITH  A  TWO  LEVEL  PRIMITIVE  EQUATIONS  MODEL 

b  y 

Eliahu  Doron  and  Abraham  IIuss 

A  bstract 

A  modified  version  of  the  primitive  equations  model  proposed  by 
Mintz  and  Arakawa  was  integrated  in  order  to  obtain  circulation  patterns  in 
a  zonal  strip.  A  hydrologic  cycle  was  included  and  its  influence  on  the  develop¬ 
ing  circulation  was  investigated.  In  the  first  integration  not  iucluding  release 
of  latent  heat  a  more  or  less  symmetric  circulation  developed,  with  a  jet  stream 
characteristic  of  summer  conditions  at  the  central  latitude.  The  strength  and 
distribution  of  the  meridional  circulation  was  also  quite  characteristic  of  summer 
in  the  real  atmosphere.  In  the  second  integration,  including  the  water  vapor 
cycle,  the  development  was  more  vigorous  and  higher  amplitudes  of  the  geo¬ 
potential  were  reached.  The  surface  pressure  gradients  were  also  more  pro¬ 
nounced  and  were  close  to  those  observed  in  the  southern  hemisphere  in  summer. 
The  resulting  precipitation  and  evaporation  amounts  were  much  too  low,  so  that 
although  an  improvement  in  the  direction  towards  a  more  realistic  simulation  of 
the  atmosphere  is  observed  in  the  second  integration,  this  improvement  is  not 
very  marked. 
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)  NTKODUCTIOW 


The  experiments  described  in  the  following  are  part  of  a 
series,  carried  out  in  order  to  assess  the  effect  of  the  hydrologic  cycle  on 
the  development  of  the  general  circulation  of  the  atmosphere.  They  are  a 
direct  continuation  of  research  carried  out  along  the  same  lines  with  a  quasi  - 
gecstrophic  model  (Huss  and  Doron,  1967).  From  the  results  of  this  earlier 
set  of  experiments  tentative  conclusions  were  drawn  as  to  the  effect  of  the 
release  of  latent  heat.  The  question,  as  to  whether  these  effects  are  real 
physical  phenomena  or  merely  derived  from  the  specific  mathematical  form¬ 
ulation  of  the  modelycould  not  be  answered  with  certainty.  One  way  of  over¬ 
coming  this  difficulty  is  to  try  to  simulate  the  same  effects  with  a  different 
model.  For  such  a  comparison  it  is  necessary  to  apply  both  models  to  the 
same  domain  and  to  include  the  same  physical  parameters,  as  far  as  possible. 

A  first  step  in  this  direction  was  made  in  the  following  ex¬ 
periments.  The  model  chosen  was  essentially  the  two  level,  primitive  equa¬ 
tions  model  developed  by  Mintz  and  Arakawa  (Langlois  and  Kwok  1969).  As 
this  model  covers  laterally  the  whole  globe  and  vertically  only  the  troposphere, 
certain  modifications  were  required.  These  included  the  application  of  a  different 
vertical  coordinate  and  the  formication  of  boundary  conditions  appropriate  for 
the  zonal  strip  chosen  for  the  quasi-geostrophic  model.  In  addition  the  original 
definition  of  heat  sources  and  sinks  was  abandoned  and  the  formulation  of  the 
quasi-geostrophic  case  was  used.  Finally,  all  the  effects  due  to  the  curvature 
of  the  globe  were  neglected  and  the  equations  were  applied  to  a  plane  strip,  as 
in  the  previous  model. 
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General  Aspects  of  the  Model 

In  order  to  investigate  further  the  influence  of  the  introduction 
of  water  vapor  and  release  of  latent  heat  on  the  general  circulation,  a  second 
more  advanced  model  was  used.  The  model  chosen  was  basically  the  two  level 
primitive  equations  model  developed  by  Mintz  and  Arakawa  (Langlois  and 
Kwok,  1968).  Certain  modifications,  which  will  be  listed  below,  were  nec¬ 
essary  in  order  to  apply  it  to  the  same  domain  chosen  for  the  previous  ex¬ 
periments.  Thus  the  vertical  coordinate  had  to  be  modified.  Furthermore, 
the  introduction  of  lateral  "walls"  as  boundaries  necessitated  the  introduc¬ 
tion  of  suitable  boundary  conditions  and  an  appropriate  definition  of  a  corres¬ 
ponding  finite  difference  scheme.  No  attempt  was  made  of  mapping  and  the 
equations  were  applied  to  a  plane.  Finally  the  original  formulation  of  heat 
sources  and  sinks  was  abandoned  and  a  simple  heating  function  was  employed, 
similar  to  the  one  used  in  the  previous  model.  Thus,  using  also  similar  initial 
conditions,  the  results  of  the  two  models  might  be  compared  The  main  ad¬ 
vantage  of  the  present  model  is  the  removal  of  the  quasi-geostrophic  restric¬ 
tion  imposed  on  the  previous  model.  A  second  advantage  would  be  that  oro¬ 
graphic  effects  could  be  introduced  in  a  straightforward  way  -  but  this  was 
not  done  in  the  following  experiments.  As  in  the  previous  model  two  levels 
were  used  in  the  vertical,  and  they  were  fixed  in  such  a  way  as  to  be  as  close 
as  possible  to  the  levels  used  previously  In  the  original  model  the  humidity 
field  was  carried  by  the  lower  level  only.  In  the  present  experiments,  water 
vapor  was  carried  explicitly  at  two  levels.  It  was  not  expected„however,  that 
this  would  lead  to  any  marked  improvements,  as  the  amount  of  water  at  the 
upper  level  must  be  at  least  one  order  of  magnitude  smaller  than  the  amount 
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at  the  lower  level.  Thus  the  release  of  latent  heat  at  the  upper  level  will 
have  only  a  small  effect.  It  was  expected  however,  that  this  modification 
would  improve  the  mechanism  of  water'vapor  transports  in  the  atmosphere. 
The  experience  gained  may  be  of  help  in  the  treatment  of  the  distribution 
of  other  materials  in  tne  atmosphere,  such  as  radioactive  substances. 

Governing  equations  and  integral  constraints 


Making  use  of  the  x,  y,  <r  ,  t  coordinate  system  introduced  by 
Phillips (1957)  the  thermo-hydrodynamic  equations  may  be  written  in  the 
form: 


■ri  v\/ 


+-  v  (7  v  +■  o~ 

^cr 


^  Ik  <  V  +■  ^  (f> 


T  O)  J; 

rr  Tier 


4-  V- 


fl) 


•+■  p-  T  o 

Tur  i  r  ~  1 


(2) 


(3) 


+  v  ,  y/Q  -v  c-  ^4  -  a/(r 

° 1  'Str  cp  "  p.^y 


o 

t 


(4) 


Cji 


C-& 

v  ITO"  • 


R/c 


P  , 


-  Q 


(5) 


where: 


v  -  Horizontal  wind  vector 
2.  -  -  Coriolis  parameter 

^  -  Vertical  unit  vector^ 

(b  -  Geopotential  / 

O  -  Potential  temperature  } 

-  Temperature 

R  -  Gas  constant  for  dry  air  , 

/ 

Cp  -  Specific  heat  at  constant  pressure 

P  -  100  cb 

o  ' 

"'T  -  Surface  pressure  . 

"  -  -±=  The  vertical  coordinate . 

rr  ;  / 

Q  r  Diabatic  heating  per  unit  mass/ 

F  -  Frictional  force  , 

Eqs.  ( 1  )  and  {  4  )  may  be  written  in  the  flux  form  by  multiply¬ 
ing  them  byrr  and  adding  (  3  )  multiplied  by  V  and  0  respectively. 

The  result  is: 

~:(Trv)  4-  V.^TTw)  4-  Tr|y<xv_  Vrr  (6) 

1  -  V  •  tt  cb  -  O 

^pr©)  +  v  +  ^Crrcro)  -  TlQ-/  ttctn^p  b  (7) 

•“p  po  ' 


6 


liy  substituting  the  temperature  for  the  potential  temperature,  eq.  (7) 
becomes: 


TTTcr;  ,  l'-ZZ£.  (8) 

fi._T  -  v  .  7tt)  -  'n_9r ,,  o. 

In  the  absence  of  frictional  forces  and  diabatic  heat  sources  the 
total  energy  -L  rr  \yz  +  rrcp7  .  the  mass-averaged  potential  temperature 

Tr  C--  ,  and  the  mass  averaged  second  moment  of  the  potential  temperature 

4  it  should  be  conserved.  Following  Mintz  and  Arakawa  our  aim 

is  to  construct  a  two  le  el  model  involving  a  vertical  finite  difference 
scheme,  which  will  maintain  the  above  mentioned  constraints.  As  the  der¬ 
ivation  and  the  final  results  are  somewhat  different  from  the  one  given  by 
Mintz  and  Arakawa.  they  will  be  discussed  in  the  following. 

The  levels  =  0,  1/4,  1/2,  3/4,  1  will  be  denoted  by  sub¬ 
scripts  0,  1,  2,  3,  4,  respectively.  Equations  (3),  (6),  (8),  will  be  ap¬ 
plied  at  levels  1  and  3.  Applying  eq.  (3)  at  both  levels  and  assuming: 


C 


-  2.7,  (_<r;  -  ‘•■Tc 


=  It r 


with  the  boundary  condition 
the  result  is  : 


9.TT  <ra 


; 


4-  ^  bTT  \V  \ 
-D-V  '  ' 


X—cn,  ~  c 


121  +  <7  . 

'  J-r 


(9) 


Next  we  write  the  last  term  of  eq.  (1 )  in  the  form: 

-  *1, 

and  applying  eq.  (1 )  at  both  levels  the  result  is: 

~  -t-  v,  «7v,  +  -V.)  -v-  |V*v,  *  4-  ~  <f>a)  =  O  , 

~i  -u  ^  -v<h.  iy%  -v.)  +■  $  ~  c  £3  ■*•$*)  =  °  .  (12 ) 

The  terms  ^,(.^3  *n  e(ls-  (12)  are  derived  by  interpolating  <r 
at  both  levels  in  the  form:  cr  -  and  o\  -  Sj. 

'CL  a. 

and  by  taking  forward  differencing  for  at  level  1  and  back¬ 

ward  differencing  at  level  3. 

Multiplying  eqs.  (12)  by  T TV,  and  tt  Vj  respectively: 

and  eqs.  (9)  by  ^  'y,2-  and  ^  respectively,  and  adding, 


we  obtain- 
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^UnV,J)  +  V-(tTT V,aV|)  frr^V.V^  +TTvf.O^|+.C<|(.^Vi  s 


o 


'  •'•  <?  it^V,lv4)-*'  ^W. -Vi  +■  TTVa.v7^a  - 


(13) 


The  fourth  terms  in  eqs.  (13)  can  be  written  in  the  form  (using  eq.  (9)): 


rxy,.C>4-|  r  V  -  4\  ^.(rrv()  -  7  .  (TT^>  V,)  4-  ^  *  2tr<M 

"V 

r-  ^  ■  TT  ^  ^  9  *  (TT  Va)  =  7  (tT^  Vaj  +  <|»a(gf  -  2ncrJ  # 

Applying  this  to  eq.  (13)  the  result  is  : 

y*[Trv.Ci^-tt>,)]  +(v,-v9  = 

-  *  v  • Qvr  +  lTr;~v  . 


o 

■>* 


2H  - 

'3-  •‘Si 


?*U  n  *&  )  +  y  ‘  C17  v*(  i  V^<t>a)j  -  C  V,  •  Wa,  r  *  TT<rx  -  fa 

s  “^x 1  ^)iM ^  >9rr)  +(<h-i)^ 

(14) 

Adding  eqs.  (14 )  and  integrating  over  the  horizontal  extent  of  the  domain 
we  get: 


r^Trtw,1  ^V^)j  ~  -  L4>1^-  <{>,)  ~£  +  ^J.CVV>  -  v*)  -  C7> 1 

-  iTTc^  -^y,-<7Tr  -<J>3Vv?TT 


(15) 


where  ( 


)  denotes  horizontal  averaging. 


Eq.  (15)  is  the  finite  difference  analogue  of  the  kinetic  energy 
equation,  when  integrated  over  the  domain 

We  now  write  eq.  (7)  at  both  levels  in  the  form: 

V7<(_TT©tVt)  +  CL'-x  G>2>a  -  o 

2_(-rr©,\  -+  <7-(XGav3)-  aire^  =0*  (16) 

o-t  ^  J  *  j 

adding  eq.  (16)  and  integrating  horizontally  over  the  domain  : 

^  v©a)j  =  O,  (17) 

Eq.  (17)  is  the  finite  difference  analogue  of  the  potential  temperature 
average,  so  that  this  constraint  is  satisfied. 

Adding  a^^T^&^and-  to  eqs.  (16) 

respectively  and  dividing  them  by  cp  and 

respectively  we  obtain  : 

2.hr<PTx)  +  ^7-Cn-CpT,  v,)  v  Xvr^cpT^  -  CpTr  ©,  %\i  &)J ^ 

^  * 

-cPire,v,  •?[(-&) “'“"I  =  Ira.-(&),>/<'eJ( 

^.C^PU)  -v-^OcpTa  V*)  ~  iTT^CpTa.  -  Cp  TT©^  t%)VlPJ  (18) 

-  CpT©^  =  lL^,)a/tp  ©x  "T*L 
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Now. 


2_ 


'  r  !  a  o  ~Qfcr~[  IGj  =  —  ( Cl  \RV  2IT 

rt  L  ^  k>)  .r  vtAc  r.  ^  s»pAiw  ^  ^ 

so  that  • 

--p-rrs,  --  -RT'  ^  / 

in  the  same  way: 

•  ipTres%[iy^]  r  V?  / 


2_rr  R-  "r  i  *5  t 


£-p  TTQi  'i t-  j 


-t,-n-e,v.  •v[(^>)aV]  = 


-  Rt  'V  .  \l~ 


I  ■  y 


pnro,vv\7[(ey^<p  ]  ---«'}'*v’') 


Applying  these  results  to  eqs.  (18)  the  result  is: 


C™ -pT.)  •*-  '7  .(-n  rpT,  Vj  +  a-TT ■ ^ CpT3  -  RT,  *•  V,  .^TT)  r 

r  aTT^.cpTa[»  -CA3MCPly 

vv7.(TrcpT,vA-  a-n-ijL  V,-  *t4(S£  +  v* . *tt)  * 

-  2.TT  cra  cpTx  [\|-)  -  \J  # 


(19) 


Finally  adding  the  four  equations  (14)  and  (19 )  and  integrating  over  the 
domain  the  result  is: 
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^-[TrLcp-r,  +iv,1  *  L  %,*■)]  =  (ftr,  V^qrr) 


Eq.  (20)  is  the  finite  difference  analogue  of  the  total  energy  equation. 

In  order  to  maintain  kinetic  energy  conservation  no  matter 
what  the  values  of  t  crx  and  v  are,  we  must  require: 

°-ri  -  ~  , 

°*Ti  -  <*> i  f 


cf 


t-  r  /  <i  \  ^ip 
■*-L  v  ^ 


-(i-V 

V2-/ 


d>  -6, 


(21) 


Eqs.  (21)  represent  the  finite  difference  form  of  the  hydrostatic  equation. 
Turning  to  eqs.  (16 )  and  writing  them  in  the  form: 


■  v 

+  V,* 

9©,' 

> 

i + 

©. 

i — i 

ifi 

-r-  I 

"Lot 

>+ 

rOTT 

^  <7  •  (-vr  V,)]  +  2-n-  <TX  Ox  r.  O, 
-*■  V7«(_tt  vO]  ~  a-TT^a.0z  “  °; 


(22) 


and  using  (9  ),  we  obtain. 


“30, 
o  r 


•■y-Oo. 


x5aC°x“©() 


J-  Va  ■  V>  -+  -  G.)  =-  O 


(23) 
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and  solving  eqs.  (21)  and  (27)  we  obtain: 


■b  - 

J_  ’  c  O.-vc1-':  . 

2-  \-v  '  ‘ 

--  C  *  *  £ )  '2.1  , 

cb  - 

-•  ;’ v*  ••  P.)  T. 

..V  ( (X  -V  0-j  T^] 

(28) 

X-  L 

<i>,  - 

1 3> 

k  j>"*5  T- 

-VC^-P*)  V\n 

We  are  now  in  a 


position  to  sum  up  the  equations  of  the  model  under  adiabatic 


and  frictionless  conditions: 


.|_(rru,j  •+-  <7  •Cttw.'V,)  Tr(.u,*u3)  <rx~^  O'''1  *  =  c/  u 


rrua)  +  -s*a  ~  17  8  "43)  ^  ~°i  1 


L(_TT  V',)  +  V‘Ctrv,  v,)  "**  *<T  (v\  *  Vs>fj  -l-  tt  |  W,  -V-  J v ~C/  ~ 


fTTV  '  v/^N  -TT(v,  +  V  HaV4.0°c5  tC/ 


- — -  U  VI 

ot  * 


^T,j  +  v.t^.v,)  -S£(3$^'<M  +  a1r4Ci)w“Ti  =  c/ 


Equations  I  -  IX  and  the  3  equations  (28)  thus  complete  the  system  of 
the  model  equations. 


FRICTIONAL  TERMS,  EDDY  DIFFUSION  TERMS  AND  D1ABATIC  HEATING 


.-is  the  vertical  diffusion  terms  and  the  surface  friction  were  es¬ 
sentially  modelled  in  the  same  way  as  sugg  i  by  Mintz  and  Arakawa,  their 
derivations  will  be  omitted  here.  The  only  difference  worth  mentioning  ,  is 
that  in  the  case  of  the  surface  stress  the  dependence  on  surface  density  was 
omitted,  i.  e.  it  was  made  to  depend  only  on  the  surface  wind  and  on  the  pres¬ 
sure.  The  surface  wind  was  extrapolated  linearly  from  the  winds  at  the  two 
reference  levels.  The  inclusion  of  changes  in  surface  density,  amounting  iu 
no  more  than  10%,  would  not  make  the  surface  ..'tress  any  more  accurate,  and 
would  take  up  more  computer  time.  The  horizontal  diffusion  terms  were  made 
proportional  to  the  Laplacian  of  the  diffused  quantity,  with  the  same  constant 
of  proportionality  for  momentum  sensible  heat  and  water  vapor. 

As  mentioned  above,  the  original  formulation  of  heat  sources 
and  sinks  was  completely  abandoned.  As  in  the  quasi-geostrophic  model  the 
only  heat  sources  taken  into  account  were  sensible  heat  flux  from  the  lower 
boundary  and  release  of  latent  heat  The  only  sinks  were  radiation  losses  to 
space  and  downward  flux  to  the  surface.  The  sensible  heat  flux  was  made  pro¬ 
portional  to  the  temperature  difference  between  air  and  ocean.  Assuming,  as 
we  did  before,  a  standard  atmosphere  lapse  rate  between  the  low>er  reference 
level  and  the  surface,  the  air  temperature  at  the  lower  boundary  became  pro¬ 
portional  to  T3.  The  heat  flux  enters  the  lowest  layer  only,  this  being  an 
additional  aspect  which  differs  from  the  previous  model.  The  radiation  losses 
were  calculated  in  the  same  way  and  were  made  equal  at  both  levels 
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ENERGY  TRANSFORMATIONS 

As  already  mentioned  above,  the  average  kinetic  energy  per 
unit  column  is  given  by  : 


14 


Multiplying  eqs.  I  -  IV  by  u  r  u  ? ,  v  ■'z  respectively  and  adding 
eq  VII  multiplied  by  ~  jG  i  "'i1  respectively,  we  obtain: 


-  c;  '  f  v,\  -  •  [3E*x  \ 

— -rr w, . ^7 4>4  -Tr^v^  ~°-T,  v,*Vtt  -RTj  v54?tt 


(29) 


Integrating  over  the  horizontal  extent  of  the  domain  : 


o  u.  _ 


~  C*.  •  V  <-  vv.  *  v  ^  LRT,  V,  v  R.T*  v3)  ,  V  rr 


(30) 


In  the  absence  of  frictional  losses  the  right  hand  side  represents 
the  transformation  of  potential  into  kinetic  energy : 

[?,  k]  =  -  5^  ^  +  C0.T,  Vt+  .  (31 ) 

Frictional  losses  are  divided  into  two  parts,  one  due  to  vertical  diffusion  of 
momentum  and  one  due  to  horizontal  diffusion.  These  losses  represent  the 
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energy  which  is  transformed  into  the  kinetic  energy  of  small  -  scale  eddies 
not  resolved  by  the  model.  For  these  two  transformations  we  obtain: 

L  r  ~  fg  OvT-  f, 

W,M  ?  rr  [  v  .  V  Vttv,)  t- 

C  being  the  eddy  diffusion  constant,  where  F^  and  Fg  are  given 
by  (Langlois  1968)  . 

F,  -  -  “^-sV175'  i 

Ft,  ~  B  ('V,  -'V-v}/tT1  ~  Hi  V^/tT  j 

B  and  D  are  constants. 

With  this  definition  of  F  and  F  the  final  result  becomes  . 

i.  u 

I>.p*l  =  +  £lv*|v*.vt/ 

with  given  by  : 

-  vV, )  . 

THE  HUMIDITY  FIELD 


(32) 


(33) 


As  in  the  previous  model  the  variable  chosen  to  represent  the 
humidity  field  was  the  specific  humidity,  q.  The  prognostic  equation  for 
the  specific  humidity  was  written  in  the  form: 

-V  V-(jrq_v)  +  ^  CO  fc)  ~  A  V V*?.)  ~Cy  (35) 


where  : 
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E  -  Specific  humidity  gain  through  evaporation, 

C  -  Specific  humidity  loss  through  precipitation^ 

A  -  Eddy  diffusion  coefficient  . 

The  formulation  of  the  evaporation  and  precipitation  mech¬ 
anism  adopted  in  the  previous  model  was  retained  also  for  the  present 
model,  except  that  the  amount  of  water  vapor  added  through  evaporation 
was  allowed  to  enter  the  lower  layer  only.  Thus  the  assumption  of  constant 
relative  humidity  with  height  was  dropped  and  no  instantaneous  adjustment 
of  the  humidity  distribution  through  evaporation  and  precipitation  was  re¬ 
quired.  As  seen  from  the  equation,  no  vertical  diffusion  of  water  vapor  was 
introduced  at  this  stage.  Adopting  the  same  vertical  differencing  as  for  tem¬ 
perature,  we  get  the  two  prognostic  equations  at  the  two  levels: 

d-  TTC, 

(,^r3  -  xtt  <r2  q2  4.  _  \  ^)sC  , 

q  was  estimated  as  a  weighted  average  of  q  and  q  .  Assuming  that  the 

£  1  u 

humidity  is  distributed  according  to  a  ^  law,  as  in  the  previous  model, 


(36) 


we  obtain: 
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with 


o  -  no  +•  ft  G* 

■  T I  •  C* 


A  +  <i  - 


Thus  • 


-> 


7?  C 


i7«H 


Experience  obtained  from  a  3-level  model  involving,  among  other 
aspects,  the  concentration  of  radon  as  a  time -dependent  variable,  showed 
that  the  poor  resolution  in  the  vertical,  with  the  assumed  time  step  may 
lead  to  negative  values  of  the  concentration.  This  was  due  to  the  fact  that 
on  the  average  the  concentration  at  the  intermediate  level  was  about  five 
times  larger  than  that  in  the  upper  level  and  about  three  times  smaller  than 
that  in  the  lower  level  when  linear  interpolation  was  used  for  q,  .  The 
vertical  advection  process  was  thus  overestimated  at  the  upper  level  and 
underestimated  at  the  lower  level.  To  overcome  this  difficulty  a  weighted 
average  was  used.  The  results  show  that  this  was  quite  adequate  and  the 
relative  humidity  never  went  below  about  30% 


BOUNDARY  AND  INITIAL  CONDITIONS 


As  mentioned  above,  the  equations  were  applied  to  the  same  domain 
as  in  the  previous  experiments.  Thus  we  require  cyclic  continuity  on  the 
eastern  boundaries.  On  the  northern  and  southern  boundaries  the  following 
boundary  conditions  were  chosen  : 


(1) 

1 SJ  *- 

-  c, 

J 

(2) 

%  C1TV)  ;v 

1 

--  '  2.oru)w  = 

(3) 

2-  (TM"' 

-D./  S  Nf 

•-  *o. 

Condition  (1 )  ensures  that  no  mass  leaves  the  domain  through  advection 
Condition  (2  )  and  (3  )  ensure  that  f  ere  is  no  diffusion  of  momentum  or 


>i'  heat  through  the  lateral  boundaries. 


The  initial  conditions  were  generated  in  the  following  way: 

First  the  ocean  surface  temperature  was  assumed  to  be  constant 
with  longitude  and  to  vary  linearly  from  27°C  at  j  =  1  to  -  25°C  at  j  =  14, 
with  a  gradient  of  4°C  per  grid  interval.  Next  it  was  assumed  that  the  air 
surface  temperature  was  equal  to  the  ocean  temperature,  so  that  initially 
no  sensible  heat  entered  the  model  atmosphere  from  the  lower  boundary. 

The  surface  pressure  was  assumed  to  be  constant  everywhere,  and  equal  to 
the  standard  pressure  of  1013  mb.  Next  the  temperature  at  level  3  was 
calculated  from  the  surface  temperature  by  assuming  a  standard  atmosphere 
lapse  rate  between  the  surface  and  level  3.  The  upper  level  temperatures 
were  assumed  to  vary  linearly  from  232°A  at  j  =  1  to  219°A  at  j  =  14, 
with  a  gradient  of  1°  per  grid  interval.  From  the  temperatures  at  both  levels 
the  geopofentials  at  both  levels  were  calculated.  At  this  stage  a  perturbation 
was  superimposed  on  the  zonally  averaged  geopotentials  of  the  following 
form  ,  -  in  r  a  ::-F  .  ^  IE*.  ,  where  Y  and  X  are  the  lengths 

of  the  domain  in  the  y  and  x  direction  respectively.  The  amplitude  chosen 
was  the  same  as  in  the  previous  model.  Next  the  temperatures  were  corrected 
to  correspond  to  the  perturbed  geopotentials  Finally  the  winds  were  assumed 
to  be  initially  geostrophic  and  were  calculated  from  the  goopotentials. 

The  Arakawa  finite  difference  scheme  was  retained  in  all  the  equations. 
The  grid  interval  used  was  the  same  as  in  the  previous  model,  but  the  time 
interval  was  reduced  to  10  minutes  to  ensure  mathematical  stability  of  the 
integration.  The  time  derivative  was  approximated  by  the  Matsuno  scheme. 


21 


FINITE  DIFFERENCE  SHCEME  AND  BOUNDARY  CONDITIONS 


As  mentioned  above  the  same  domain  chosen  for  the  previous 
experiments  was  retained.  It  contained  15  grid  points  in  the  x  and  y 
directions  respectively.  The  staggered  grid  configuration  proposed  by 
Arakawa  was  used.  As  it  was  thought  advisable  to  write  the  finite  difference 
equations  in  a  formulation  which  would  allow  a  straightforward  application 
to  the  computer  program,  the  different  variables  at  the  four  different  points 
were  assigned  the  same  indices  (see  Figure  below). 
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and  eq  VIII  becomes: 


t 

,  r  1 

1 

1 

i  ;  7T  Lu‘.i  "  U‘.i 

~  UvH,j 

,  4.  .  1 

1 

-v  \J  -  \j  ■ 

'•]  M 

-  v-i*. 

4  'A. in  J  . 

(39) 

The  boundary  conditions  for  the  last  two  eqs.  become: 

'J-  ~  V;  m-  =  o 

for  every  i  and  k  ,  where  k  denotes  the  level. 

The  last  condition  ensures  conservation  of  mass  in  the  form: 

£  V;  •  -(nv)  =  O  . 

M 

Next  we  turn  to  the  divergence  terms  in  eqs.  I  -  VI.  They  are 
written  in  the  form: 


V'  X.  -  +•  Uivt(i-I  -V 


-|K 
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-,+vsi. 

■V 
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4-  vbi  +  l 
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-  Vf-., 
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tKi-*  4  -V  (40) 
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With  x  denoting  either  u  or  v. 

To  ensure  the  condition  (  V  •  (_tr  /  v)  =  0  the  following 

s 

formulation  is  used  near  the  boundaries: 


(-w‘,i  +ui+«,» +ut+i,0 

-Tr;.l(,Cui,,  +u1-.1|,)(aill  +u;. +  f  vU 

+  3  VUCUV  *  +  |  +  ) 

(ttuv);,*  =g^  (.«;,«  +u;*l(,r)(u;|ir 

~"i'-t|wvCus>r  vu;.,tn*)(ui,if  ^  U  v_,  -  g  t^(_  ^  S«v  v  ^  C-((n 

^  2  *•  Ul+l,^)  +  ^  Cui,»T  +UH,uv)]  .  (41) 


The  v  equations  are  treated  in  a  different  way  as  with  v  =  0  on  the 
boundaries,  the  equations  are  not  applied  there.  This  requires  a  different 
formulation  of  the  scheme  at  a  distance  of  one  grid  interval  from  the  boundary, 
in  the  form. 
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^  v)  ,a 


•  g  (.  ^  v*  *  ^  i'A-  ^  ^V3  *■  ^*1 ,1}  ,■  CV‘ Vv  *VitvO 


-^-!v.'^m"v  *  ^-1,11  +  W-  ;.|,uv  +  ^i-',i3)(vi. 


|ift-Av.‘-i  >a; 


+  I  V2>  +  *  Vsi4.~  V.,^)  a  Vt'  +  ,  |(^ 

~  i(UV3  4  ^V;.,il(V+V;.|(ii)(v;|llvvV;.v^)  J  ,  <42) 


The  application  of  this  scheme  ensures  conservation  of  momentum  in  the 
flux  terms.  Whether  it  also  conserves  energy  does  not  appear  to  be  certain, 
but  it  is  expected  that  even  if  it  does  not,  the  error  will  be  very  small,  par¬ 
ticularly  when  a  more  extensive  domain  with  more  grid  points  will  be  used. 

The  vertical  flux  terms  in  eqs.  I  -  IV  are  written  in  the  form: 

^  ’H  (.x’  ■*  x*>)  <r  2  ]  l  •  - 

'  i^ii  *  **))(  +  +  (43) 

with  X  denoting  either  u  or  v. 

The  coriolis  parameter  f  is  defined  at  the  same  points  as  u  and  v  so  that: 


0  -  4t-TT;,.1 ^SJ  . 


$v)si  ;  ',TT;-Il‘j +usiH  +  ^£-'.1-,)  Ki  • 


In  the  u  equation  on  the  boundaries  it  is  of  course  omitted,  by 
virtue  of  the  condition  v  =  0,  so  that  this  presents  no  problem  as  far  as 
conservation  of  energy  is  concerned. 
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Next  we  write: 


'  *Ki  +ttm-'  4>si-«  'k-u  'Tr:-i,i-i  ^  y 

^f):„  •= 


o 

°Y 


aas^1’'  k-i,j  " TT*-.i -i  ^ i  j-i  k-i.i-i),  (45) 


and  finally: 
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+■  taTM-'  (-*•!  *')  C^iii-'  «•  -',j-0]  j 

[Ctr-4.)  =  4rJCaTM  -4ij  ^Tt-„  -  4,-JC^  -trl(j.() 

4  CftTi-',j  -  +aTc.,-.,  -4>C-,,i-.)C^c-tti  #  (46) 

On  the  boundaries  we  write  in  the  u  equations,  using  the  condition  £  ^4)  ;  C 

'  t  / 

;,,«  'TT;-v* 

fl.,!  T 

'  -L— [<•“•' V'V‘K'1  +  ^-Vft  >t-tr:-v*3(47) 


mi 
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Next  we  write: 


ias 


Om  <Ki  +7Tsi-'  4>si-i  -'"Hi  4> i-i .a  ^  * -i , i -• )  y 


<^L,i  i  j-i 
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On  the  boundaries  we  write  in  the  u  equations,  using  the  condition  ^  (jt  -  q 
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These  formulations  ensure  that  no  momentum  is  generated  by 
the  diffusion  terms. 

The  vertical  diffusion  and  surface  friction  terms  present  no 
problem  as  they  only  involve  quantities  defined  at  the  same  grid  point. 

Now  we  turn  to  eqs.  V  and  VI: 

The  horizontal  flux  terms  are  written  in  the  form: 


-*■  -Vc,j  (Ti'.j  +  •  (53) 

With  V- ,  =.  V:  ,r  =  0 »  this  definition  is  valid  also  for  the  flux  terms  near 
the  boundaries  so  that  at  j  =  1  the  fourth  term  is  omitted,  and  at  j  =  14  the 
third  term  is  omitted. 


For  the  vertical  flux  terms  we  first  substitute  for  Tg  its  value 


from  eq  IX ,  so  that : 


=Si,j  (Tij  . 


The  third  terms  in  eqs.  V  and  VI  present  no  problem  and  they  are  written 
as: 


Q.T  OJTl 

7? 


(55) 


with  the  definition  of  given  by  (38). 

Special  care  must  be  taken  with  the  remaining  term,  ~  Tv  9tr, 
This  term  represents  part  of  the  energy  transformation  from  potential  energy 
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into  kinetic  energy.  Thus  if  total  energy  is  to  be  conserved,  it  should  have 
the  same  form  as  the  corresponding  term  in  the  kinetic  energy  equation, 
which  is  defined  at  a  different  grid  point.  To  overcome  this  difficulty  the 
total  energy  is  defined  at  the  points  where  T  is  defined.  The  kinetic 
energy  is  then  calculated  as  an  average  of  its  values  at  the  four  surrounding 
u  points. 

Thus  the  calculation  of  the  terms  kT>y.  ^tT  proceeds  in  the 
following  way:  first  it  is  defined  at  the  points  where  u  is  defined,  by  using 
the  definition  of  the  terms  and  [o.T in  eqs.  I  -  VI. 


Next  we  define: 


(>Tv  ^]C(i 


=  c.i  . 


(56) 


Finally  the  value  of  this  term  at  the  points  where  T  is  defined 
is  calculated  as  an  average  of  its  values  at  the  four  surrounding  u  and  v 
points.  This  procedure  presents  no  special  problems  near  the  boundaries  as 
^ r  v  =  0  on  the  boundaries  and  R-T u  has  already 

been  defined  there. 

It  remains  now  only  to  define  the  time  differencing  scheme. 

For  this  purpose  we  use  the  definition  of  (rru)  and  (_tt  v) 
given  by  (48)  and  (51)  and  also  define: 


l0 


(57) 


For  the  time  differencing  the  so-called  Matsuno  scheme  is  used  in  the 
following  way; 
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Given  the  equation  :  -  ^(<*1,  *1/  ■■■  ,**)  (58) 

W  I 

the  first  step  is  an  ordinary  forward  difference  scheme  of  the  form: 


+♦1 


-  X 


M  f  L< 


.*) 


(59  > 


and  the  second  step  is  : 
..++1 


*'  -.st  t  (*]"' ,  , ... , 


(60) 


The  two  equations  (28)  for  <f>'^  and  4>  V  are  calculated  at  each  point 
by: 


‘Ki  -■  ilu+«)r.v  *  (B.+ (4)^*3, 
4>'i 


(61) 


The  third  of  the  equations  (28)  for  has  already  been  used  to  transform 

the  last  terms  of  the  eqs.  I  -  IV,  so  that  cfs2-  does  not  appear  explicitely 
in  the  program. 
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Analysis  of  results 

Two  integrations  were  performed.  The  integration  not  involving 
release  of  latent  heat  will  be  referred  to  as  Experiment  I,  and  the  experiment 
in  which  water  vapor  is  carried  as  Experiment  II.  In  Experiment  I  the  model 
equations  were  integrated  for  a  period  of  37  days,  and  the  integration  was  dis¬ 
continued  when  a  quasi-steady  state  was  attained  and  no  new  information  could 
be  gained  by  prolonging  the  integration.  Consequently  emphasis  will  be  placed 
in  the  following  on  the  conditions  prevailing  on  the  last  days  of  the  integration. 

The  early  stages  of  development  were  characterized  by  a  rapid 
growth  of  the  amplitude  of  the  disturbance  at  both  levels,  similar  to  the  develop¬ 
ment  observed  in  the  quasi -geostrophic  model  (Hubs  1967)  referred  to  above. 

It  is  worth  while  mentioning  that  this  growth  stopped  on  the  third  day  and  a  de¬ 
crease  in  amplitude  was  observed  between  the  third  and  seventh  days  (Fig.  1). 
Afterwards  the  amplitudes  started  to  grow  again  until  they  reached  their  steady 
state  value  on  about  the  20th  day.  The  temporary  decrease  in  amplitude  was 
observed  in  all  the  experiments  with  the  quasi-geostrophic  model,  though  start¬ 
ing  at  a  later  date,  i.  e.  on  the  12th  day.  It  thus  appears  that  a  tendency  towards 
the  development  of  an  index  cycle  is  observed  in  all  experiments.  This  tendency, 
however,  is  damped,  in  most  cases,  as  the  integration  is  continued.  The  reason 
for  this  damping  is  not  clear  from  the  results.  It  might  be  due  to  the  very  crude 
formulation  of  the  surface  friction  and  vertical  diffusion  terms,  which  were  ex¬ 
pected  to  have  a  damping  effect.  On  the  other  hand  it  might  also  be  due  to  the 
finite  difference  formulation  of  the  model  equations,  as  some  difference  schemes 
are  known  to  have  marked  damping  effects. 

In  the  quasi-geostrophic  experiment  involving  water  vapor,  a 
marked  index  cycle  developed,  which  was  not  damped  until  the  end  of  the  integ¬ 
ration  on  the  67th  day.  Thus  it  seems  that  in  the  models  strong  differential 
heating  is  required  in  order  to  maintain  the  index  cycle  and  counteract  the 
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damping  effect,  whatever  its  origin  may  be.  In  the  present  case  the  values  of 
the  sensible  heat  flux  in  the  south  were  lower  than  in  the  previous  experiments 
and  thus  the  index  cycle  appears  to  have  been  damped  very  fast. 

In  addition,  the  present  model  contains  lateral  eddy  diffusion 
terms,  which  have  a  damping  effect  and  thus  it  was  expected  that  more  pro¬ 
nounced  differential  heating  would  be  required  in  order  to  maintain  the  index 
cycle  The  values  of  the  amplitudes  of  the  geopotential  waves  reached  towards 
the  end  of  the  integration  were  about  half  of  the  corresponding  values  obtained 
in  the  quasi-geostrophic  case.  This  result  was  expected  and  is  a  consequence 
of  the  summer  conditions  imposed  on  the  model. 

The  mean  zonal  wind  distribution  (Fig.  2)  is  characterized  by 
a  maximum  westerly  flow  at  the  upper  level  at  the  central  latitude  and  easter¬ 
lies  at  lower  levels  near  the  boundaries.  However  the  intensity  of  the  jet  was 
weaker  than  in  the  quasi-geostrophic  model.  The  maximum  value  obtained  was 
about  25  m  sec  1  which  is  characteristic  for  summer  conditions  in  the  real 
atmosphere.  Summer  conditions  were  obtained  in  the  model  by  reducing  the 
surface  temperature  meridional  gradients  from  4.  5°  per  grid  interval,  in  the 
quasi-geostrophic  case,  to  4°  per  grid  interval  in  the  present  case.  This  of 
course  leads  to  a  much  weaker  thermal  wind  and  thus  to  a  weaker  jet.  Thus, 
with  a  proper  choice  of  the  surface  temperature  gradient  it  is  possible,  even 
with  the  crude  heating  function  assumed  in  the  model,  to  simulate  either  winter 
or  summer  conditions  at  least  as  far  as  temperature  and  wind  conditions  are  con¬ 
cerned.  One  undesirable  result  was  the  appearance  of  stronger  easterlies  near 
the  northern  boundary  as  compared  with  the  belt  of  easterlies  in  the  south.  This 
seems  to  be  a  result  of  thto  assumption  of  a  constant  coriolis  parameter,  whereas 
in  the  quasi-geostrophic  case  this  strong  restriction  is  partly  offset  by  assuming 
a  non  vanishing  ^  term  in  the  vorticity  equation.  This  restriction  is  not  an  in¬ 
herent  aspect  of  the  present  model,  as  a  variable  coriolis  parameter  can  easily 
be  introduced.  In  that  respect  the  present  results  are  only  of  a  preliminary  nature 
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and  a  variable  coriolis  parameter  will  be  introduced  In  future  experiments. 

A  schematic  presentation  of  the  mean  meridional  circulation  is 
shown  in  Fig.  7.  The  qualitative  picture  agrees  quite  well  with  the  one  observed 
in  nature,  with  two  direct  cells  in  the  northern  and  southern  parts  of  the  region 
and  an  indirect  cell  in  the  middle.  The  southern  Hadley  cell  is  of  course  too 
small  compared  with  the  corresponding  Hadley  cell  in  the  real  atmosphere, 
probably  due  to  the  presence  of  a  "wall"  as  the  southern  boundary. 

Its  northern  edge  is  situated  at  the  right  location  at  about  30°  N 
and  the  corresponding  meridional  and  vertical  wind  speeds  have  the  correct 
magnitude,  appropriate  for  summer  conditions  (Mintz  and  Lang  1955).  The  central 
indirect  cell  is  also  too  small  in  its  latitudinal  extent  and  its  northern  edge  is 
situated  too  far  to  the  south,  thus  contributing  to  an  incerease  in  the  size  of  the 
northern  direct  cell.  This  again  seems  to  be  a  result  of  the  assumption  of  a  con¬ 
stant  coriolis  parameter  together  with  the  restriction  of  the  integration  to  a  small 
domain.  Thus  it  appears  that,  with  the  shortcomings  mentioned  above,  the  summer¬ 
time  mean  meridional  circulation  is  reasonably  well  simulated  by  the  model. 

The  mean  zonal  surface  pressure  is  presented  in  Fig.  8.  Again 
as  noted  in  the  case  of  some  other  variables,  the  qualitative  distribution  is  quite 
realistic.  Thus  the  tropical  low  pressure  belt,  the  subtropical  high  pressure  belt, 
the  low  of  the  mid  latitudes  and  the  polar  highs  are  simulated  quite  well  in  the 
correct  latitudes.  However,  the  pressure  difference  between  the  low  and  high 
pressure  belts  is  too  large,  if  compared  to  conditions  prevailing  in  the  northern  hemi¬ 
sphere  and  too  small  if  compared  to  those  in  the  southern  hemisphere.  In  the 
northern  hemisphere  of  the  real  atmosphere  the  pressure  distribution  is  very 
much  affected  by  the  ocean  continent  contrasts,  as  most  of  the  land  masses  are 
concentrated  there.  This  effect  is  smoothed  out  when  zonal  averages  are  computed, 
making  the  mean  meridional  pressure  gradients  much  smaller  than  those  of  the 
southern  hemisphere. 
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It  should  be  noted  here  that  the  initial  surface  pressure  was 
chosen  to  correspond  to  the  mean  pressure  at  sea  level  observed  in  the  atmo¬ 
sphere.  With  the  exclusion  of  the  mountians  in  the  model,  this  means  that  the 
model  contains  more  mass  than  the  real  atmosphere.  This  is  expected  to  lead 
to  exaggerated  values  for  the  fluxes  of  momentum,  heat  and  water  vapor.  One 
way  of  overcoming  this  difficulty  is  by  assuming  an  initial  surface  pressure 
corresponding  to  the  real  mass  of  the  atmosphere  as  done  by  Smagorinsky 
(1965).  However,  this  was  not  done  in  the  present  experiments. 

The  variation  with  time  of  the  total  kinetic  energy  is  given  in 
Fig.  5.  As  with  the  amplitude  of  the  geopotentials,  traces  of  an  index  cycle 
are  observed  at  the  early  stages  of  the  integration.  Starting  on  the  15th  day 
a  constant  decrease  of  kinetic  energy  is  observed,  lasting  until  the  end  of  the 
integration.  As  mentioned  above,  an  instantaneous  heat  balance  was  imposed 
on  the  model,  so  that  there  was  no  net  heating  or  cooling  of  the  whole  atmo¬ 
sphere.  The  energy  losses  due  to  surface  friction  and  vertical  anrl  lateral  dif¬ 
fusion  could  not  be  replenished  from  any  other  sources.  It  was  thus  expected 
that  the  total  energy  in  the  model  would  decrease  with  time,  thus  limiting  the 

validity  of  the  prediction  for  a  long  period.  The  observed  decrease  in  kinetic 

-2 

energy  from  the  15th  to  the  37th  day  was  about  100  Kj.  m  which  was  about 
4%  of  the  total  kinetic  energy.  This  again  motivated  the  decision  to  discontinue 

— the  integration  before  this  effect  became  too  serious  The  total  kinetic  energy 

_2 

at  the  end  ol  the  integration  was  about  2200  Kj  m  This  value  is  somewhat 
higher  than  the  average  values  obtained  by  Oort  (1964)  but  lies  within  the  limits 
of  error  quoted  in  his  paper  This  high  value  could  be  explained  by  the  fact  that 
we  had  too  much  mass  in  the  model  and  these  high  values  were  obtained  even 
though  the  wind  velocities  were  quite  realistic 

Evidence  of  the  damping  of  the  index  cycle  is  most  clearly  obser¬ 
ved  in  the  values  of  the  transformation  of  potential  energy  into  kinetic  energy, 
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ay  a  function  of  time  (Fig.  4).  A  series  of  oscillations  is  observed  with  a  period 

of  about  7  days,  and  a  decreasing  amplitude  as  the  integration  preoceeds,  until 

total  damping  of  the  oscillation  takes  place  at  the  end  of  the  experiment.  The 

*2  -1 

value  of  the  transformation  on  the  last  day  was  about  130  Kj.  m  day  , 

which  is  smaller  by  a  factor  of  two  as  compared  with  the  value  presented  by  Oort 
(1964).  The  kinetic  energy  losses  due  to  lateral  eddy  diffusion  were  rather  low, 
and  never  exceeded  15  Kj.  m  ^  day  \ 

A  comparison  of  the  kinetic  energy  changes  computed  from  values 

on  two  consecutive  days,  with  the  energy  changes  computed  from  the  energy  trans- 
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formations,  reveals  a  difference  of  the  order  of  10  Kj.  m  day  .  With  the  care 
taken  in  the  model  to  ensure  conservation  of  energy,  it  was  concluded  that  this 
difference  could  be  attributed  only  to  the  deficiency  (from  the  point  of  view  of 
energy  conservation)  of  the  time  scheme  employed  in  the  model.  Inspecting  the 
daily  energy  changes  in  the  model  we  observe  that  this  fictitious  energy  loss  ac¬ 
counts  for  about  10%  of  the  change.  The  cumulative  effects  might  thus  be  quite 

serious,  in  the  long  run.  Combining  the  kinetic  energy  losses  due  to  all  three 
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effects,  we  obtain  a  value  of  155  Kj  m  day  for  the  total  kinetic  energy  sink, 

-2  -1 

as  compared  to  a  value  of  about  200  Kj.  m  day  obtained  by  Oort. 

We  thus  conclude  that  the  energy  cycle  obtained  in  the  model  was 
somewhat  too  weak,  a  fact  which  is  consistent  with  the  damping  effect  observed. 
The  dominant  kinetic  energy  losses  were  caused  by  the  surface  friction  and  the 
vertical  eddy  diffusion.  The  losses  due  to  lateral  diffusion  were  an  order  of  mag¬ 
nitude  smaller.  This  differs  markedly  from  the  results  obtained  by  Smagorinsky 
(1965)  where  the  lateral  diffusion  accounted  for  a  loss  of  about  twice  as  much 
as  the  vertical  diffusion.  Oorts  results  do  not  give  seperate  values  for  these 
two  effects  for  the  atmosphere  so  it  could  not  be  decided  which  of  the  results  is 
closer  to  reality,  though  from  the  results  given  by  Phillips  (1956) ,  it  appears 
that  both  effects  should  have  the  same  order  of  magnitude. 
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Next  we  turn  our  attention  to  the  heat  balance  components  (Fig.  6) . 

_2 

The  net  heating  had  a  marked  maximum  at  j  =  3  with  a  value  of  about  0. 1  Kj  m 

-1  -2  -1 
sec  and  a  minimum  at  j  -  12  with  a  value  of  about -0.  08  Kj  m  sec  .  Qual¬ 
itatively  this  agrees  quite  well  with  the  distribution  given  by  Rasool  and  Prabhakara 
(1965),  for  the  summer  season  in  both  hemispheres.  The  maximum  heating  in 

August  in  the  northern  hemisphere  is  located  at  about  15°N  with  a  value  of  about 
-2  -1 

0. 15  Kjm  sec  ;  Our  maximum  thus  lies  too  far  to  the  north  and  the  values  are 
somewhat  too  low.  The  minimum  on  the  other  hand  is  located  too  far  south.  This 
seems  to  be  a  result  of  the  choice  of  a  too  small  domain  and  the  effect  of  the  "walls". 
The  gradients  of  differential  heating  obtained  are  thus  quite  realistic,  which  is 
surprising  with  the  simple  heating  function  employed  in  the  model.  The  southern 
"wall”  however  seems  to  arrest  the  heating  processes  there,  and  thus  in  the  south¬ 
ern  part  of  the  region  the  gradients  of  heating  are  exaggerated  as  there  is  almost 
no  heating  near  the  boundary.  As  a  result  of  this  the  total  excess  of  heating  in  the 
southern  region  over  the  northern  region  is  much  too  low  and  consequently  the 
heat  transport  to  the  north  is  also  too  small.  The  present  experiment  carries  no 
water  vapor,  which  is  expected  to  account  for  a  large  part  of  the  heating  in  the 
south  and  thus  we  conclude  that  the  inclusion  of  this  process  is  essential  for  a 
correct  distribution  of  heat  sources  and  sinks  in  the  model. 

The  radiational  losses  obtained  were  rather  low  compared  with 
values  obtained  by  Astling  and  Horn  (1964),  from  Tiros  II  measurements.  The 
satellite  data  contained  of  course  the  total  outgoing  long  wave  radiation  and  thus 
included  also  that  part  coming  from  the  lower  surface  through  the  atmospheric 
"window".  Nevertheless  the  low  values  obtained  in  the  model  have  a  marked  effect 
on  the  distribution  of  the  heating  with  height,  making  the  lower  level  heating  at  least 
an  order  of  magnitude  stronger  than  the  upper  level  heating. This  is  in  contrast  to 
the  previous  model  where  the  heating  was  distributed  evenly  in  the  column.  The 
situation  in  the  atmosphere  probably  lies  between  these  two  extremes. 
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Inspecting  the  poleward  flux  of  zonal  momentum  at  both  levels 

the  following  distribution  is  observed:  strong  poleward  fluxes  appear  at  the  upper 

2  -2 

level  in  the  southern  part  of  the  region,  with  a  maximum  value  of  8.  9  m  sec 

2 

at  j  =  6 .  Maximum  southward  fluxes  at  the  upper  level  with  a  value  of  13.  9  m 
-2 

sec  were  observed  at  j  =  10.  These  values  are  low,  compared  to  the  summer 
values  obtained  by  Holopainen  (1967).  The  northern  maximum  is  situated  too  far 
south  as  a  result  of  the  southward  shift  of  the  edge  of  the  central  indirect  merid¬ 
ional  cell,  commented  on  above. 

The  mean  zonal  poleward  heat  transport,  on  the  other  hand,  shows 
a  better  qualitative  agreement  with  the  situation  in  the  real  atmosphere,  as  pre¬ 
sented  by  Rasool  and  Prabhakara  (1965).  Thus  a  small  region  of  southward  fluxes 
appears  near  the  southern  boundary,  at  j  =  1  and  2.  Further  north,  northward 
fluxes  prevail,  reaching  a  maximum  value  at  j  =  7  and  then  decreasing  until 
southward  fluxes  are  encountered  again  at  j  =  15.  The  values  presented  in  the 
paper  referred  to  above,  were  computed  for  latitude  circles,  so  that  with  the 
different  length  of  each  latitude  circle  in  the  atmosphere,  as  compared  to  our 
planar  zonal  strip,  actual  comparison  of  values  is  difficult.  As  the  upper  level 
temperatures  obtained  in  the  model  are  somewhat  too  high  and  the  meridional  wind 
components  have  the  right  magnitude,  it  appears  that  the  values  obtained  for  the 
poleward  flux  of  heat  are  somewhat  too  high. 

Experiment  II 

The  second  experiment,  with  the  model  including  a  water  vapor 
cycle  and  release  of  latent  heat,  was  continued  for  a  period  of  32  days.  As  before, 
a  quasi-steady  state  was  attained  and  further  integration  seemed  pointless.  As  it 
is  desired  to  assess  the  influence  of  the  release  of  latent  heat  on  the  development  of 
the  circulation  patterns,  the  water-vapor  budget  will  be  discussed  first  and  the 
discussion  of  synoptic  manifestation  will  follow  later. 
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The  choice  of  a  smaller  meridional  gradient  of  the  sea  surface 
temperatures  led  of  course  to  lower  air  temperatures  in  the  southern  part 
of  the  region,  as  compared  with  the  previous  model.  Consequently  the  water 
vapor  in  the  south,  where  most  of  the  humidity  is  expected  to  be  concentrated, 
was  much  lower  in  the  present  case.  For  example,  the  amount  of  preceipit- 
able  water  per  unit  column  at  j  =  1  was  about  37  mm  in  the  pres  nt  case,  as 
compared  to  47  mm  at  the  same  latitude  in  the  previous  model.  This,  of 
course,  had  a  marked  effect  on  the  precipitation  and  evaporation  in  the  south 
The  distribution  with  latitude  of  the  mean  zonal  relative  humidity,  on  the  other 
hand,  was  similar  to  that  obtained  in  the  previous  experiments.  (Fig.  9).  Thus 
a  maximum  in  the  relative  humidity  appeared  at  the  southern  boundary,  a 
minimum  at  j  =  5  and  a  second  maximum  close  to  the  northern  boundary  of  the 
region. 


As  in  the  previous  model,  the  meridional  gradient  of  relative 
humidity  was  exaggerated,  as  the  relative  humidities  in  the  south  and  in  the 
north  were  too  high,  and  the  relative  humidity  in  the  dry  belt  was  too  low, 
compared  to  the  values  presented  by  Kaurwitz  and  Austin  (1944) 

Thus  it  seems  that  the  shortcomings  of  the  humidity  distribution 
in  the  present  experiment  result  from  too  low  temperatures  in  the  south.  This 
can  easily  be  corrected  in  future  experiments  designed  to  simulate  summer 
conditions  by  assuming  a  higher  average  sea  surface  temperature,  while  re¬ 
taining  the  same  meridional  gradients.  This  will  not  change  the  initial  gradients 
of  tho  geopotentials,  and  the  same  initial  wind  distribution  will  also  be  retained, 
while,  on  the  other  hand  higher  values  for  the  specific  humidity  will  be  obtained 
throughout  the  region. 

The  inspection  of  the  zonally  averaged  distribution  of  the  precip¬ 
itation  and  evaporation  amounts  reveals  the  same  shortcomings.  While  the 
qualitative  distribtution  is  essentially  correct,  with  two  precipitation  maxima, 
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one  at  the  southern  wall  and  the  other  somewhat  to  the  north  of  the  central 
latitude  and  an  evaporation  maximum  at  j  =  3  -  the  total  amounts  are  much 
too  low  Moreover,  the  southern  maximum  in  precipitation  is  lower  than 
the  mid-latitude  maximum,  contrary  to  the  situation  in  the  atmosphere,  where 
the  tropical  maximum  is  dominant,  especially  during  summer.  The  same 
effect  was  observed  also  in  the  previous  model  and  is  due  to  the  presence  of 
the  southern  boundary  wall  and  to  the  exclusion  of  the  tropical  regions. 

The  excess  of  precipitation  over  evaporation,  which  determines 
in  the  real  atmosphere  the  net  release  of  latent  heat,  has  in  the  model  a  max.- 
imum  positive  value  at  about  j  =  11  and  a  maximum  negative  value  at  j  =  3. 

In  the  model,  however,  this  variable  does  not  determine  the  net  release  of 
latent  heat,  as  it  was  assumed  that  the  amount  of  heat  needed  for  evaporation 
is  drawn  from  the  ocean,  which  is  assumed  to  have  an  infinite  heat  capacity, 
and  not  from  the  atmosphere.  Therefore  the  net  release  of  latent  heat  is  always 
positive  in  the  model.  This  agrees  with  the  results  obtained  by  Rasool  and 
Prabhakara  (1965)  for  the  summer  season  in  both  hemispheres.  In  the  atmo¬ 
sphere  this  means  that  there  exists  an  appreciable  transport  of  latent  heat 
between  the  two  hemispheres,  an  effect  which  cannot  be  in  a  model  simulating 
only  a  part  of  one  hemisphere. 

In  conclusion  it  may  be  stated  that  though  the  qualitative  distri¬ 
bution  of  the  elements  of  the  water  vapor  cycle  was  simulated  quite  well,  the 
various  mechmisms  were  underestimated.  Consequently,  it  should  be  ex¬ 
pected  that  the  total  effect  of  the  release  of  latent  heat  on  the  circulation  pat¬ 
tern  in  the  model  would  not  be  very  marked. 

Bearing  this  result  in  mind,  we  now  turn  our  attention  to  the 
synoptic  manife  '-tions. 

The  mean  zonal  flow  onthe  last  day  of  the  integration  is  pre¬ 
sented  in  Fig.  3.  The  only  difference  from  the  result  of  experiment  I  is 


39 


the  appearance  of  stronger  easterlies  in  the  north  and  weaker  (almost  vanish¬ 
ing)  easterlies  in  the  south.  This  undesirable  feature  is  even  more  pronounced 
in  the  present  experiment.  It  is  thus  concluded  that  a  variable  coriolis  para¬ 
meter  is  essential  for  obtaining  more  realistic  circulation  patterns,  even  in 
models  not  containing  tropical  regions.  This  could  be  seen  even  in  some  ex¬ 
periments  with  the  quasi-geostrophic  model  where  the  (*,  term  has  been  drop¬ 
ped  (Huss  1964). 

The  rate  of  growth  of  the  amplitudes  of  the  geopotential  distur¬ 
bance  (Fig.  1),  on  the  other  hand,  became  much  more  rapid  when  precipitation 
Btrated  on  the  fifth  day.  Thus  a  value  of  about  130  g.  p.  m.  was  reached  at  the 
upper  level  on  the  tenth  day,  as  compared  to  a  value  of  100  g.  p.  m.  on  the  same 
day  in  the  previous  experiment  The  amplitude  remained  higher  in  the  present 
experiment  during  the  whole  period  of  integration  -  although  at  the  end  of  the 
period  the  difference  was  not  very  marked  and  amounted  only  to  about  10  g.  p.  m. 
at  the  upper  level.  The  damping  of  the  index  cycle  also  proceeded  at  a  some¬ 
what  slower  rate,  supporting  our  contention  that  stronger  heating  might  sustain 
the  index  cycle. 

Evidence  of  a  slower  damping  is  observed  also  in  the  changes 
with  time  of  the  total  kinetic  energy  (Fig.  5)  At  the  early  stages  a  marked 
increase  in  the  kinetic  energy  was  observed  right  after  precipitation  started, 
and  at  least  three  appreciable  oscillations  were  observed  before  total  damping 
took  place  and  the  energy  started  to  decrease  continuously,  as  it  did  in  the  pre¬ 
vious  experiment.  The  period  of  oscillation  remained  the  same,  i.  e.  about  7 
days. 

The  energy  cycle  on  the  other  hand  was  much  more  pronounced 
than  in  experiment  I.  Thus  a  very  rapid  increase  in  both  the  trans¬ 
formation  and  the  transformation  was  observed  after  the  start  of  precip- 
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itation.  The^p,^  transformation  reached  a  value  of  about  140  Kj  m  day  on 
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the  third  day  as  c  <m pared  to  a  value  of  about  90  Kj  m  day  on  the  same  day 
in  experiment  I.  The  LP.dJ  transformation  oscillated  between  high  and  low 
values  until  the  end  of  the  integration,  with  a  larger  amplitude  than  in  experi¬ 
ment  I  However,  the  oscillations  were  still  very  weak  and  are  hardly  char¬ 
acteristic  of  the  index  cycle  oscillations  observed  in  nature.  The  average  value 

-2  -1 

of  the  transformation  of  potential  into  kinetic  energy  was  about  160  Kj  m  day 
This  value  is  still  somewhat  too  low,  compared  with  the  values  presented  by 
Oort  (1964).  It  might  thus  be  concluded  that  a  more  realistic  distribution  of  the 
water  vapor  cycle  elements  would  have  an  even  greater  effect  toward  increas¬ 
ing  the  energy  cycle  elements  and  would  bring  them  closer  to  the  observed  values 

An  inspection  of  the  heat  balance  components  (Fig.  o)  reveals  that 
the  net  heating  in  the  southern  part  of  the  region  was  almost  the  same  as  in 
experiment  I,  and  the  only  effect  of  the  release  of  latent  heat  there  was  to  some¬ 
what  reduce  the  gradients  of  differential  heating.  In  the  northern  legions  where 
the  precipitation  maximum  was  located,  the  effect  was  more  pronounced  The 
release  of  latent  heat  in  the  north  resulted  in  the  direct,  heating  of  the  lower  level 
As  a  result  the  temperature  difference  between  the  atmosphere  and  ocean  suY- 
face  increased,  leading  to  higher  values  of  cooling  due  to  the  sensible  heat  /lux 
from  the  atmosphere  into  the  ocean.  Thus  the  net  result  of  the  release  of  latent 

heat  in  the  north  was  to  increase  the  net  cooling,  instead  of  reducing  it.  The 
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maximum  value  of  net  cooling  was  about  0. 13  Kj  m  sec  as  compared  to  a 

-2  -1 

value  of  about  0.  08  Kj  m  sec  obtained  in  experiment  I.  As  the  amount  of 
heating  in  the  south  was  the  same  in  both  experiments  we  conclude  that  the  dif¬ 
ferential  heating  was  larger  in  the  present  experiment,  a  fact  which  accounts 
for  the  higher  amplitudes  observed  and  the  slower  damping  of  the  index  cycle. 

The  surface  pressure  distribution  was  also  qualitatively  the  same 
as  in  the  first  experiment,  with  the  low  and  high  pressure  belts  located  at  the 
same  latitudes.  The  meridional  gradients  of  the  surface  pressure,  however. 
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were  somewhat  larger  than  in  the  first  experiment.  The  distribution  was  tl 
closer  to  the  situation  in  the  southern  hemisphere,  so  that  as  far  as  surfac 
pressure  is  concerned  the  inclusion  of  the  hydrologic  cycle  leads  to  an 


improvement. 
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CONCLUSION 

Reviewing  the  results  of  the  last  two  experiments,  as  comparer! 
with  the  results  of  the  experiments  with  the  quasi-geostrophic  model,  the 
following  conclusions  may  be  drawn1 

The  mam  qualitative  features  of  the  general  circulation  were 
simulated  quite  well  in  all  experiments.  These  features  include  the  meridion¬ 
al  circulation  consisting  of  three  cells,  the  central  jet  stream,  and  tl.  i  more 
or  less  adequately  simulated  mean  zonal  flow.  In  fact,  these  features  were 
obtained  in  some  models  even  without  the  introduction  of  diabatic  heating.  On 
the  other  hand,  it  has  to  be  remembered  that  the  initial  distributions  of  the 
geopotential  and  temperature  fields  were  already  quite  close  to  tl  v  observed 
distributions,  at  least  as  far  as  zonal  averages  are  concerned.  It  thus  seems 
that'  the  maintenance  of  the  general  circulation,  at  least  for  relatively  short 
periods,  is  not  very  much  influenced  by  the  heating  function  employed.  On 
the  other  hand  the  generation  of  the  features  mentioned  above,  starting  from 
an  isothermal  atmosphere  at  rest,  would  depend  very  much  on  the  heating. 

We  thus  conclude  that  even  a  simple  model  with  a  simple  heating  function 
might  be  useful  for  certain  purposes.  For  example,  such  a  model  might  be 
used  to  get  a  rough  estimate  for  the  concentration  with  latitude  of  radioactive 
fallout  from  an  atomic  explosion.  For  predictions  of  changes  in  the  general 
circulation,  such  as  the  index  cycle  or  the  seasonal  variations,  it  appears 
that  the  inclusion  of  a  realistic  heating  function  is  essential.  In  that  respect 
it  seems  that  the  inclusion  of  the  hydrologic  cycle  goes  quite  a  way  towards 
making  the  heating  function  more  realistic.  In  addition  it  provides  useful 
information  about  large  scale  precipitation  areas. . 

One  essential  factor  left  out  in  the  experiment  was  the  ocean- 
continent  contrast.  This  factor  must  be  included  in  actual  forecasts  of  changes 
in  the  general  circulation  as  the  zonal  distributions  of  the  various  field  vari¬ 
ables  affecting  the  weather  were  of  course  not  simulated  at  all  by  our  models. 
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AREA  AVERAGED  KINETIC  ENERGY  TRANSFORMATIONS  IN  kjmJday' 
RUN  WITHOUT  WATER  VAPOR 


AREA  AVERAGED  KINETIC  ENERGY  TRANSFORMATIONS  IN  kj  m^day'’ 
RUN  WITH  WATER  VAPOR 

UO. - r - , - 1 - i - t - r 


Fig.  4.  Mean  kinetic  energy  transformations.  P,  K  ]  =  transformation 
of  potential  into  kinetic  energy,  K,  F  =  transformation  of 
kinetic  energy  into  surface  friction.  ^_K,  A1  =  transformation 
of  kinetic  energy  into  eddy  diffusion.  Experiments  I  and  n 


AREA  AVERAGED  KINETIC  ENERGY  PER  UNIT  COLUMN  IN  kj  nv2 
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Fig.  5.  Area  averaged  kinetic  energy.  Experiments  I  and  II 
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ZONALLY  AVERAGED  HEAT  BALANCE  COMPONENTS 
ON  36,h  DAY  IN  kj  rrWsec*' x10r 
RUN  WITHOUT  WATER  VAPOR 


ZONALLY  AVERAGED  HEAT  BALANCE  COMPONENTS 
ON  32nd  DAY  IN  kj  nr?  sec1  x  102 
RUN  WITH  WATER  VAPOR 


Fig.  6.  Heat  balance  components.  Experiments  I  and  H 


Fig.  8.  Mean  zonal  surface  pressure.  Experiments  I  and  II 


MEAN  ZONAL  RELATIVE  HUMIDITY  IN  •/.  ON  32Nn  DAY 
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INCLUSION  OP  LATERAL  VISCOSITY 
AND  THE  APPLICATION  OP  A  MATSUNO  SCHEME 
IN  A  TWO  LEVEL  MODEL  OP  THE 
GENERAL  CIRCULATION  OP  THE  ATMOSPHERE 

b  y 

A.  HUSS  and  M.  NAVON 

Abstract 

Equations  of  a  two  level  qua si -hydrostatic  and  quasi- 
geostrophic  model,  have  been  numerically  integrated  in  order 
to  study  the  evolution  of  general  circulation  patterns. 

A  lateral  friction  term  was  included  in  a  model  whose 
main  features  are  described  by  Hues  (1967). 

The  eddy  viscosity  coefficient  was  assumed  to  be  cons¬ 
tant.  Two  different  integrations  were  performed.  The  first  one 
involving  an  integration  of  k2  days  was  based  on  a  model  which 
did  not  include  water  vapor. 

In  the  second  experiment  which  was  integrated  for  65  days 
a  humidity  distribution  was  introduced  and  the  lateral  friction 
term  was  also  included. 

The  calculated  distributions  of  the  field  variables  were 
compared  with  the  results  obtained  in  the  1967  report  mentioned 
above  and  with  actual  distributions  based  on  real  data.  In  both 
the  experiments  a  considerable  smoothing  of  the  var:'  us  field 
variables  was  obtained  and  considerable  differences  were  observed, 
compared  to  the  previous  experiments. 
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In  the  second  experiment  the  smoothing  destroyed  the 
pronounced  index  cycle  which  appeared  in  a  previous  exper iment. - 

In  the  third  experiment,  the  Matsuno  time  differencing 
.cherne  w  ■  3  introduced  instead  of  the  leap-frog  scheme,  while 
no  lateral  friction  was  included. 

The  numerical  integration  was  carried  for  62  days.  A 
smoothing,  similar  to  that  obtained  in  the  first  two  experiment 
was  observed,  but  with  a  more  realistic  distribution  of  the 
field  variables. 
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INTRODUCTION 

a )  The  Lateral  Friction  Term 

De  i' In  in  £  an  eddy  as  a  parcel  of  atmospheric  air  capable 
of  retaining  some  measure  of  individuality  and  moving  as  a 
separate  entity  in  the  surrounding  fluid  -  we  shall  take  the 
eddy  instead  of  the  molecule  as  the  fundamental  part  -  so 
that  entities  3uch  as  eddy-viscosity  or  eddy-conductivity 
will  be  the  macroscopic  counterpart  of  the  molecular  viscosity 
and  conductivity. 

It  is  obvious  that  any  finite-difference  description 
cannot  describe  explicitly  those  scales  of  motions  which  are 
too  small  compared  to  the  grid  size. 

The  energy  in  the  numerical  models  tends  to  accumulate 
in  the  shorter  wave  range  since  it  cannot  cascade  to  sufficiently 
small  scales  and  dissipate  by  molecular  viscosity. 

Eddy  diffusion  terms  are  thus  used  in  order  to  keep 
the  energy  of  the  short  waves  small  -  and  to  prevent  compu¬ 
tational  instability. 

Under  the  influence  of  the  lateral  eddy  diffusion  the 
shorter  the  wavelength  of  the  perturbation,  the  more  pronounced 
is  the  smoothing  effect. 

Yet  the  shortcoming  of  such  procedures  is  that,  when 
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integrations  are  made  with  lateral  friction  terms  -  solutions 
appear  to  be  too  smooth  if  large  enough  diffusion  coefficients 
are  employed  to  give  computational  stability  and  to  filter 
small  scale  perturbations. 

b )  Review  of  numerical  values  adopted  for  the  lateral  friction 
coefficient 

Approximating  the  turbulent  stress  terms  after  Boussinesq. 
(1877)  by  analogy  with  viscous  stresses  we  obtain: 

(  ^denoting  averaging,  u,  v  being  components  of  velocity 
in  the  horizontal  plane  and  K  defined  as  the  lateral  friction 
coefficient.  Very  little  is  known  about  the  value  of  K.  Most 
authors  assume  K  to  be  a  constant;  others  allow  it  to  vary 
with  time,  height  etc. 

Defant  (1921 )  in  discussing  his  "large  Austausch  process" 

found : 

1/  k  /o  2  -t 

n  y  =  S»  tO  cm  •  .sec 

Richardson  (1926)  found  an  empirical  law  appropriate 
for  estimating  the  dispersion  of  clusters  in  the  atmosphere. 

This  involved  an  eddy  diffusion  coefficient  depending  on  the 
size  of  the  eddies  (or  clusters): 

>Vs 

K  =  0.2  l 
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For  instance,  ii’  we  take  1  to  be  equal  to  the  grid  spacing 
adopted  in  our  model  : 

-  7 

l  -  =  JfOO  Km  -  !0  cm. 

sve  obtain  : 

a  //»  3  2  _/  5  2-/ 

K  *  O.c.^f./O  =  0.3  ,  to .  /O  =J.7'>0  cm-se-c  =  /.7-/0  m  sec 

Since  the  first  formulation,  Richardson’s  law  received  some 
theoretical  justification  from  Kolmogoroff ’ s  similarity 
hypothesis  and  Batchelor  (1950)’ s  work. 

Leith  (196l),  in  a  general  circulation  experiment,  took 
0  value  of  K  -  10  m*  sec/  and  carried  his  experiment 

for  seven  months. 

Malk-o 3  and  Witt  (1959)  used  Richardson’s  relation  for 
a  3mall  scale  convection  experiment  with  K-  0,5  m  -sec 

Smagorinslcy  (1958)  uses  an  eddy  diffusion  coefficient 
varying  with  space  and  time  of  the  form: 

K  =  (kl)2  CDD  ,  k  =  ■££-  1  =  grid  spacing 

[D]  =  two  dimensional  deformation  tensor. 

K  increases  for  smaller  scale  variations  and  for  larger  velocity 
amplitudes.  He  justified  the  choice  of  such  a  coefficient  by 
similarity  considerations  applied  to  v/ind  tunnel  experiments. 

Lilly  (1962)  used  a  similar  coefficient,  which  also 
depended  on  static  stability. 


In  a  numerical  experiment  dealing  v/ith  diffusion  and 
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fallout  of  stratospheric  radioactive  materials  Davidson, 

Friend  and  Seitz  (19 66,  Tellus)  used: 

A=  K  =  K  -  3*/0  cmS sec1  -  2* /O^m  .  &ec[ 

Phillips  in  his  well  known  general  circulation  experi- 

c  p 

ment  (1996)  took  A  =  ICr  m  /sec,  Matsumoto  (i960)  in  a  5  level 
geostrophic  model  assumed  the  same  value. 

Machta  (1955)  assuming  Fickian  diffusion  of  the  dust 
cloud  of  radioactive  debris  from  atomic  explosions  found: 

K-  A  *  O.S'/0S*n.stc. 

On  the  other  hand  Sutton,  (1953)  and  Priestly  by  assuming 

turbulent  regimes  (free  and  forced  convection)  arrive  at 

6  h  p  -1 

values  of  K  =  10  -  10  cm  sec  and  the  same  value  results 

from  Taylor's  statistical  considerations. 


Koo-Tehen-Chao  shows  an  equivalence  between  the  smoothing 

/  2 

method  defined  by: 

and  the  introduction  of  the  lateral  eddy  viscosity  term, 

J  -  being  the  relative  vorticity,  v2  -  finite  difference 
Laplacian  }  lc  -  the  smoothing  coefficient. 

He  obtained: 


-_A 


At 

-  the  time 

step 

Ah 

-  the  grid 

spacing 

A 

-  the  eddy  viscosity 
coefficient 

-  58  - 


8  2  —1 

Do8s  (1956)  assumed  A  =  7.10  cm  sed  .  When 

As  =  300  km  Ai  =  2k  hr  k  =  1  . 

10  2  -1 

Sigryggsen  assumed  A  =  2.10  cm  sec 

S  ^  800  ion  At  =  12  hr  and  thus  k  =  ^  , 

With  At  =  1  hr,  k  =  1/1 3  A  s  =  300  km  we  get 

9  2-1 

:?roin  the  above  relation:  A  =  5.10  cm  sec 


In  an  asymptotic  development  of  the  diffusion  in  the 

Ekman  Layer,  Csanady  (1969)  obtained  a  formula  which  gives 

5  2  -1 

an  eddy  diffusion  coefficient  of  10  m  sec  .  His  formula  is: 
=  M-ir.. .  -  -  S0Sm2*et' )  Uj  =  ySm-&  i0 

It  is  now  generally  accepted  that  A  varies  bety/een 

9  2  —1  5  2  —1 

10^  cm  sec  and  10^  cm  sec 


Governing  equations 

The  model's  equations  and  the  boundary  conditions  will 
refer  to  the  0,250,  500,  750  and  1000  mb  levels  which  will  be 
denoted  by  subscripts  0,1,2,3>U  respectively. 

The  quasi- geostrophic  assumption  is  implied.  The  basic 
dependent  variables  are  the  stream  function  at  the  two  levels 
250  mb  and  750  mb,  denoted  ip,  ,  and  .  Because  of  the 

quasi-geostrophic  assumptions  they  are  related  to  the  geo¬ 
potential  by  means  of  : 

(1)  vp  --4- 

Jm 

for  mean-coriolis  parameter  = 


-  59  - 


/row  the  stream" func tions ,  the  non'd ivergent  component  of 
x'clocities  is  obtained: 


i.  *  Vf 


The  vorticity 


7 


is  given  by: 


lc  =  unit  vector 
along  p-axis 

V  =  horizontal  del 
operator  on  a 

pressure  surf;-<. 

surface 

V  vV 


The  continuity  equation  is: 


(2) 


V 


2cj 

dp 


-  -  vertical  p-  velocity. 

Through  the  hydrostatic  relation  T  can  be  obtained  from 
the  stream-function; 

(3)  ~Tp-  *  ^as  c a) stunt 

The  2  prognostic  equations  are  written  in  terms  of  potential 
vorticities,  denoted  q.j  and  q^. 

The  thermodynamic  energy  equations  cun  be  written: 


(4) 


H 


9  '2A 

TF  'bp 


+  \yiv 


g  4  . 

dp 


din  -0 

-~9p 


-  potential  temperature. 

H  -  the  heat  added  to  unit  mass  in  unit  time. 
Cp  -  specific  heat  at  constant  pressure. 
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The  static  stability 


(5)  c r=2^- 

ctp  op 

is  assumed  to  be  a  function  of  pressure  only. 

Dividing  (4)  by  and  introducing  (5)  and  (3)  we  obtain/ 


(6) 


d 

~?T 


W 


CO  d  _ 


JS_  JL 

Jm,  'Cf»  fi 


Introducing  a  finite  difference  approximation  of  by 

dp 

centred  differences  between  p,  and  ps  and  multiplying  by  P2  we 
obtain  : 


&(?>-rj=-  wfy-%)+  zo?***  + 

jm. 


Hg  is  now  separated  into  a  general  adiabatic  term  H0  and  a 
term  due  to  diffusion: 


(8) 


H2  =  H2a  +  H2d; 


with 

(S) 

where  Ar  is  the  lateral  eddy  diffusion  coefficient  for 
sensible  heatj  thus  ^2d  proportional  to  the  Laplacian 
of  the  temperature  at  500  mbs. 


Equation  (7)  assumes  the  form: 
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In  the  equations  of  motion  "both  surface  friction  and  lateral 
viscosity  are  incorporated. 


The  frictional  stress  X  due  to  the  presence  of  a  lower 
"boundary,  is  assumed  to  vary  according  to; 

(11 )  X 


with: 

02) 

Thus  the  vorticity  equations  at  the  two  levels  and  p2 

will  be; 


(13) 


if'  -/  *-(%;)& 


Av  -  lateral  eddy  viscosity  coefficient. 

f>  -  mean  Rossby  parameter  -  1 .6  X  1011  m”1sec’’1> 

<j  -  acceleration  of  gravity  -  9.81  m  sec”2. 

The  vertical  velocities  in  (13)  are  eliminated  by  substituting 
centered  differences  for  with  the  assumption  of  the 

boundary  conditions  <£/,  =  CJ^ -O  and  substituting 

for  from  the  thermodynamic  energy  equation  : 

(I'd  [£  (W  - %) 'Wfy %)] ' 
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We  assume: 


(15)  Plv  ^  At  =  105 m.&ec  u 


Thus,  we  obtain  the  final  form  of  the  prognostic  equations 
in  terms  of  q^  and  q~: 

(16)  fgfs-N 

f  *v?* 


with 


(17)  f,  ->lfl- fs) 


i-  Jatimnk1 

C  />2* 


In  deriving  (15)  and  (16)  we  assumed 

\y  -  V\  ± 
z  2 

consequently 

^ -v(  f,  -ft--  %ftf,  -  %>%^(\ frfs) 

(t> )  The  dlabatic  heating  due  to  influx  of  sensible  heat  from 
the  bottom  of  the  region. 

The  dlabatic  heating  in  the  model  is  exactly  the  same 
as  assumed  by  Hubs  (1967). 

The  amount  of  heat  entering  unit  column  through  its  lower 
boundary  situated  at  1000  mb,  Qq  is  made  proportional  to 


the  temperature  difference  "between  the  surface  and  the 
air  just  above,  and  to  the  geostrophic  wind: 

(18)  ft  =A»I^I  PJ  ~T"  ) 


Tg  -  the  prescribed  surface  temperature 
Th  -  air  temperature  at  lOOOmb 
1000  mb  wind  speed 
-  a  constant. 


Tg  -  is  constant  with  time  and  Tit-  was  extrapolated  from  the 
temperature  at  the  500  mb  level. 


(19) 


_ zfm _ 

(  k/>3  “  t”/3/) 


(%-%). 


Assuming  a  standard  atmosphere  lapse-rate  between  the 
500  mb  level  and  the  lower  boundary: 


Similarly  we  obtain  the  stream  function  at  1000  mb  level 

(20a)  %  =  9./7x/0STz 


and  the  geostrophic  surface  winds 


(2°t) 


From  (18)  ,  (20),  (20b)  we  obtain  the  value  of  Qs,  Assuming 
Q  to  be  evenly  and  distributed  throughout  the  whole 


column  of  air: 
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Release  of  latent  heat  nas  introduced  in  the  second  experiment. 

The  net  outgoing  radiation  was  asnumud  to  be  constant  over  the 
whole  region  and  an  Instantaneous  heat  balance  was  maintained 
at  each  moment  of  the  integration’s  time.  In  the  first  experiment, 
due  to  the  sy*.i.atry  features  of  the  initial  conditions  and 
the  structure  of  the  model,  this  resulted  in  zero  net  out¬ 
going  radiation.  In  the  second  experiment,  due  to  the 
introduction  of  the  humidity  field  and  the  release  of  latent 
heat,  the  symmetry  tvas  lost  and  the  outgoing  radiation  did 
not  vanish. 

DEVELOPMENT  OP  THE  W  EQUATION  TAKING  INTO  ACCOUNT 
THE  EFFECTS  OF  LATERAL  EDDY  VISCOSITY  AND  THE 

LATERAL  EDDY  DIFFUSION  COEFFICIENT  FOR  HEAT 

,/e  obtain  by  substracting  16  (b)  from  16  (a) 

< » >  ir(iru)-$(%,p) y(\ '%.?,) -?(*> - - 

which  may  be  rewritten  in  the  form 

(23 )  §t  - %Y2x?t(%-%) y &,%)-?&’*)-?&/-'$)- 

Substituting  for  “  Vj} 

dynamic  energy  equation,  i.e.  Eq.  10: 


from  the  thermo 
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h  2>g( v>  ■  %)*  j^tjr  Hio. ->■ 

(24)  ■+e-*R'7Z(H'r%Yg(t<,'/'i)-g(is'%)- ?(*>-**)  - 

~3*-[<-- ),]' vV«  -  f$  +ftv\ -*vu 


Taking  the  Laplacian  of  iiq.  10  yields: 

( 25 )  TT vg( % .%Y*f  A *fr/H^+ 


Subtracting  (25)f rom(2U)  v/e  obtain: 

(z6)  0 = 

-gfa.  %yy*,  -  %)  -j  ^  [§pr( ¥P)31^%  *  AA-  - 

-  ft ir  *?J - y/tf ,  <&}-£& v  4- 3-  -Avyftv?i 

It  is  easily  seen  that  the  terms  involving  A  cancel 
each  other. 


Thus,  the  effects  of  lateral  eddy  viscosity  and  lateral 
eddy  diffusion  of  heat  do  not  change  the  tv>  -equation, 
provided  we  assume  that  the  coefficients  are  equal: 

f\ v  -  Ar  =■  A  =  10  ™Z  sec.  . 

Thus,  v/e  obtain  the  co-equation  in  the  form: 
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(?7)  r,)-g (%■?>)- .*)- 

A  sii.dl:  •  .a  ression  v;ns  used  by  :  atsunoto  (i960)  in 
his  5-level  iuasi-;reostrophi c  model. 
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Boundary  conditions  and  initial  conditions 

The  region  oi'  integration  extends  from  x  =  0  to  x  as  X 
and  from  y  -  0  to  y  =  -Y  .  The  following  averages  are 
defined: 

Off*  >*■;  Off  ^ 

o  ° 

The  following  "boundary  conditions  were  applied  in  the  two 
experiments : 

Cyclical  continuity  of  field  variables  in  the  x-direction, 
i.e.  if  oC  denotes  a  variable: 

(26)  otU)^(o) 

Vpj  =  Vg  =  0  j  S  and  IT  being  the  southern  and  northern 
boundaries  respectively.  With  our  geostrophic  assumption, 
the  stream  function  will,  therefore,  be  constant  along  each 


of  the  boundaries. 

(29) 

(30) 

(IH’fiHf0 

(31) 

tf 

In  order  to  obtain  a  single  valued  solution,  a  single 
value  of  the  stream-function  at  some  grid  point  is  specified 
arbitrarily  (conditions  (30) and  (31  )  have  to  be  modified  in 
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a  finite  difference  formulation  to  ensure  conservation  of 
potential  vorticity  at  each  level).  The  addition  of  the 
lateral  friction  teivn  does  not  interfere  in  the  conser¬ 
vation  of  potential  vorticity  at  each  level. 

<»)  •  Af  jVf*» 

O  0 

According  to  (17)  and  the  boundary  conditions  (28)and(29): 
(33) 

O  O 

Thus  the  total  potential  vorticity  must  be  conserved,  and 
care  must  be  taken  that  this  condition  be  satisfied  in  the 
finite  difference  formulation.  In  the  lateral  friction 
term  the  conventional  4  point  Laplacian  was  used  at  all 
internal  points. 

On  the  boundaries  v/e  write: 

(3h>  2  ?<■> -(Asf  (<iz~ 


The  remaining  terms  were  approximated  by  finite  difference 
schemes  as  described  in  the  report  mentioned  above.  Thus  the 
scheme  for  calculating  the  new  values  of  the  potential  vor¬ 
ticity  on  the  boundaries  becomes: 


(3 


)  2  fa>*  ~  2  t3i'°  If'  9/>  (AS)4(^M  ~  2  faCj,  Q.o  + 
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Initial  conditions 

The  prescribed  surface  temperature,  constant  with  time, 
varies  linearly  in  the  y  direction  with  a  gradient  of  4.5° 
per  grid  interval.  The  temperature  was  assumed  to  he  constant 
along  the  x  direction,  and  had  a  value  of  274°  K  at  the 
central  latitude.  Initially  Tg  =  T4  and  from  the  prescribed 
surface  temperatures  the  initial  temperature  at  500  mb  was 
obtained  as  described  above.  The  stream  function  at  1000  mb 
is  initially  assumed  to  vanish.  From  (20a)  we  obtain  and 

from  (19)  we  obtain 

An  identical  perturbation  is  superimposed  on  the  zonal 
distribution  obtained  in  this  way  both  at  p1  and  ut  p^, 
having  sinusoidal  shape  with  a  wavelength  of  621 6  Km  in  the 
x-direction  and  an  amplitude  of  8.6x1  Cmm /sec  at  the  central 
latitude,  decreasing  sinusoidally  towards  O  at  the  southern 
and  northern  boundaries.  The  numbers  of  grid  points  in  both 
x  and  y  direction  was  15 »  including  boundary  points. 

The  grid  spacing  was  AS  =  444  ion,  the  time  step  being 
A  t  =  1  hour. 

Energy  transformations 

An  insight  into  the  energetics  of  our  system  will  be 
provided  by  the  energy  transformations. 
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The  area  averages  of  the  various  energy  components 
will  "be  denoted  by  K,  K',  P  and  P'  respectively,  lower 
indices  denoting  the  corresponding  pressure  levels. 

If  oC  is  any  field  variable  (scalar  or  vector)  we 
will  write:  e*C  —  ©<£  oC 


where  (  )'  denotes  the  deviations  from  the  zonal  mean. 

V/e  shall  start  with  the  vorticity  equation  at  the  250  mb 
level  (the  procedure  for  the  750  mb  level  being  identical). 


Averaging  zonally  we  obtain: 

(37)  *•>($), i*+/lvt7‘ 


The  term  -  Vj p  vanishes  because  of  the  geostrophic  approx¬ 
imation  and  cyclical  continuity. 

<\/ 

By  multiplying  (37)  by  and  integrating  the  product 
over  the  total  area  of  the  Treasure  surface,  we  obtain  the 
rate  of  change  of  the  kinetic  energy  of  the  zonal  flow. 

Only  the  lateral  friction  term  will  be  treated  here, 
other  terms  having  already  been  treated  in  a  previous  report. 

The  term  resulting  from  lateral  viscosity  is: 

(3M  /)J 


where  L,  is  the  domain's  length  in  the  x  direction 


Integrating  by  parts,  and  applying  the  boundary  condition, 
this  term  becomes 

(39)  ~f\ 

Thus,  including  all  terms: 


(40 


>  I  i  y^y-  4;  ?, 


and  the  rate  of  change  of  the  kinetic  energy  of  the  zonal  flow 
in  r-  pressure  interval  o rdcL  becomes: 


And  similarly  at  the  750  mb  level: 


If  we  take  the  two  levels  to  be  representative  of  the 
upper  and  lower  half  of  the  model  atmosphere  respectively, 
the  rate  of  change  of  the  zonal  kinetic  energy  will  be  given 
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and  the  corresponding  average  per  unit  pressure  interval: 


A  similar  procedure  is  employed  to  obtain  the  rate  of  change 
of  the  zonally  averaged  potential  energy.  V/e  begin  by  zonally 
averaging  the  thermodynamic  energy  equation: 


Z  jm  J'”  f*  d 

Multiplying  it  by  (V,  ~  and  integrating  over  x  and 

y  we  obtain: 


±(Y, eg «*(?,-%)+ 

^cr  ~ ft)] 


Jm  ‘-A 

Multiplying  by 


J*r_.hE 

eiP’$  '  9 


we  obtain  the 


rate  of  change  of  the  zonally  averaged  potential  energy  per 
unit  column: 


9P 

(47)  'Zt 


V-Kfir 


-  ft 


<9 

i, 


The  vorticity  equation  for  the  eddies  is  obtained  by  subtract- 
ing  (36)  from  (37): 


(48) 


^Vtz- 


-  73  - 


The  rate  of  change  of  the  mean  kinetic  energy  of  the 
eddies  is  obtained  by  multiplying  (43)  by  and  integ¬ 

rating  the  products  over  x  and  y.  For  the  lateral  eddy- 
viscosity  term  we  have  after  integration  by  parts,  and  ap¬ 
plication  of  the  boundary  conditions: 


(49)  vZj'  dc  =  -  dC 


Finally  we  obtain  : 


—  2 — L.  y..j&  -  J—m  v ') _ dpL  u/'cj'  — 

3 1  '  3 1  3  £  ~  9  <  9y '  '  '  '  9  f>2  T1 '  2 


(50)  dt  ’•  Vt  £  £  #  ’ey 

9C 


~  -  f(}!) 

And  similarly  for  the  750  mb  level.  As  above: 

(51>  &«'=*£  &(«) 


For  deriving  the  rate  of  change  of  the  potential  energy  of 
the  eddies,  we  subtract  (45)  from  the  energy  equation: 


(52)  2- 2 )--v„v(i'r't'sy^v (?;-%)+  + 

-  ~Tc  W?a+AvV/-ri) 

Jm  Lls 

Multiplying  (52)  by  (  \f/,  -  ^  )  and  integrating  we 

.  2J2-  2L-  2s. -  2fu/,~  W// )C-_  dm  -  (J  \  + 

(53)  9 i  'at  cra  gvri  rs)  -  c'zp2^,  t/z  '  '  s' 


v»r??  «■  *  v2r?9- 

+  -^r  2222) +  ^ps%  H‘ 
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So  that  using  Phillips*  notation  (1956)  we  have  the  following 
transformations: 
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PHS  DgVgLOPg.lEI-T  OF  THE  CIRCULATION 

The  first  experiment  to  be  denoted  as  Experiment  I, 
does  not  include  water  vapor,  and  is  similar  to  Experiment 
I  by  Huss  (1967),  except  that  lateral  viscosity  has  been  in¬ 
cluded.  Like  in  the  corresponding  experiment  mentioned  above, 
tie  disturbance  was  found  to  grow  initially  at  an  exponential 
te.  Later  on  this  growth  ceased  and  the  appearance  of  hi$ier 
harmonics  was  observed. 

In  fig.  (l)  we  observe  the  variations  of  the  amplitudes 

cf  the  first  wave  components  of  the  stream  function  at  the 

central  latitude,  as  a  function  of  time.  The  amplitude  of 

the  first  component  at  both  levels  reached  its  first  maximum 

after  about  13  days,  the  numerical  value  attaining  9.3 X  10  m  /sec 

at  250  mb,  almost  half  the  value  obtained  by  Huss  in  the  cor- 

7  2 

responding  experiments,  which  v/as  1.7  X  10  m /sec.  As  we 
-hall  frequently  refer  to  Huss'  corresponding  first  ex¬ 
periment  we  shall  call  it  thereafter  the  "previous  run", 
without  any  other  specification. 

After  reaching  the  maximum  the  amplitude  decreased  for 

5  days,  then  continued  to  rise,  reaching  a  second  maxima  at 

7  2 

the  25th  day  amounting  to  1 ,hZ  X  10  m  /sec  at  the  250  mb  level. 

No  pronounced  oscillations  as  in  the  previous  run  were 
observed.  After  a  minimum  on  the  27th  day  of  1.05  X  107  m2/sec 
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/L  rj 

at  250  mb  and  of  1.3  X  10  m  /sec  at  the  750  mb  level,  the 

amplitudes  continued  to  rise  until  the  35th  day.  Following 

this  date  the  amplitudes  at  both  levels  remained  more  or 

7  2 

less  stable  with  2.65  x  10  m  /sec  at  the  250  mb  level  and 
6.4  x  10u  nv/sec  at  750  tub,  retaining  those  values  until 
the  end  of  the  run  on  the  42nd  day.  Generally  speaking  the 
latest  values  attained  seemed  to  be  half  the  average  values 
of  the  previous  run  during  the  corresponding  stages.  The 
smoothing  is  due  to  the  action  of  the  lateral  friction  term, 
now  introduced. 

Fig. (2) shows  us  the  amplitudes  of  the  odd  harmonic 
components  of  the  disturbance  at  the  central  latitude  of 
the  upper  reference  level. 

The  component  with  wave  number  (I)  initially  imposed 
on  the  flow  remains  dominant  throughout  the  run.  Compared  to 
the  previous  run,  a  considerable  smoothing  is  observed  - 
in  the  amplitudes  of  the  shorter  waves.  Still,  the  components 
with  odd  wave  numbers  are  much  greater  than  those  with  even 
wave  numbers.  On  the  10th  day  only  the  3rd  component  exceeds 
10  m  /sec  while  the  5th  has  a  value  of  about  10  m  /sec. 

The  7th  component  reaches  this  value  only  on  the  25th  day. 

On  the  20th  day  the  fundamental  wave  (no.  I)  has  a 

6  2 

minimum  consisting  of  6.4  X  10  m  /sec  while  the  3rd  har- 

5  2 

monic,  the  value  of  7  .  10^  m  /sec,  the  5th  a  value  of 
7.5k  10^  m2/sec  and  the  7th  a  value  of  5.5  X  10^  m2/sec. 
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Vie  observe  that  they  are  one  or  two  magnitudes  smaller  than 
the  fundamental  wave.  Thus,  the  noise  level  v/as  considerably 
reduced  and  lowered  in  the  present  model  compared  with  the 
previous  run: 

On  the  last  day  of  the  run  the  3rd  component  attained 
the  value  of  2  X,  ^0'>  m2/sec,  the  3th  3*2  X  10^  m2/sec  and 
the  7th  3X  10  m2/sec. 

Fig.  3  shows  us  the  amplitudes  of  the  various  harmonic 
components  of  the  disturbance  with  even  wave  numbers  at  the 
central  latitude  of  the  upper  reference  level. 

The  amplitudes  of  the  even  harmonics  never  exceeded  the 
h  2 

value  of  10  m  /sec  and  the  highest  value  in  the  run  was 
attained  by  the  2nd  component  on  the  34th  day;  4X  10^  m2/sec. 
The  even  harmonics,  amplitudes  were  one  or  two  magnitudes 
smaller  than  the  amplitudes  of  the  corresponding  odd  com¬ 
ponents. 

On  the  42nd  day,  the  2nd  harmonic  attained  1  .4  X  10^rn2/sec, 
the  4th  attained  7X102  m2/sec  and  the  6th  4x  102  rn2/sec. 

This  indicates  a  drastic  smoothing  if  we  compare  our 
values  to  those  of  the  previous  run,  thus  confirming  our 
expectations  as  to  the  effect  of  the  inclusion  of  the  lateral 
eddy  viscosity  term. 

The  latitude  pressure  distribution  of  the  poleward  eddy 
momentum  transports  is  shown  for  4  selected  days  in  fig.  6. 

"he  values  of  uV  must  be  antisymmetrical  around  the 
central  latitude  -  therefore  the  poleward  flux  vanishes  there. 
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<~n  t'i*  10tr.  duy  the  flux  is  consistently  toward  the 

central  latitude  at  the  2  levels  and  at  all  latitudes.  Ab- 

o  .  9 

solute  maximum  values  were  of  the  order  of  11.5  m /sec  com- 

9  •*> 

pared  to  20  r:  /sec  in  the  previous  run  for  the  corresponding 

2  2 

period,  on  the  15th  day  maxima  of  bO  m  /sec  were  attained, 
the  highest  values  during  the  entire  run 

On  this  date  small  regions  in  the  extreme  South  and 
"orth  were  characterized  by  transports  toward  the  boundaries. 

On  the  20th  day  the  distribution  of  the  fluxes  was  quite 
Jiffere.it  indicating  a  negative  transport  in  the  south  and 
positive  transports  in  the  north.  Their  maximum  magnitude 
was  11.2  m  /sec  .  In  small  regions  in  the  extreme  north  and 
south,  the  flux  was  directed  away  from  the  boundaries  in 
opposite  direction  of  the  hulk  of  the  flux  in  each  half  of 
the  atmosphere.  Such  an  opposite  transport  appeared  also  in 
the  central  region . 

On  tne  yOth  day  the  fluxes  dropped  to  smaller  magnitudes 
but  the  characteristic  feature  of  positive  transports  in  the 
southern  half  and  negative  ones  in  the  northern  half  reappeared. 

2  2 

The  maxima  are  10.5  m  /sec  at  1  »  10  at  the  upper  level. 
The  transports  opposite  in  direction  to  the  bulk  of  the  flow 
disappeared  along  the  lateral  boundaries. 

On  the  both  day  values  increased  again  reaching 
2  2 

35.1  n  /sec  at  the  upper  level. 
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Transports  in  opposite  directions  are  again  observed  at 
the  lower  reference  level  at  3  =  6£-  and  d  =  9a  • 

The  change  with  time  of  the  various  components  of  the 
specific  kinetic  energy  is  shown  in  fig.  9. 

The  zonal  kinetic  energy  decreased  at  the  upper  reference 

2  2 

level  from  the  value  of  536  m  /sec  on  the  first  day  to  the 
value  of  387.5  on  the  40th  day.  The  decrease  was  accompanied 
by  small  oscillations.  A  tendency  of  stabilization  can  be 
perceived  around  the  42nd  day,  which  is  the  last  day  of  the 
present  run. 

At  the  lower  level  a  slight  decrease  of  the  zonal  kinetic 
energy  is  observed  until  the  22nd  day,  followed  by  a  slight 
increase.  Towards  the  end  of  the  run,  a  tendency  towards 
stabilization  may  be  indicated. 

At  the  upper  level,  the  eddy  kinetic  energy  increased 

2  2 

while  oscillating  and  attained  its  highest  maximum  of  137  m  /sec 
on  the  3<3th  day. 

The  eddy  kinetic  energy  at  the  lower  level  was  v^ry  small. 

It  appears  to  increase  somewhat  with  time  and  attains  a  value 
2  2 

of  10  m  /sec  on  the  37th  day,  thereafter  stabilizing  at  this 
value. 

The  total  averaged  kinetic  energy  decreased  very  slightly 

p  p 

during  the  formative  stage  from  282  m  /sec  on  the  16th  day 

2  2 

and  then  more  abruptly  to  248  m  /sec  on  the  22nd  day.  After 
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this  date  the  energy  increases  again,  attaining  a  maximum 
of  280  m  /sec  on  the  37th  day.  This  is  followed  "by  almost 
constant  values  and  apparently  a  steady  state  has  been  ob¬ 
tained. 

The  variation  of  the  zonally  averaged  flow  at  1000  mb 
as  a  function  of  latitude  and  time  is  shown  in  figure  5. 

The  development  was  considerably  slower  than  in  the 
previous  experiment.  A  value  of  4.20  meter/sec  v/as  observed 
at  the  central  latitude  on  the  16th  day,  and  this  v/as  the 
highest  value  of  the  entire  run. 

The  boundaries  between  the  central  latitude  and  the 
regions  of  the  easterlies  which  developed  were,  relatively, 
more  smooth  during  the  entire  course  of  the  integration, 
being  located  in  the  South  between  j  =  4  and  j  =  and 
in  the  Uorth  between  j  =  10  and  ,1  =  12. 

The  speeds  in  the  belts  of  the  easterlies  and  in  the 
belt  of  the  westerlies  were  considerably  lower  than  in  the 
previous  run. 

The  maximum  wind  speed  decreased  from  4.20  n/sec  at 
the  central  latitude  on  the  16th  day  to  0.22  nv/sec  on  the 
22nd  day,  then  rose  again  and  reached  a  maximum  on  the  27th 
day,  of  3.71  m/sec,  again  decreasing  to  1.60  m/sec  on  the  32nd 
day  and  rising  again  to  a  local  maximum  of  3*^5  n/sec  on  the 
38th  day.  These  variations  might  indicate  a  cycle 
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of  11  days  at  1000  mb,  but  we  cannot  draw  such  a 
conclusion  from  only  42  days  of  integration  and  it  seems 
that  afterwards  this  clear-cut  cycle  disappears  due  to  the 
increasing  effect  of  the  eddy  diffusion  terra,  as  it  is  seen 
in  the  amplitude  of  the  disturbance. 

The  highest  speed  of  easterlies  appeal’s  betweeu  the  36th 
-  42nd  days,  reaching  speeds  not  exceeding  2  m/sec  (the 
highest  speed  being  2,13  m/sec  on  the  38th  day  =  the 
same  day  on  which  the  highest  local  maximum  of  3.65  m/sec 
was  observed  in  the  westerlies ) . 

The  local  maxima  of  easterlies  appear  on  the  same  days 
as  the  maxima  in  the  westerlies.  Their  intensities  being 
however  lower. 

Various  energy  transformations  are  shown  in  figures  7, 

8,8*. 

The  highest  value  of  {p,Pf$  amounting  to  4.2  x  10  ^  kj  ra 
sec-1  cb”1  was  attained  on  the  37th  day.  This  was  the  third 
maxima  during  the  entire  run,  the  first  being  observed  on  the 
14th  and  the  second  on  the  25th  days,  with  values  of  1.88  x 
10"3  and  3.14  x  10  3  kj  m  sec  cb  respectively.  During 
the  course  of  the  run  the  values  were  found  to  drop  down  to 
0.12  x  10-5  on  the  19th  day.  In  this  run  much  lower  values 
were  attained  than  in  the  previous  run. 


The  conversion  of  heat  into  zonal  potential  energy  is 
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consistently  positive  and  into  eddy  potential  energy,  con¬ 
sistently  negative. 

The  positive  conversion  in  this  experiment  was  larger 
than  the  negative  conversions. 

Friction  causes  a  reduction  of  kinetic  energy.  The 
variations  in  the  tranformation  of  kinetic  energy  into 
friction  are  not  exactly  in  phase  with  the  other  trans¬ 
formations,  attaining  maxima  on  the  16th,  27th  and  39th  days. 
Thus  having  a  phase  lag  of  two  days  relative  to  {p,  p'>. 

The  absolute  values  are  smaller,  by  a  factor  of  3  than  the 
energy  gained  by  heating.  Frictional  losses  are  shown  ac¬ 
cording  to  a  scale  reduced  by  10  compared  to  the  other 
transformations.  Their  values  are  much  lower  than  those 
obtained  in  the  previous  run.  The  highest  value  of  {K,F  } 
was  attained  on  the  39th  day  having  a  value  of  7.5  x  10~^ 
compared  to  24  x  10  the  absolute  maximum  attained  in  the 
previous  run. 

During  this  run  most  of  the  eddy  kinetic  energy  was 
transformed  into  kinetic  energy.  Oscillations  that  indicated 
small  negative  values  were  obtained  on  the  22nd  and  the  33rd 
day  but  the  negative  values  were  very  small  compared  to  the 
dominant  positive  ones. 

In  fig.  8  the  frictional  losses  of  kinetic  energy  due 
to  lateral  eddy  viscosity  are  shown. 


The  values  of  these  transformations  are  smaller  than 


-  83  - 


the  values  of  all  other  transformations  taking  part  in  the 
energy  exchanges. 

-7  -2  -1 

Starting  with  127  X  10  kj  m  sec  on  the  first 

day  of  the  run  j  K,  A  ]■  decreased  until  the  15th  day,  when 

-7  -2  -1 

it  reached  the  values  of  32  X  10  k;j  m  sec  .  A 

local  maximum  was  attained  on  the  18th  day,  a  second  one  on 

the  29th  day  and  values  appear  to  become  stable  on  the  40th 

-7  -2  -1 

day,  v/ith  43  X  10  kj  m  sec  .  The  values  of  the  first 

—7  —7  —p  —1 

two  maxima  were  40  X  10  and  41  X  10  kj  m  sec  res¬ 
pectively. 

The  frictional  losses  of  eddy  kinetic  energy  rose  slov/ly 

-7  -2  -1 

from  the  first  day  to  a  maximum  of  40  X  10  kj  m  sec 

—7  —2  —I 

on  the  16th  day.  A  second  maxima  of  44X10  kj  m  sec 

was  observed  on  the  27th  day,  and  a  higher  maximum  of  63  X 
-7  -2  -1 

10  kj  m  sec  was  reached  on  the  37th  day.  While  in  the 
first  half  of  the  run,  i.e.  until  the  24th  day  the  values 
of  dominate,  in  the  second  half  i.e. 

from  the  24th  day  until  the  end  of  the  run  on  the  42nd  day 
the  values  of  dominate.  Towards  the  end 

of  the  run  a  tendency  towards  stabilization  can  be  observed 
in  the  values  of  both  transformations. 

The  zonally  averaged  diabatic  heating  as  a  function 
of  latitude  is  shown  in  fig.  io  for  the  10th,  20th,  30th  and 
40th  days. 

The  magnitudes  wer e  considerably  lower  than  in  the  previous 
run  and  the  distribution  was  also  much  smoother. 
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or.  the  10th  nay  the  heating  was  still  very  small,  almost 
imnercc  Jtible. 

On  the  20th  day  it  was  not  much  higher,  but  on  the  30th 
day  the  r-jner.-.l  pattern  could  be  easily  observed. 

The  region  of  most  extreme  heating  was  in  the  south 
arr-und  5  =  2.  Proceedin'-  to  the  North,  the  vaLues  became 
negative  between  5  =  k  and  j  =  5*  No  diabatic  heating  was 
obtained  at  the  central  latitude,  due  to  the  constraints  of 
the  model. 

bmall  positive  values  v?ere  obtained  hotiever  between 
5=8  and  5  =  11%.  Intense  cooling  was  observed  in  the 
north  with  Its  highest  value  at  5  =  Ik. 

No  simplification  of  the  distribution  ocurred  on  the 
UOth  day  -  contrary  to  the  previous  run. 

.faile  the  values  were  lower,  the  same  pattern  persisted, 
so  that  most  of  the  heat  was  supplied  to  the  atmosphere  at 
the  very  low  latitudes  while  the  sink  was  located  in  the 
extreme  north. 

The  zonally  averaged  zonal  flow  is  shown  in  fig.  h  as 
a  function  of  latitude  and  pressure  for  k  selected  days 
at  intervals  of  10  clays. 

The  same  tendencies  ns  in  the  earlier  run  towards  the 
formation  of  a  "jet"  at  the  central  latitude  at  the  formative 
stage  of  the  run  are  observed.  (The  jet  was  slightly  "split"). 
On  the  10th  day  the  velocity  at  p^  at  the  central  latitude 
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was  36.1  m/sec  while  the  maxima  located  at  3  =  6,  and  3  =  11 
are  36.7  ny'sec.  A  less  pronouncfed  similar  effect  was 
observed  at  the  lower  level.  On  the  20th  day  the  splitting 
remained  evident  -  the  velocity  at  the  central  latitude 
was  37  rr/sec  and  maxima  of  40.1  m/3ec  were  located  at  j  =  6 
and  0  =  10. 

On  the  30th  day  the  splitting  tended  to  disappear.  The 
velocity  at  at  the  central  latitude  was  41.4  rn/sec  while 
the  2  maxima  of  42.3  m/sec  were  observed  near  it  at  3  =  7 
and  3=9.  On  the  40th  day  a  single  well  developed  central 
jet  could  be  observed  with  a  speed  of  48.2  m/sec  .  The  de¬ 
velopment  of  the  easterlies  proceeded  at  a  sLower  rate.  On 
the  20th  day  easterlies  were  found  slightly  below  730mb  level. 
Their  lateral  extent  was  constant  during  the  run,  not  penetrating 
beyond  3=4  and  3=12  , 

On  the  40th  day  easterlies  penetrated  up  to  the  500  mb 
level . 

The  maxima  of  the  speeds  of  the  easterlies  did  not 
exceed  during  the  entire  course  of  the  run,  the  value  of 
2  m/sec.  The  highest  maximum  was  attained  by  the  wester¬ 
lies  on  the 40th  day  reaching  a  value  of  48  m/sec. 


Description  of  the  development  in  the  second  Experiment 

The  amplitude  of  the  first  harmonic  of  the  stream 
function  at  the  central  latitude,  at  the  two  reference 
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levels  is  shown  in  fig.  12.  The  stage  of  initial  growth  was 

*/ 

similar  to  the  corresponding  stage  in  experiment  I  until 
the  13th  day.  The  maximum  attained  on  the  13th  day  had  a  value 
of  9.8  X  10°  m2/sec  at  the  250  mb  level,  and  of  6  X  10^  m2/sec 
at  the  750  nib  level.  These  maxima  are  almost  half  of  the  cor¬ 
responding  maximum  values  obtained  on  the  13th  day  in  the 
previous  run. 

No  small  amplitude  oscillations  were  observed  in  this 
experiment. 

The  pronounced  variation,  which  may  presumably  be  identi¬ 
fied  as  an  "index  cycle”  which  appeared  in  the  previous  run  - 
was  also  not  present  in  this  run,  although  one  long  period 
oscillation  was  observed,  with  a  period  of  10  days.  This  may 
pei’haps  be  a  remnant  of  the  former  index  cycles,  all  later 
cycles  having  been  obliterated  by  the  strong  effect  of  the 
lateral  eddy  viscosity  term.  The  18th  day  presents  the  lowest 
value  of  this  oscillation,  with  2.8  X  10°  m2/sec  (at  250  mb 
and  1.9  X  10  m /sec  at  750  mb).  The  minimum  of  the  present 
run  is  approximately  1/3  of  the  corresponding  values  in  the 
previous  run. 

The  following  minimum  was  observed  on  the  28th  day, 

7  2  62 

with  1 .4  X  10  m  /sec  and  1.6  X  10  m  /sec  respectively  at 
the  250  and  750  mb  levels. 

The  750  mb  minimum  on  the  28th  day  is  somewhat  lower 
than  the  corresponding  value  at  the  beginning  of  the  same  cycle. 
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a  detail  that  characterized  also  the  first  cycle  in  the 
previous  run. 

Prom  the  23  Lh  day  on  -  our  amplitudes  depart  from 

the  "index  cycle"  pattern  which  was  typical  of  the  previous 

run.  We  observe  a  slow  hut  constant  growth  at  both  levels 

at  a  roughly  exponential  rate,  until  the  51th  day  of  the 

7  2 

integration  when  the  values  attained  were  5.04  10  m  /sec 

7  P 

at  the  250  nib  level  and  2.0  K  10'  m^/sec  at  the  750  mb  level. 
From  the  51st  day  until  the  end  of  the  run,  on  the  65th  day, 
the  amplitudes  at  both  levels  became  more  or  less  stabilized 
at  the  values  attained  on  the  51st  day  -  indicating  consider¬ 
able  smoothing  of  the  first  harmonic  of  the  disturbance  at 
both  levels.  Yet,  those  stabilized  values  corresoond  roughly 
to  the  mean  amplitudes  obtained  at  both  levels  in  the  previous 
run. 

In  figs.  13,14,15,16  the  change  with  time  of  the 
amplitudes  of  the  various  harmonic  components  of  the  dis¬ 
turbance  at  the  central  latitude  of  the  upper  reference  level 
are  shown.  Figs.  14  and  16  show  the  components  with  odd  and  even 
wave  numbers  respectively  in  the  previous  run,  while 
figs.  13  and  15  show  the  corresponding  components  in  the  present 
experiment. 

In  the  previous  run  components  with  even  wave  number,? 
attained  higher  values  than  those  with  odd  wave  numbers.  In 
the  present  ran,  a  considerable  smoothing  was  observed  in  the 
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case  of  both  odd  and  even  wave  components.  In  this  run  the 
components  with  odd  wave  numbers  also  attained  higher  values. 

:.’ig.  14  show 3  the  components  with  odd  wave  number  of 
exocri'  ent  III  of  the  previous  run.  The  3rd  harmonic  ex¬ 
ceeds  105  m2/nec  on  the  4th  day,  the  5th  harmonic  on  the 
3th  day  and  the  7th  harmonic  on  the  11th  day. 

In  fig.  13  the  corresponding  odd  wave  number  components 
of  the  present  experiment  are  shown.  Here  the  3rd  harmonic 
exceeds  ICp  m2/sec  on  the  4th  day,  the  5th  harmonic  on 
the  9th  day  and  the  7th  harmonic  on  the  15th  day. 

In  experiment  III  of  the  previous  run,  the  3rd  harmonic 

5  2 

attained  a  first  maximum  of  1.3  X  10  m  /sec  on  the  16th  day 
while  in  the  present  run  on  the  same  day  the  3rd  harmonic 

r  o 

attained  3.1  X  10^  m  /sec.  The  same  harmonic  attained  its 
highest  value  of  5.5  X  10°  m2/sec  on  the  51st  day  in  the 

g 

previous  run,  and  it  attained  its  highest  value  of  3.0  X  10 
ni  /sec  on  the  present  run  on  the  56th  day. 

6 

The  5th  harmonic  attained  a  first  maximum  of  2.1  X  10 

2 

m  /sec  on  the  20th  day  in  the  previous  run  while  in  the 
present  experiment  its  first  maximum  9.5  X  10  m  /sec  is 
attained  on  the  22nd  day.  The  highest  value  attained  by  the 
5th  harmonic  in  the  previous  run,  2.3  X  10  m  /sec,  was  ob¬ 
served  on  the  54th  day  while  in  the  present  run  the  highest 
value  of  the  same  harmonic  6  X  10  m  /sec  was  attained  on 
the  60th  day.  The  7th  harmonic  attained  a  first  maximum  of 
4.3  X  10m  /sec  on  the  22nd  day  in  the  previous  run,  while 


-  89  - 


in  the  present  experiment  the  same  harmonic  had  its  first 

h  p 

maximum  on  the  22nd  day  also,  with  a  value  of  9.2  X  10  m  /sec 

5  2 

The  highest  value  attained  "by  the  ?th  harmonic,  8  \  10  m  /sec 

was  observed  in  the  previous  experiment  on  the  48th  day  while 

5  2 

in  the  present  model  the  highest,  value  of  2.8  X  10  m  /sec  was 
observed  on  the  58th  day. 

Pig.  16  shows  the  components  with  even  wave  number  of 
Experiment  III  of  the  previous  run,  while  fig.  15  shows 
the  corresponding  components  of  the  present  experiment. 

Y.'e  see  that  the  higher  harmonic  components  are  more 

attenuated  than  the  lower  components.  In  the  previous  ex- 

6  2 

periment  the  2nd  component  attained  a  value  of  10  m  /sec 
at  the  first  noticeable  maximum  on  the  22nd 
day  while  in  our  experiment  it  attained  a 

u  2 

maximum  of  6.5  x  10J  m  /sec  on  the  27th  day.  The  highest 

7  2 

value  of  1.2  \  10  m  /sec  was  attained  by  the  2nd  component 
in  the  previous  experiment  on  the  4-0 th  day,  while  the  highest 
value  in  the  present  run  was  8.5  \  m2/sec. 

In  the  previous  experiment  the  4th  harmonic  attained  the 

6  2 

first  maximum  and  of  1.8  \  10  m  /sec  on  the  22nd  day  while 

on  the  same  day  the  same  harmonic  in  our  present  experiment 

L  2 

attained  the  value  of  9.5  x  10  rn  /sec. 

The  highest  value  attained  by  the  4th  harmonic  in  the 

6  ? 

previous  experiment  of  4.2  X.  10  m  /sec  was  observed  on  the 

44th  day,  while  in  the  present  experiment  the  highest  value 
5  2 

was  6X10  m  /sec  and  it  was  observed  on  the  55th  day. 
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In  the  previous  experiment  the  6th  harmonic  attained 

5  2 

its  first  maximum  7.5  x  10  m  /sec,  on  the  23rd  day,  while 

U  2 

in  the  present  run  it  attained  a  first  maximum  of  k  X  10  m  /sec 
on  the  26th  day.  It  attained  its  highest  value  of  3.9  X  10^ 
m2/see  in  the  previous  run  on  the  59th  day.  While  in  the 
present  run  it  attained  its  highest  value  of  1  .2  X  ICr  m  /sec 
on  the  57th  day. 

Thus  again  the  effect  of  the  lateral  viscosity  is  strongly 
marked. 

.’lie  variation  of  the  zonally  averaged  zonal  component  of 
the  wind  at  1000  mb  as  a  function  of  latitude  and  time  is 
shown  in  fig.  17. 

The  same  fundamental  pattern  subsists,  westerlies  dom¬ 
inating  the  central  region,  while  the  easterlies  were  con¬ 
fined  to  the  north  and  the  south. 

Speed  maxima  are  lower  than  in  the  previous  run  and  the 
distribution  differed  considerably. 

Maxima  of  the  westerlies  are  observed  in  the  present 
model,  on  the  16th,  27th  and  on  the  46th  day.  The  index  cycle 
was  no  longer  a  prominent  feature  of  the  developing  circulation 
although  traces  of  it  may  perhaps  be  detected  during  the 
earlier  stages  of  the  run. 

Between  the  maxima  of  the  westerlies  their  speeds 
go  down  to  25,3  of  the  peak  values  between  the  first  two  maxima. 
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and  to  50/o  of  the  peak  values  "between  the  second  and  the 
third  maxima.  The  values  of  the  maxima  attained  by  the 
westerlies  on  the  16th,  27th,  43th  and  64th  days  were 
4.99  m/sec,  5«07  n\/sec,  6.47  m/sec  and  3.10  m/sec  res¬ 
pectively.  They  are  by  50,u  lower  than  the  corresponding 
maxima  obtained  in  the  previous  run. 

Intensification  of  the  easterlies  in  the  south  was 
practically  in  phase  with  the  changes  of  the  westerlies  in 
the  central  belt. 

All  the  three  maxima  of  the  easterlies  in  the  north 
coincided  in  time  with  the  maxima  of  the  westerlies  in  the 
central  belt.  In  the  south  only  two  maxima  were  observed, 
on  the  16th  and  43th  day  during  which  the  westerlies  also 
had  reached  maximum  values.  The  first  maximum  of  the  easter¬ 
lies  in  the  north  had  a  higher  value  than  the  corresponding 
maximum  of  the  easterlies  in  the  sovith.  On  the  16th  day,  the 
maximum  speed  attained  in  the  north  was  located  at  j  =13 
the  extremal  velocity  was  1.65m/sec  while  the  corresponding- 
value  in  the  south  at  0=3  was  1.31  m/sec. 

The  second  speed  maxima  in  the  north  at  the  same  latitude 
of  d  =  13  indicated  a  velocity  of  2.06  n/sec  while  in  the 
south  the  corresponding  value  at  d=3  (which  was  not  a  speed 
maximum)  attained  3.53  m/sec. 

The  third  speed  maximum  of  the  easterlies  in  the  north 
was  obtained  at  the  same  latitude  of  j=13  on  the  43th  day 
with  a  velocity  of  5.75  m/sec  while  in  the  south  a  simultaneous 
maximum  was  observed  at  j=3  with  4.73  iVsec. 
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Compared  with  the  previous  experiment  the  values  of  the 
easterlies  were  much  lower. 

A  northward  shift  of  the  central  belt  is  no  longer  ob¬ 
tained  in  the  present  experiment  although  at  the  end  of  the 
run  such  a  tendency  appeared  to  develop.  Consequently  all  the 
maxima  in  the  westerlies  were  observed  at  the  central  latitude 
j  =  8  and  only  on  the  64  th  day, one  day  before  the  conclusion 
of  the  present  run  the  local  maximum  of  the  westerlies  was 
located  at  j  =  9  • 

An  additional  aspect  of  the  structure  of  the  zonal  flow 
i s  presented  by  the  pressure  latitude  distribution  of  LL  for  4 
selected  days,  as  shown  in  fig.  18. 

I  have  chosen  the  4  days  which  were  the  more  significant. 
Some  of  the  latest  stages  of  the  run  were  selected,  as  in 
experiment  I.  They  resembled  the  pattern  ali*eady  obtained  in 
that  experiment. 

After  the  48th  day  -  the  velocity  in  the  axis  of  the  jet 
reached  42.5  m/sec  and  on  the  54th  day,  31.2  m/sec.  The  cor¬ 
responding  values  on  the  60th  day  and  68th  day  were  of  27.7 
m/sec  and  27.5  n/sec.,  respectively. 

The  most  intense  velocities  in  the  easterlies  in  the 
south  were  4.7  ir/sec  on  the  40th  day,  4.0  ny'sec  on  the  54th  day, 
3.5  m/sec  on  the  60th  day  and  2.8  ny'sec  on  the  65th  day.  No 
indication  of  an  index  cycle  such  as  in  the  previous  experiment 
could  be  detected.  In  the  North  the  easterlies  pehetrated 


-  93  - 


vertically  beyond  the  500  mb  level  on  the  48th  day.  On  the 
other  3  days,  chosen,  they  extended  only  to  a  limited  extent 
beyond  the  750  mb  level.  Also  in  the  north  the  highest  speeds 
were  located  near  the  surface.  The  largest  speed  in  the 
easterlies  in  the  north  among  those  4  days,  amounted  to  -5.8 
m/sec  on  the  48th  day.  The  corresponding  value  on  the  54th  day 
was  -4.3  m/see..  These  two  values  were  higher  than  those  at¬ 
tained  by  the  easterlies  in  the  south  on  the  3ame  days.  On  the 
6oth  day,  the  easterlies  in  the  north  attained  -2.5  m/sec  and 
on  the  65th  day  -2.2  m/sec.,  i.e.,  lov/er  than  the  corresponding 
values  in  the  south  on  the  same  dates. 

On  the  60th  and  65th  days,  the  vertical  axis  of  the  wind 
speed  maxima  in  the  belt  of  the  westerlies  tilted  with  height 
towards  the  south  from  3=9  at  1000  mb  to  3=6  at  250  mb 
on  the  60  th  day,  and  from  j  =  9  to.  3  =  5  at  the  same 
levels  on  the  65th  day.  A  small  tilt  of  the  axis  towards  the 
north  from  3=8  at  1000  mb  level  to  j  =  9  on  the  250  mb 
level  was  observed  on  the  54th  day. 

Pig.  19  shows  us  the  patterns  of  the  zonally  averaged 
northward  eddy  transport  of  momentum.  The  same  days  quoted 
above  were  chosen.  Here  also,  lilce  in  the  previous  experiment, 
large  variations  from  one  day  to  another  v/ere  observed.  On 
the  48th  day,  transports  toward  the  south  prevailed  over  the 
northern  half  of  the  region  and  a  part  of  the  southern  one. 
Northward  transports  were  observed,  in  the  southern  region, 
and  at  250  mb  in  certain  parts  in  the  northern  region.  A 
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2  2 

maximum  of  southward  transports  -  43.3  m  /sec  v/as  obtained 

on  the  UCth  day  at  the  upper  level  at  j  =  10.  The  maximum 

2  2 

of  the  northward  transport,  31  .2  m  /sec  at  j  =  5  on  the  same 

i j 

day,  was  observed  also  at  the  upper  level.  On  the  54th  day  the 

pattern  -.ins  additional  complexity.  The  northward  positive 

values  extend  now  well  into  the  northern  half  of  the  region, 

attaining  a  maximum  value  of  43*0  rn  /see  at  the  upper  level 

at  j  =  6.  The  southward  negative  values  attain  a  maximum  of 

-  23  n2/sec2  at  the  upper  level  at  j  =  12.  At  the  750  mb 

level  on  the  same  day  the  situation  is  reversed.  In  the  northern 

half  a  patch:  of  positive  northward  values  was  observed  with 

a  maximum  value  of  15.1  m  /  see  at  j  =  10.  In  the  southern 

half  of  the  region  a  patch  of  negative  southward  transports 

v/as  observed,  with  -15.3  m  /see  at  j  =  5.  On  the  60th  day 

the  distribution  resembles  the  pattern  observed  on  the  54th 

day.  The  southward  transports  attained  at  the  upper  level  on 

2  2 

this  day  the  maximum  of  -27.3  in  /sec  at  j  =  12  and  the 

2  2 

northward  maximum  of  45 .3  m  /sec  at  $  =  6. 

The  same  pattern  persisted  on  the  65th  day.  The  north- 

2  2 

ward  transports  attained  the  higher  value  of  53.0  m  /sec 

at  j  =  5  and  a  maximum  negative  transport  of  -33.0  m  /sec 

at  0  =  10.  At  the  750  nib  level  the  same  inversed  situation 

2  2 

was  observed.  A  maximum  northward  transport  of  12.6  m  /sec 
was  observed  at  j  =  12  while  in  the  southern  region  a  maximum 

p 

negative  transport  of  -11.3  m  /sec  was  obtained  at  3=5. 

The  mean  zonal  heat  balance  for  the  same  4  representative 
days  is  shown  in  fig.  20. 
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Two  maxima  in  the  sensible  heat  flnx  can  be  distinguished. 
One  in  the  South  wqs  located  . on  the  48th  day  (which  is  a  max¬ 
imum  day)#at  3  =  3.  Then  it  shifted"  southward  to  3  =  2  during 
the  3  additional  representative  days.  The  other  in  the  north 
was  located  systematically  at  j  =  14.  This  may  have  been  a 
result  of  the  fact  that  in  the  present  model  no  clear  cut 
index  cycle  was  observed. 

The  net  heating  and  cooling  distribution  in  the  present 
model  has  the  general  features  of  the  low-index  situation  in 
the  previous  run,  i.e.,  only  two  subregions  are  observed: 
heating  taking  place  between  the  southern  boundary  and  3  =  ^0 
for  the  48th  day,  between  the  southern  boundary  and  3=9 
for  the  54th  day,  between  the  southern  boundary  and  3  = 
for  the  60th  and  the  65th  day.  Thus  the  boundary  between  the 
two  subregions  shifts  somewhat  as  time  goes  on.  In  the  present 
numerical  experiment,  maximum  heating  was  observed  on  the 
48th  day,  at  3=2,  amounting  to  18.7  *  102  k3  m"*2  sec**1. 

WO  m" 

The  most  intense  cooling  on  the  same  day  was  -21.5x10  k3  m  sec 
near  the  pole,  at  3  =  14.  This  corresponded  to  observations 
which  indicate  a  pronounced  sink  near  the  pole.  On  the  54th 

day  the  maximum  of  the  heating  attained  the  value  of  2l|XlO“2k3 
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m  sec  at  3=2  in  the  south  and  the  maximum  value  of 
cooling  of  -23.5  *  10 2  k3  m  2  sec  1  at  3  =  14  in  the  north. 

On  the  60th  day  the  maximum  appeared  at  the  same  latitudes, 
of  the  corresponding  values  being  26.5x10~2k3  nf^sec”1  and 
-23X1 02  kj  m”2  sec”1. 

On  the  65th  day  the  maximum  values  were  observed  at  the 
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same  latitudes  attaining  the  values  of  22.2x10  k3  m  sec 
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and  -21 .5x10""^k j  m  ^sec-1 

Fig.  21  shows  the  march  of  the  various  kinetic  energy 
components  with  time. 

The  variation  of  the  energy  of  the  zonal  flow  at  the  250mb 
level  did  not  vary  periodically  as  in  the  previous  run  hut  was 
characterized  "by  a  slov/  hut  steady  decrease.  The  only  local 
maximum  was  observed  on  the  28th  day.  After  this  day,  there 
was  a  slow,  hut  continuous  descent  until  the  42th  day.  After 
this  day,  the  decrease  became  more  intense  until  the  60th  day 

on  which  a  tendency  towards  stabilization  began  to  be  manifested. 

2  -2 

Starting  from  5 33  m  sec  the  zonal  flows  energy  decreased  to 
250  m  /sec  '  on  the  oOth  day  stabilizing  itself  at  this  value 
until  the  end  of  the  run  on  the  65th  day.  From  the  53**d  day 
onwards,  the  eddy  kinetic  energy  at  the  250  mb  level  became 
greater  than  the  kinetic  energy  of  the  zonal  flow  at  the  3ame 
level . 

At  the  750  mb  level  the  variations  were  less  clear  cut. 
the  K750  comP°nervfc  had  a  slight  maximum  on  the  28th  day, 
coinciding  in  time  with  a  slight  maximum  in  ^250*  Another 
local  maximum  was  observed  on  the  42nd  day,  and  thereafter 
K ytjQ  decreased  slowly  until  the  end  of  the  run.  The  energy 
of  the  perturbations  exceeded  that  of  the  zonal  flow  at 
750  mb  from  the  48th  day  and  it  increased  slowly  up  to  the 
59th  day  where  it  appeared  to  stabilize  itself  around  the 
value  of  67.5  m  sec  . 
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The  mean  total  kinetic  energy  |  K  +  K*  ^  attained  a 
small  minimum  on  the  22nd  day  and  *•  rnnximum  on  the 
28th  day.  After  a  small  decrease  it  “began  to  rise  from  the 
32nd  day  and  grew  steadily.  Attaining  a  maximum  value  of 
395.5m  sec  ^  on  the  58th  day,  v/hen  a  tendency  towards 
stabilization  could  be  observed,  the  values  of  the  mean 
total  kinetic  energy  oscillated  around  the  value  of  387.5 
nf'sec  until  the  end  of  the  run  on  the  65th  day. 

Compared  with  the  previous  rum  the  march  of  the  various 
components  of  the  kinetic  energy  is  smoother  with  less  oscil¬ 
lations  and  shows  a  tendency  towards  stabilization  at  the 
end  of  the  run. 

l-’ig.  22  shows  us  the  transformation  of  zonal  potential 
energy  into  eddy  potential  energy  and  the  frictional  losses. 

No  cycle  appears  in  the  variation  of  [p.P'J  with  time.  Its 
intensity  was  much  lower  compared  to  the  previous  run,  at¬ 
taining  similar  values  only  towards  the  end  of  the  run.  A 
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first  maximum  of  2.2  X  10  kj  m  sec  was  observed 
on  the  14th  day.  A  second  maximum  was  observed  on  the  27th  day, 
with  the  value  of  0.45  X  10  ^  kj  m  ^  sec”1. 

After  a  minimum  on  the  29th  day  continuous  growth  was 

—It  —2  —1 

observed,  leading  to  a  maximum  of  0.95  X  10  kj  m  sec  on 
the  53rd.  day.  From  this  day  on  the  transformation  decreased 

until  the  end  of  the  run,  reaching  on  the  last  day,  the  65th  day, 
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the  value  of  0.73  X  10  kj  m  sec  .  The  frictional  losses 
of  the  eddies  were  quite  small.  They  did  not  exceed  10  kirn  sec 
until  the  47th  day  of  the  integration.  From  this  day  on 
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a  gradual  increase  of  the  fricticnul  losses  of  the  eddies 

_  k  _p  _■* 

was  observed,  reaching  a  maximum  value  of  27X10  kj  m  sec 

on  the  56th  day.  After  that, a  gradual  decrease  followed  so 

—5  —2  —2 

that  on  the  6r>th  day  the  value  of  2.3  X  10  lc3  m  sec 
was  obtained. 

The  frictional  losses  of  the  zonal  flow  dominate  those  of 
the  eddies  until  the  47th  day.  From  that  day  on  a  gradual 
decrease  was  observed  until  the  end  of  the  run.  two  maxima 
of  the  [F,K]  transformation  on  the  15th  and  26th  day  lag 
one  day  behind  the  corresponding  maxima  of  the  [P,P'J  trans¬ 
formation.  On  the  27th  day  the  transformation  reached  the 

value  of  1 .0  X  10  ^  lcj  m“2  sec”1  and  on  the  47th  day  1 .1  X 
-5  -2  -1 

10  kj  m  sec  .  At  the  end  of  the  run  the  intensity  of 

_p 

the  transformation  became  very  small:  0.2  X  10  kj  m  sec  . 

The  area  averaged  precipi table  water  as  a  function  of  time  is 
shown  in  fig.  23-  The  broken  line  corresponds  to  a  constant 
value  of  15  mm  which  is  the  average  w inter  amount  quoted  by 
Peixoto  (1956).  The  water  content  rose  from  the  initial  stage 
until  it  reached  a  maximum  of  16.3  mm  on  the  23rd  day.  After  this 
day,  a  tendency  towurd  lower  values  appeared,  with  much  smaller 
oscillations  than  in  the  previous  model.  No  index  cycle  was 
observed  and  the  slight  oscillations  superimposed  on  the 
general  trend  had  variable  periods.  On  the  49th  day  the  water 
content  reached  the  value  of  the  winter  average  given  by  Peixoto 
and  continued  to  decrease  until  the  end  of  the  run,  when  on 
the  65th  day  the  water  reached  a  value  of  13»7  mm. 
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The  area  averaged  precipitation  and  evaporation  as  a 
function  of  time  is  shown  in  fig.  2h.  At  the  start  of  the 
integration  -  no  evaporation  was  observed, because  of  the 
vanishing  surface  wind  speed  at  1000  mb.  The  evaporation 
grew  vapidly  -  while  the  precipitation  started  only  on  the 
10th  day.  Until  the  21st  day  precipitation  was  less  than 
evaporation.  Prom  the  21st  day  on  the  evaporation  and  the 
precipitation  had  a  similar  distribution, their  values  being 
close  one  to  the  other;  the  precipitation's  values  slightly 
exceeding  the  evaporation  values.  A  first  maximum  of  the 
precipitation  was  observed  on  the  16th  day,  with  0.5  nm/day. 

A  second  maximum  of  the  precipitation  was  observed  on  the 
27th  day,  with  1.23  mnv/day.  Precipitation  was  slightly  higher 
than  evaporation  in  the  vicinity  of  the  maxima  while  at  the 
vicinity  of  tha  minima  tho  situation  wuo  rovorcod.  iter  a 
minimum  of  the  precipitation  on  the  30th  clay,  with  0.8s  r-;:/day, 
the  precipitation  values  rose  steadily  with  small  oscillations 
superimposed  on  the  general  trend  until  they  reached  the  highest 
value  of  the  entire  run  on  the  55th  day  with  2. Si  mm/day.  This 
value  is  lower  than  the  maximum  obtained  in  the  previous  run, 
amounting  to  3. 56  iraiv/day.  The  evaporation  in  the  present  model 
on  the  same  day  was  2.77  ran/day.  After  the  55th  '-n.y  a  trend 
towards  lower  values  was  observed,  again  accompanied  by  small 
oscillations  on  the  65th  day.  The  precipitation  reached  the 
value  of  2.75  mm/day  and  the  evaporation  the  value  of  2.61  nun/day 

The  latitudinal  distribution  of  the  relative  humidity  as 
a  function  of  time  is  shown  in  fig.  25. 
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during  the  first  14  days  the  relative  humidity  increased 
constantly,  over  the  entire  region.  On  the  ISth  day  some 
general  I’eatuivs  which  were  retained  during  the  whole  course 
of  the  ran ,  became  already  evident,  The  extreme  south  "beyond 
j  =  2\  was  very  humid,  the  hui.iidity  exceeding  9C,j.  ITorth  of 
j  =  3  a  Veit  of  relatively  dry  air  was  located  in  which  the 
values  occasionally  went  down  to  70^  and  even  less.  The 
smallest  relative  humidity  values  of  64v->  were  obtained  on  the 
41st  day.  A  northward  shift  of  the  dry  belt  such  as  in  the 
previous  run  could  hardly  be  detected,  ho  irregularities?  in 
the  ?C  humidity  belt  were  observed  in  the  present  model, 

\7':.:ie  they  were  quite  pronounced  in  the  previous  run,  part¬ 
icularly  north  of  the  dry  belt.  The  dry  belt  remained  in  the 
southern  half  of  the  region,  and  did  not  penetrate  beyond  the 
central  latitude. 

The  changes  in  the  relative  humidity  within  the  dry 
belt  were  quite  pronounced  during  the  early  stages  of  the  run, 
but  no  indications  of  an  index  cycle  were  observed.  Another 
feature  worth  mentioning  is  the  appearance  of  a  small  dry 
region  in  the  extreme  north  with  values  of  80#.  This  was  ob¬ 
served  during  the  43rd  and  the  48th  days,  at  j  =  14'2‘*  The  region 
in  excess  of  90,j  in  the  northern  half  was  larger  and  more 
compact,  in  the  present  model,  with  a  high  degree  of  uniformity. 
In  the  extreme  south  the  line  of  90^  values  was  Quite  uniform, 
with  few  fluctuations.  Towards  the  end  of  the  run  the  humidity 
was  generally  higher  than  in  the  previous  model,  the  values 
in  the  last  10  days  of  the  run  (after  the  56th  day)  being 
everywhere  higher  than  80;u. 
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In  fig.  27  the  frictional  losses  of  kinetic  energy  due 
to  lateral  eddy  viscosity  are  shown. 

Comparing  the  results  in  the  present  experiment  (llo.  2) 
with  our  first  experiment  sore  striking  features  are  evident: 
the  frictional  losses  of  eddy  kinetic  energy  through 
lateral  eddy  viscosity  attain  much  higher  vaLues  in 
the  course  of  the  run. 

In  tliis  experiment  harmonic  components  with  high  wave 
numbers  had  considerably  hifjier  amplitudes  and  thus  the  trans¬ 
formation  O  ,  aJ  was  more  intense.  On  the  other  hand,  the 
frictional  losses  of  zonal  kinetic  energy  followed  the  same 
pattern  as  in  the  first  experiment  but  during  the  last  20 
days  of  the  run  (i.e.,  from  the  45th  to  the  65th  day) 
considerably  lower  values  of  Jk,  aJ  were  obtained. 

The  ["iC1,  aJ  transformation  followed  as  mentioned  above  a 
different  evolution  in  the  present  run  exhibiting  a  totally 
different  march  in  the  last  half  of  the  present  run. 

Starting  slowly  from  zero  on  the  first  day,  j^K',  aJ 
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attained  a  first  maximum  of  46X10  ltj  m  sec  on  the  loth  day. 

A  second  maximum  was  observed  on  the  27th  day,  58X1 0~7  lcj  m~2sec~] 
much  higher  than  the  corresponding;  value  of  the  same  tranformaticn 
on  the  same  day  in  the  first  experiment.  After  some  small  oscil¬ 
lations  which  lasted  until  the  37th  day,  the  values  of  jjC* ,  iTj 
rose  constantly  and  very  rapidly  to  attain  231X10~7  lcj  m“2  sec"1 
on  the  59th  day  and  from  then  until  the  end  of  the  run  on  the 
65th  oay  its  values  oscillated  around  an  average  value 
of  215X10“7  kj  m“2  sec-1 . 
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The  frictional  losses  of  zonal  kinetic  energy  due  to 

lateral  eddy  viscosity  £  IT,  a]  started  from  a  value  of 
-7  -2  -1 

127  X  10  kj  rn  sec  on  the  first  day  of  the  run  and 

decreased  steadily  until  the  15th  day  when  the  transformation 
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attained  a  value  of  3^*10  kj  m  sec  .  A  local  maximum 
was  obtained  on  the  19th  day,  a  second  one  on  the  29th  day  and 
a  third  one  on  the  43 ru  day,  the  corresponding  values  being 
41X10~7,  45X10-7  and  49*10“7  kj  rrf2  sec"1. 

From  tiie  43rd  day  until  the  end  of  the  run  a  gradual 

decrease  is  observed,  the  value  of  IX  Aj  declining  to  12% 
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10  kj  m  sec  on  the  64th  day. 

The  zonally  averaged  poleward  eddy  flux  of  water  vapor 
i 8  shown  in  fig.  26. 

On  the  40th  day  of  the  run  the  values  were  very  close 
to  those  observed  in  the  real  atmosphere  by  Peixoto  (1958). 
On  the  other  4  representative  days  -  due  to  the  absence  of  an 
index  cycle  -  a  similar  pattern  was  observed,  with  values 
almost  double  of  Peixoto’ 3  values. 

The  largest  poleward  flux  was  located  at  the  low  levels, 
the  maxima  being  a  little  south  of  the  central  latitude.  On 
the  40th  day,  the  maximum  was  5.27 gm.m.kg  '  see"1  obtained 
at  j  =  6.  On  the  48th  day  the  maximum  was  9.66gnum.kg  sec" 
at  j  =  7.  On  the  59th  day  the  maximum  was  11. 60  at  j  «  5,  on 
the  60th  day  it  was  11.04  at  j=5,  and  on  the  65th  day  it  was 
of  10.0  at  j=5.  l/e  notice  that  towards  the  end  of  the  run 
there  is  q  tendency  of  the  poleward  fluxes  to  reach  maxima  in 
south. 
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Fluxes  with  inverse  sign  were  located  on  the  northern 
border  of  the  region,  and  they  appear  on  each  representative 
day.  The  maxima  at  the  upper  level  (the  250  mb  level)  were 
smaller  by  one  order  of  magnitude  than  those  located  at  the 
lower  level,  but  they  appear  consistently  at  the  same 


latitudes 
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INTRODUCTION  OP  THE  MATSUNO  SCHEME 

Introduction 

Equationo  of  a  two  level  quasi -hydrostatic  and  quasi- 
geostropliic  model  were  integrated  again,  this  time  using  the 
Matsuno  (Euler-Backward)  two  level  time  difference  scheme 
instead  of  the  centered  time  difference  scheme  (the  leap¬ 
frog  scheme)  used  in  the  first  two  experiments  discussed 
above,  one  without  and  one  with  the  introduction  of  lateral 
eddy  viscosity. 

This  scheme  has  been  shown  to  have  a  damping  effect  on 

high  frequency  waves  (low-pass  filter)  but  it  is  only 
conditionally  stable. 

Matsuno  (1966)  mentions  that  his  scheme  reduces  the  amp¬ 
litudes  of  longitudinal  horizontal  waves  in  addition  to  those 
of  gravity  waves.  As  gravity  waves  were  eliminated  from  our 
model  by  the  geostrophic  assumption,  only  the  effect 
of  the  scheme  on  the  longitudinal  horizontal  waves  v/ill  be 
observed. 

Experiment  No.  3,  was  similar  in  all  details  to  Ex¬ 
periment  No.  1,  mentioned  above,  without  any  lateral  friction, 
but  with  time-differencing  according  to  the  Matsuno  scheme 
instead  of  the  leap-frog  scheme.  The  integration  v/as  carried 
for  62  days. 


The  calculated  distribution  of  the  field  variables  was 
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compared  with  the  results  of  2 x.  So.  1  by  Muss  (1967)  mentioned 
above  and  also  with  the  two  experiments  already  discussed  in 
this  report. 

Smoothing  of  the  field  variables  was  observed  in  this 
experiment  mainly  during  the  last  stages  of  the  integration. 

The  first  40  days  of  the  62  days  integration  showed  n'  serious 
departure  from  the  results  of  Muss's  Ho.  1  Exp.  (1967). 

The  Matsuno  scheme  is  as  follows:  (Katsuno,  1966): 

¥/e  are  given  the  partial  differential  equations  : 

(i)  ff~=  W  (i.j-A  8....N) 


The  F's  involving  linite  differencing  operators,  then  Matsuno's 
two  level  scheme  may  be  written: 


(a) 

Ob) 


fa*"-u'*=  F(u)  M 


v/here  '  U***  *  is  an  intermediate  value  in  the  calculations. 
Let  us  assume  that  equations  (l)  are  of  a  pure  oscillatory 

where  Oi  is  the 
frequency  of  the  i-th  component  of  normal  modes.  Let  us 
consider  a  single  component,  and  omit  the  subscript. 


type,  i.e.  icJitli 


Defining 


SI  =  CO  Ab 


the  amplification  factor  for  the 
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f.'atsuno  scheme  will  be: 

A  =  {  in  -  Cl 


and: 


Thus  the  scheme  is  conditionally  stable.  LL  ^  will  be  damped  if 
|fl|  1  ,  neutral  if  |xi|  =  1  and  unstable  if 


M  =- 1 .  To  obtain  the  value  of  fl  at  which  strongest 


damping  occurs,  we  differentiate 


A 


zjL-  -  e.-fi  o 

3  J2  ^  2(7  -  sia  *  SL^~  r  hSL*  +  Sl*y*  u 


Thus  CJ  .j.--— 
cm  ^ 


In  a  meteorological  problem,  such  as  the  problem  under  dis¬ 
cussion,  the  modes  of  highest  frequencies  are  normally  assumed 
to  be  undesireable  noise  effects,  generated  by  various  processes. 
Any  filtering  device,  introduced  into  the  model  either  explicitly, 
or  implicitly,  should,  therefore  have  the  property  of  smoothing 
out  the  higher  frequencies,  and  leaving  the  lower  frequencies 
intact.  Damping y if  present,  should  therefore  increase  monotonously 
with  frequency.  This  would  not  occur  in  the  '.'.at sun o  scheme  if 
J1  is  allowed  to  increase  beyond  —  .  Therefore,  the 

scheme  will  be  stable,  and  will  also  act  like  a  high-cut  filter, 

if  the  time  step  will  be  chosen  so  that:  /\t  -=^~.  j _ 

nJT1  cjmcul) 


CJ, 


being  the  highest  frequency  in  the  system. 


5 

-] 


< 
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It  is  easily  seen  that  the  scheme  acts  a?  a  filter  bp 

(e+O* 

eliminating  LL  from  the  equations  defining  the 

scheme . 


f*+i)  .  /  / 5>  .  j  \  /*>  .  2/65/2 

u.  -  u  u,  icju.  At  =  6cJu.  At  -  ^  d.  At  S 


By  comparing  the  resulting  expression  with  the  Taylor  ex- 

2  2 

pansion  we  note  that  the  scheme  introduces  a  term. 


which,  in  a  way,  acts  like  a  diffusion  term  in  time. 


The  Development  of  the  Circulation 

In  fig.  29  the  variations  of  the  amplitudes  cf  the  first 
harmonic  wave  components  of  the  stream  function  at  the  central 
latitude,  at  250  nib  and  at  750  mb;as  a  function  of  time  are 
shown.  Curing  the  first  30  days  of  the  integration,  no  serious 
departure  from  the  results  of  Exp.  1  of  Kuss  was  observed. 

The  amplitude  of  the  first  component  at  both  reference 

levels  reached  its  first  maximum  on  the  14th  day,  amounting  to 
7  2 

1  .5  X  10  m  /sec  at  250  mb.  The  evolution  from  now  on  until 
the  30th  day  was  similar  to  that  observed  in  the  Experiment 
described  by  Hu.ss,  with  a  time-lag  of  two  days. 

On  the  32nd  day,  the  amplitudes  at  the  upper  level  of 

7  2 

reference  were  2.7  X  10  m  /sec,  as  compared  with  the 

7  2  ?  o 

3.4*  10  m  /sec  in  Huss’s  experiment  ana  2.5  X  10'  m  /sec 

in  the  experiment  with  eddy-viscosity. 
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During  the  latest  stages  of  the  run  no  oscillations  were 

observed.  The  dominating  feature  was  a  smooth  and  gradual  rise 

of  the  amplitudes  at  both  reference  levels  reaching  on  the 

7  2 

50th  day  values  of  5.28  X  10  m/sec  at  the  250  mb  level 
7  2 

and  2.20  X  10  m  /sec  at  the  750  mb  reference  level. 

From  the  50th  day  until  the  end  of  the  run,  on  the  62nd 
day,  the  amplitudes  remained  about  equal  to  the  values  of  the 
50th  day.  These  stable  values  were  twice  as  high  as  those  ob¬ 
served  in  the  experiment  including  lateral  eddy  viscosity,  i.e. 
2.65  X  107  m^/sec  at  the  250  mb  level  and  6.4  X  10^  m^/sec 
at  the  750  mb  level  (stabilization  taking  place  already  on 
the  42nd  day  of  the  run,  with  eddy  viscosity).  This  appears 
to  be  due  to  the  different  smoothing  involved  in  the  two  models. 

Pis*  (30)  shows  us  the  amplitudes  of  the  odd  harmonic  com¬ 
ponents  of  the  disturbance  at  the  central  latitude  of  the  upper 
reference  levels.  The  component  with  wave  nurriber  (I)  is 
dominant  throughout  the  entire  run  -  as  in  the  earlier  ex¬ 
periments.  A  prominent  feature  of  the  present  run  is  that  the 
magnitudes  of  the  components  with  even  wave  numbers  were  extremely 
small.  Thus  the  noise  due  to  non-physical  effects  is  eliminated 
by  the  use  of  the  Hatsuno  scheme  leaving  the  interactions  bet¬ 
ween  the  odd  wave  number  components  which  are  mainly  of  physical 
nature. 

The  3rd  component  exceeds  10^  m  /sec  on  the  6th  day,  the 
5th  on  the  11th  day  and  the  7th  on  the  19th  day.  Thus  in 
the  beginning  this  run  exhibited  the  same  features  as  those 
which  were  characteristic  of  Hubs' s  experiment  I. 


5 

The  3rd  component  attained  a  first,  maximum  of  9*5  X  10 

2 

m  /eec  on  the  17th  day,  and  after  a  series  cf  small  oscillations 
stabilized  itself  from  the  48th  day  onwards  at  a  value  of 
2.5  X  107  m2/ sec . 

6  2 

The  5th  component  attained  a  first  maximum  of  10  m  vsec 
on  the  21st  day  and  after  a  series  of  oscillations  began  to 

K  O 

stabilize  at  the  value  of  3*5  X  10^  m  /sec. 

U  2 

The  7th  component  reached  a  maximum  of  6.5  X  10  m  /sec  on 
the  22nd  day,  then  underwent  a  series  of  large-amplitude  oscil¬ 
lations.  Towards  the  end  of  the  run  the  values  oscillated 

C  r\  1 1  O 

between  10^  m  /  sec  and  2.0  \  IQ4,  m/sec,  The  amplitude 

of  these  oscillations  decreased  somewhat  towards  the  end  of 

the  run.  The  amplitudes  of  the  components  were  much  larger, 

generally  two  or  even  three  times  larger  than  these  observed 

in  the  experiment  with  lateral  viscosity.  Fig.  31  shows  us 

the  amplitudes  of  the  various  harmonic  components  of  the 

disturbance  with  even  wave  numbers  at  the  central  latitude 

of  the  upper  reference  level.  Here,  a3  in  the  experiment  with 

lateral  eddy-viscosity,  the  amplitudes  never  exceeded  a  value 

j,  2 

of  10  m  /sec  and  their  magnitudes  were  the  smallest  in  all 
the  three  experiments.  On  the  42nd  day,  for  instance,  the 

2  2  p  p 

second  component  attained  7x10  m/sec^the  4th  4X10  m /sec 
2  2 

and  the  6th  10  m  /sec.  On  the  last  day  of  the  present 

run,  the  S2nd  day,  the  2nd  component  attained  a  value  of  2.8x10^ 
2  2  2 

m  /see,  the  4th  4X10  m  /sec  and  the  6th  component  attained 
2  2 

7  X  10  ra  /sec.  This  drastic  damping  of  the  even  harmonics 
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which  is  more  severe  than  in  the  corresponding  experiment 
with  eddy-viscosity,  is  reflected  "by  the  high  degree  of  sym¬ 
metry  "between  the  northern  and  southern  halves  of  the  field 
variables  around  the  central  latitude,  which  was  conserved 
perfectly  until  the  end  of  the  run. 

The  variation  of  the  zonally  averaged  zonal  flow  is 
shown  in  fig.  32  as  a  function  of  latitude  and  pressure  for 
h  selected  days.  . 

During  the  earlier  stages  of  the  run,  i.e.,  until  the 
30th  day  the  patterns  were  quite  similar  to  those  obtained  in 
Husa’s  experiment  Ho.  1,  including  a  tendency  toward  the 
formation  of  a  split  jet.  Thus  on  the  10th  day  the  highest 
velocities  at  the  upper  reference  level  37.3  iVsec  were  located 
at  j  =  11  and  j  =  5»  while  at  the  central  latitude  a  vaLue 
of  35.7  m/sec  was  observed.  On  the  20th  day  no  trace  of  the 
splitting  of  the  jet  was  observed  and  the  highest  velocity 
within  the  single  well-developed  central  jet  was  of  44  ir/sec. 

On  the  same  day  easterlies  were  observed  at  the  lower  boundary. 
On  the  30th  day  the  central  jet  already  had  a  maximum  value 
of  48.6  m/sec  at  the  250  mb  level.  The  easterlies  extended 
now  beyond  the  750  mb  level  penetrating  up  to  j=  5  in  the 
south  and  up  to  j  =  1 1  in  the  north.  On  the  40tL  day  the 
central  jet  had  become  slightly  weaker  with  a  maximum  velocity 
of  45  m/sec. 

The  corresponding  value  of  the  maximum  obtained  in  IIuss's 
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experiment  I,  was  of  38*7  m/sec.  The  easterlies  had  a  similar 
extent  as  those  of  the  30th  day. 

On  the  50th  day  further  decrease  in  the  velocities  was 
observed  with  a  very  slight  tendency  towards  splitting  of  the 
jet,  the  velocity  at  the  central  latitude  at  250  mb  being 
30  Vsec,  while  at  j  =  9  and  j  =  7  it  was  30.2  m/sec. 

The  extent  of  the  easterlies  was  smaller,  but  they  still 
penetrated  a  little  beyond  the  750  mb  level. 

On  the  60th  day  of  the  run,  the  weakening  of  the  velocities 
became  even  more  pronounced,  as  we 11  as  the  splitting  of  the 
jet.  The  velocity  at  the  central  latitude  at  250  mb  was  24.2  m/sec 
v/hile  the  maxima  of  26.6  m/sec  were  located  at  j  =  10  and 
.1=6. 

The  easterlies  were  found  only  below’  the  750  mb  level, 
end  their  lateral  extent  was  also  reduced:  they  did  not  penetrate 
inward  beyond  j  s  12  and  j  =  4. 

The  variation  of  the  zonally  averaged  flow  at  1000  mb  as 
a  function  of  latitude  and  time  is  shown  in  fig.  34.  As  could  be 
inferred  already  from  fig.  32,  the  development  followed  almost 
the  same  pattern  as  the  corresponding  distribution  of  Huss's 
experiment  Ho.  1  until  the  30th  day  with  slight  changes  and 
a  time  lag  of  2  days. 

The  highest  velocity  was  obtained  at  the  central  latitude 
on  the  17th  day  amounting  to  G.4  n/8ec  as  compared  to  the 
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corresponding  value  of  9.8  n^/cec  observed  in  Hubs’ 8  ex¬ 
periment  on  the  15th  day. 

The  boundaries  between  the  central  region  and  the  region 
of  the  easterlies  were  smoother  in  the  present  experiment  and 
especially  during  the  last  stages  of  the  run  (days  40-61 ) 
the  boundary  line  has  fev/  sinuosities. 

The  boundary  line  was  located  between  j  =  5  and  j  =  7 
in  the  south  and  between  d  =  9  and  j  =  11  in  the  north. 

After  a  first  maximum  of  the  surface  wind  speed  at  the 
central  latitude  on  the  17th  day,  a  minimum  of  5.0  n^/sec  was 
observed  on  the  20th  day,  followed  by  another  maximum  of 
8.2  nv/sec  on  the  23rd  day.  This  in  turn  was  followed  by  a 
gradual  decrease  down  to  3.3  m/sec  on  the  30th  day.  The  next 
maximum  of  5.7  ir/sec  was  observed  on  the  43rd  day. 

Thereafter,  until  the  end  of  the  run  a  continuous  decrease 
of  the  surface  wind  speed  was  observed,  the  westerlies  at  the 
central  latitude  decreasing  down  to  2.8  m/sec  on  the  62nd  day. 

Easterlies  with  values  exceeding  2.5  n\/sec  were  observed 
between  the  16th  -  58th  day,  except  for  the  26th  day. 

The  highest  velocity  attained  by  the  easterlies,  4.7  ny'sec, 
was  observed  on  the  43rd  day  simultaneous  v/ith  a  maximum  in  the 
westerlies . 

Local  maxima  in  the  easterlies  were  observed  on  the  17th 
and  23rd  days  -  at  j  a  3  in  the  south  -  also  coinciding  in 
time  with  maxima  of  the  westerlies  indicating  a  strong  correlation 
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between  the  intensities  of  the  easterlies  and  westerlies. 

The  latitude-pressure  distribution  of  the  poleward  eddy 

•j 

momentum  flux  is  shown  in  fig.  33  for  4  characteristic  days. 

On  the  10th  day  the  transport  was  towards  the  central 
latitude  at  both  levels  and  resembled  the  pattern  obtained 
by  IIuss  for  the  same  day. 

On  the  20th  day  trends  in  opposite  directions  appeared  at 
the  lower  reference  level. 

The  maximum  values  of  the  fluxes  directed  towards  the 

2  -2 

central  latitude  amounted  to  23.0  m  sec  ,  An  additional  local 

2  —2 

maximum  of  20.4  m  sec  was  observed  further  tov/ards  the  center. 
On  the  30th  day  the  values  observed  had  increased  and  the 
maximum  magnitude  was  40.8  m  sec  .  Trends  in  opposite  direc¬ 
tions  were  located  in  the  extreme  South  and  North  with  the  flux 
directed  toward  the  boundaries. 

On  the  40th  day  the  values  dropped,  the  maximum  observed 
2  —2 

being  28.6  m  sec  .  At  the  lower  reference  level  positive 
transports  extended  over  the  whole  northern  half  of  the 

2  — p 

region  with  a  maximum  of  4.3  m  sec  . 

Symmetrical  negative  transports  appeared  in  the  south. 

A  totally  different  distribution  was  observed  on  the 
50th  day.  At  the  upper  level  in  the  north  small  negative 
transports  were  observed, maxima  having  dropped  to  the  smallest 
value  observed  during  the  entire  run,  -8.3  m2/sec2  - 
while  at  the  lower  reference  level  positive  transports  attained 
values  of  11.8  m2/Sec2  at  j  s  11. 
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The  positive  values  in  the  North  extended  "beyond  the 

500  mb  level.  On  the  60th  day  a  similar  distribution  was 

observed.  At  high  latitudes  in  the  north,  momentum  transports 

reached  -25  m'/sec  .  At  the  lower  reference  level  positive 

transports  were  observed  in  the  North  and  negative  transports 

in  the  South  with  ±  9*9  m  /sec  .  Their  extent  Y/as  more  limited, 

and  these  fluxes  remained  below  the  500  mb  level.  The  change 

with  time  of  the  various  components  of  kinetic  energy  is  shown 

in  fig.  56.  The  zonal  kinetic  energy  at  the  upper  level  decreased 

from  an  initial  value  of  536  rn  /sec  to  a  value  of  414  m  /sec^ 

on  the  35th  day.  Then  the  values  dropped  more  abruptly  to  reach 

264  m  /sec"  on  the  52nd  day  and  appeared  to  become  stable  at  a 
2  2 

value  of  248  m  /sec  for  the  last  days  of  the  run. 

At  the  lower  level,  ^750  firs-t  increased  slightly  until  the 
25th  day  and  then  a  gradual  constant  decrease  was  observed 
until  the  end  of  the  run. 

A  reversed  situation  was  observed  in  the  case  of  the  eddy 

kinetic  energy.  At  the  upper  reference  level,  %■' 250  incl’e»sed 

2  2 

from  zero  to  a  maximum  of  156  m  /sec  on  the  21st  day.  After 
that  a  sharp  decrease  down  to  60  m  /sec  on  the  25th  day  was 
observed.  From  this  day  on  K'2ro  increased  rapidly.  On  the 
47th  day  the  kinetic  energy  of  the  eddies  became  equal  to 
the  kinetic  energy  of  the  mnal  flow.  A  stable  value  of 
378  m  /sec  was  attained  on  the  56th  day.  The  kinetic  energy 
of  the  eddies  at  750  mb  increased  beyond  the  kinetic  energy 
of  the  zonal  flow  at  this  level  on  the  44th  day. 
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The  total  kinetic  energy  K  +  K'  increased  during  the 

formative  stages  of  the  run  from  280  m  /sec  on  the  first  day 
2  2 

to  322  m  /sec  on  the  21st  day  of  the  run,  decreasing  down 
to  269  m  /sec  on  the  29th  day.  After  this  day  continuous 
growth  of  the  total  kinetic  energy  up  to  the  52nd  day  was 
observed.  On  this  date  a  value  of  345  m  /sec  was  observed  with 
little  further  change  observed  until  the  conclusion  of  the  run. 
Various  energy  transformations  are  shown  in  fig.  35. 

The  transform  {p.p-}  attained  its  first  marked 

-5  -2 

maximum  of  7.0  X  10  kj  m  sec  on  the  21st  day  and  another 

—5  —2  —1 

pronounced  maximum  of  7.40  X  lO^kjm  sec  towards  the  end 
of  the  run,  on  the  47th  day.  The  values  fluctuated  during  the 
first  half  of  the  run,  dropping  occasionally  down  to  1.4x  10"”^ 
k j  m“^  sec”^ .  During  the  second  half  of  the  run  the  transformation 
increased  up  to  the  maximum  of  the  47th  day  mentioned  above  and 
then  a  gradual  decrease  with  a  tendency  towards  stabilization 
near  the  end  of  the  run,  was  observed.  The  conversion  of  heat 

into  zonal  potential  energy.  (  "Q,  P"  }•  was  nositive,  its  variations 

>  3 1 

being  much  smoother  in  the  present  run  than  in  Huss's  experiment. 

It  attained  its  maximum  on  the  50th  day  amounting  to  6  X  10-^ 

—2  —1 

lc;j  m  sec  thereafter  decreasing  gradually  until  the  end  of 
the  run.  Conversion  of  heat  into  eddy  potential  energy  ,  P'  | 
was  negative  throughout,  its  magnitudes  being  much  smaller  than 
those  of  .  The  largest  magnitude  -  3.25  X  10~J  kj  m~2sec'’1 

was  observed  on  the  51st  day. 

The  transformation  of  zonal  and  eddy  kinetic  energy  into 
friction  are  shown  in  Fig.  35. 
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After  the  42nd  day  {iC' ,  P’}  exceeded  P  },the  first 
increasing  and  the  second  decreasing  in  time  from  the  "begin¬ 
ning  of  the  run,  corresponding  to  the  "behaviour  of  &25O  an<^ 

K*  .  towards  the  end  of  the  run,  however,  each  of  those 
2!)U 

transformations  seemed  to  stabilize  at  a  constant  value. 

The  zonally  averaged  diabatic  heating  as  a  function  of 
latitude  is  shown  in  fig.  37  for  6  representative  days.  On 
the  10th  day  very  small  heating  was  observed.  On  the  20th 
day  the  region  of  most  intense  heating  was  located  in  the 
south  at  j  =  2  with  8x10“2kj  m“2sec"1  .  Prom  j  =  6  up  to 
the  central  latitude  cooling  was  observed.  Due  to  the  nature 
of  the  model  no  diabatic  heating  could  result  at  the  central 
latitude.  Small  positive  values  were  found  between  3*8  and 
j  =  10  and  cooling  is  observed  north  of  j  =  10  with  a  max¬ 
imum  at  j  =  14. 

On  the  30th  day  the  value  of  maximum  heating  dropped  to 

—2  — p  _-t 

6.4X10  Icjm  sec  and  the  zone  of  negative  values  in 

the  South  extended  from  j  =  4z  to  j  =  8.  On  the  40th  day 

-2  -2  -1 

the  heating  became  more  intense  with  maximum  of  12X10  !cj  m  sec 
On  the  50th  day  the  most  extreme  values  of  cooling  and  heating 
during  the  entire  run  were  observed  amounting  to  1162X10  m 
sec  .  The  zone  of  negative  values  in  the  south  disappeared 
so  that  the  entire  southern  half  was  a  region  of  heating 
while  the  entire  northern  half  was  a  region  of  cooling. 

The  same  pattern  persists  until  the  end  of  the  run,  the 
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value  of  the  maximum  heating  in  the  south  at  j  =  2  dropping 

“P  —1 

to  14.1  X  1C;  Icj  m  sec  on  the  60th  day.  We  notice  that 
in  all  field  variables  a  tendency  of  attenuation  and  smoothing 
towards  the  end  of  the  run  can  be  observed. 

In  figs.  38,  39,  40  the  zonally averaged  meridional 
velocities  Vd  at  250  mb  in  cm  sec  and  thei r  variation 
v/ith  latitude  and  time  is  shown.  The  mean  meridional  flov/ 
was  obtained  from  the  computed  values  of  the  vertical  velocity 
at  500  mb.  The  values  of  the  lower  information  level  are 
equal  in  magnitude  but  opposite  in  sign  to  those  of  the  250  mb 
upper  level. 

Fig.  38  presents  the  values  obtained  in  Exp.  !\o.  1  with 
eddy-viscosity.  Until  the  6th  day  a  single  Perrel  cell  \ras 
observed; from  this  day  on,  the  familiar  distribution  of  3 
cells  in  vertical  planes  emerged, wi th  a  central  Ferrel  cell, 
surrounded  by  2  direct  Hadley  cells. 

On  the  15th  day  values  of  -29  cm  sec  were  observed  in 
the  center  of  the  Perrel  cell  at  the  central  latitude.  This 
distribution  persisted  until  the  end  of  the  run  on  the  42nd 
day.  The  first  maximum  was  obtained  on  the  15th  day,  a  second 
maximum  was  obtained  on  the  21st  day,  the  values  obtained  at 
the  central  Perrel  cell  being  -35  cm  sec  and  the  third  and 
last  maximum  was  observed  on  the  36th  day  v/ith  -41  cny'sec  at 
the  central  latitude. 

Tv/o  points  are  perhaps  worth  being  mentioned.  First  -  the 
3  cellular  character  persisted  almost  during  the  entire  run, 
except  during  the  first  6  formative  days. 
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Second  -  the  maximum  mean  meridional  velocity  values 

obtained  in  Exp.  No.  1  were  only  about  half  those  obtained 

previously  by  Huss  in  1965  when  no  lateral  viscosity  was 

introduced.  n,he  positive  values  in  the  Hadley  cells  did  not 

*•  i 

exceed  ;  cm  sec  at  most. 

Fig.  39  corresponds  to  Exp.  No.  2,  including  also  eddy- 
viscosity  in  which  water-vaoor  was  introduced  as  an  additional 

j 

.1'  i  el  d  va  ri  abl  e . 

Again  until  the  7th  day  only  the  indirect  Ferrel  cell 
appeared; occupying  all  the  region.  In  this  experiment  no 
symmetry  around  the  central  latitude  was  expected  once  pre¬ 
cipitation  began. 

As  in  Exp.  I  the  familiar  3  cellular  pattern  persisted 
until  the  end  of  the  run.  A  first  maximum  in  the  central  Ferrel 
cell  was  obtained  on  the  15th  day  at  the  central  latitude, 
with  -33  cm/sec,  a  second  maximum  on  the  26th  day,  with  -60  cn^/sec, 
and  a  third  on  the  48th  day  with  -60  cm/sec.  This  was  also  ob¬ 
served  on  the  51st  day  at  the  same  latitude.  From  the  51st  day 
on,  the  values  of  the  mean  meridional  velocity  decreased,  drop¬ 
ping  to  a  minima  of  -37  cm/sec  on  the  67th  day,  at  the  central 
latitude  (in  the  Ferrel  cell).  Until  the  40th  day  the  Southern 
Hadley  cell  was  larger  than  the  Northern  one,  the  highest 
values  rising  up  to  10  cm/ sec  at  3=3.  From  the  40th  day  on, 
the  northern  Hadley  cell  became  more  prominent,  with  a  maximum 
velocity  of  +21  cn\/sec  on  the  48th  day.  This  was  followed  by 
a  decline  until  the  end  of  the  run  on  the  65th  day.  While  the 
values  observed  in  this  run  we re  higher  than  the  corresponding 
values  in  Exp.  1,  no  marked  asymmetry  developed. 
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Pig.  No. 40  shows  the  pattern  obtained  in  Experiment 
No.  3>  with  the  Iilatsuno  time  differencing  scheme.  In  this 
experiment  a  striking  feature  emerges  immediately,  the 
persistence  of  the  3  cellular  structure  from  the  beginning 
of  the  run  until  the  end  on  the  62nd  day.  A  first  maxima 
of  -40  cn/sec  was  observed  on  the  15th  day  at  1=7  and 
j  =  9  in  the  central  Perrel  cell.  Another  maximum,  the 
highest  attained  in  the  run,  v/as  observed  on  the  22nd  day, 
with  a  value  of  -65  cn/sec  at  the  central  latitude. 

Another  local  maximum  was  observed  on  the  42nd  day  with 
-44  cn/sec  at  the  central  latitude  and  the  following,  last 
one,  was  observed  on  the  52nd  day,  with  -43  err/ sec  at  the 
central  latitude.  Prom  this  day  until  the  end  of  the  run  the 
magnitudes  of  the  mean  zonal  meridional  velocity  in  the  central 
Perrel  cell  decreased,  reaching  on  the  62nd  day  -32  cm/sec  at 
the  central  latitude. 

The  positive  values  in  the  Kadley  cells  never  exceeded 
17  cn/sec,  a  value  attained  on  the  36th  day.  It  appears  thus 
that  the  use  of  the  Matsuno  scheme  results  in  a  meridional 
circulation  which  most  closely  resembles  the  meridional 
distribution  in  the  atmosphere,  for  all  three  experiments 
described  here. 

In  figs  41,  42,  the  time  averaged  mean  geostrophic 
poleward  flux  of  sensible  heat  at  500,  for  exp.  No.  1 
and  Exp.  No.  2  is  shown. 
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10  —1—1 

The  units  ere  in  10  Icj  sec  cb 

In  fig.  43  the  observed  time  .averaged  mean  geostrophic 
polev/ard  flux  of  sensible  heat  at  300  mb  is  given,  as  obtained  by 
Mintz  (1955) from  daily  circumpolar  synoptic  charts  (in  the 
Northern  Hemisphere,  in  winter). 

As  the  size  of  the  grid  interval  used  in  our  model  is 
444  lan,  it  ro ug'nly  corresponds  to  4  degrees  at  45  °N. 

Bearing  this  in  mind,  as  well  as  the  limitations  of  the 

model,  the  two  figures  can  be  compared  .  In  Exp.  I,  the  two 

maxima  of  the  sensible  heat  flux  were  located  at  j  =  10  and 

10  ~1  — 1 

3=6  respectively,  with  values  of  3»3  X  10  kj  sec  cb 

The  values  decreased  symmetrically  towards  the  south  and 

10  —1—1 

the  north  down  to  0,2o  X  10  kj  sec  cb  at  j  =  14  and 

3  =  2.  At  the  central  latitude,  i.e.,  at  j  =  8  a  value  of 
10  -1-1 

2.9  X,  10  kj  sec  cb  was  observed. 

Tig.  42  shows  the  corresponding  distribution  obtained  in 
Exp.  Ho.  2.  The  mean  geostrophic  polev/ard  flux  of  sensible 
heat  at  500  mb,  was  time  averaged  between  the  30th  and  the  60th 
days  of  the  integration. 

Due  to  the  introduction  of  water  vapor  no  symmetry 

around  the  central  latitude  was  observed.  A  single  maximum 

10  —1—1 

was  observed  at  j  =  9  of  4.5  x  10  kj  sec  cb  ,  higher  by 

50,j  than  the  corresponding  maximum  in  Exp.  No.  1  . 


Values  decreased  more  sharply  towards  the  north  than  towards 
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the  south,  down  to  a  value  of  0.21  X  1010  lcj  sec”1  cb”1  at 

10  _-f 

0  =  11+  and  of  0.27  X  10  k3  sec  1  ch  at  3  =  2. 

Fig*  43  shows  us  the  time-averaged  mean  poleward  eddy 
flux  of  sensible  heat  at  500  mb  as  observed  by  Mintz. 

This  is  a  wintertime  average  (Jan-Feb  1949).  Y/e  notice 
a  single  maximum  of  3*41  X  1010  kj  sec”1  cb  1  at  55°  IT. 

The  decrease  toward  the  north  is  not  more  rapid  than  the 
decrease  towards  the  south. 

In  broad  outline  quite  a  good  correspondence  with  Hintz's 
observational  values  was  obtained.  It  appears  that  in  Exp.  No. II 
the  heat  flux  was  somewhat  overestimated. 

Exp.  I  gave  realistic  magnitudes  but,  due  to  the  nature  of 
the  model,  presented  two  symmetrical  maxima. 

In  fig.  44  the  kinetic  energy  losses  due  to  the  use  of 
Matsuno's  time  differencing  scheme  are  shown. 

These  losses  were  computed  at  the  end  of  each  day  of  the 
integration  between  the  23rd  and  24th  time  steps  following  the 
formula: 

(K+  K')/Ctl,-(K  -(K,F),5)At-fK  +  KV  =  Mi 

Until  the  15th  day  the  scheme  is  taking  energy  out  from  the 
model,  but  during  the  1 6-1 8th  days  and  also  during  the  period 
of  the  20th  to  the  26th  day  the  model  is  gaining  energy  from 
the  scheme's  application  .  The  values  were  fluctuating  between 
0.3  “  0.6  m  /sec  per  day.  Towards  the  end  of  the  run  this  rate 
decreases  and  negative  values  are  again  occasionally  observed. 

It  is  to  be  stressed  that  In  order  to  gain  further  insight  -  a 
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much  more  detailed  analysis  is  to  be  effectuated. 

Yet,  even  this  simple  analysis  shows  that  Matsuno's  scheme 
has  t  as  a  rule,  a  damping  effect,  and  on  the  whole  is  causing 
kinetic  energy  losses. 

As  an  illustration  we  have  computed  daily  changes  of  the 
total  kinetic  energy  and  received  values,  the  magnitude  of 
which  was  of  the  same  order  as  the  values  above  multiplied 
by  the  number  of  time  steps  per  day. 
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CONCLUSIONS 

Two  snoothing  techniques  were  applied  in  our  model. 

The  results  lead  us  to  the  conclusion  that  the  Matsuno 
scheme  (Euler-backward  scheme)  appears  to  give  "better-  results. 

The  second  order  diffusion  term,  which  was  first  included 
in  the  model  was  very  effective  in  damping  waves . 

The  damping  of  the  short  wave  lengths  prevented  accumulation 
of  energy  at  those  scales  -  thus  stabilizing  the  solution  - 
and  preventing  non  linear  instability. 

Yet,  the  results  seem  to  indicate  that  larger  wave  lengths 
were  also  damped.  Thus,  in  the  first  experiment,  no  significant 
development  was  observed  after  the  40th  day,  while  in  the  second 
experiment  (which  included  a  humidity  field)  -  a  steady  state 
was  attained  after  the  60th  day.  The  choice  of  a  constant  coef¬ 
ficient  of  lateral  diffusion  (A  =  10^  cm2sec*”^  =  10^  m^sec*"^  in 
our' case)  which  is  quite  common  in  the  numerical  experiments 
simulating  the  general  circulation,  is  derived  from  some  large- 
scale  turbulence  theories  by  Richardson,  Batchelor  and  Kolmo- 
g or off . 

Theory  and  observation  indicate  that  the  lateral  diffusion 
coefficient  has  a  strong  dependence  on  the  vertical  coordinate, 
and  changes  also  with  latitude  and  season.  As  pointed  out  by  R.P. 
Kurgatroyd  (1969)  these  values  can  range  from  2X10 ^  -  20x10^m2sec”1 
depending  on  height,  latitude  and  season. 
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Ti'erefore  a  more  refined  and  differentiated  smoothing 
would  he  obtained  if  different  coefficients  of  diffusion 
v/ould  be  employed  at  the  different  levels  of  the  model. 

our  results,  derived  from  the  use  of  the  second  order 
diffusion  ten.!  of  the  Fickian  type  -  just  open  the  question 
of  the  practical  upper  limit  of  any  long-term  numerical 
integration  -  which  uses  diffusion  or  diffusion-like  terms 
for  smoothi  ng. 

An  upper-limit  of  several  weeks  (up  to  SO  days)  imposes 
itself, as  it  seems  that  afterwards  no  significant  changes  can 
occur. 

In  our  third  experiment  the  Fatsuno  time  differencing 
scheme  v/as  used,  acting  implicitly  as  a  smoothing  device. 

It  introduced  a  diffusion -like  second  term  in  time,  but  it 
had  a  more  selective  damping  action  -  filtering  mainly  the  high 
frequency  waves  -  provided  At  is  chosen  so  that  the  scheme 
will  remain  computationally  stable. 

As  pointed  by  Y.  Kurihara  ( "1 965 )  this  gives  a  relation 
between  A  t  and  the  shortest  wavelength  we  can  treat. 

An  even  more  highly  selective  damping  would  presumably 
result  with  the  use  of  a  3“level  time-differencing  scheme  as 
pointed  out  by  Aurihara  (1965).  This  scheme  is  known  under 
the  name  of  the  Fodified  liuler  Backward  Iteration. 

In  the  third  experiment,  the  results  after  a  62  days 
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integration  sho\7  again  a  typical  smoothing  effect  -  and  the 
development  reached  a  steady  state  at  the  end  of  the  integration 
period. 

However  in  this  experiment,  a  more  realistic  distribution 
of  the  field  variables  was  obtained  -  and  they  had  a  higher 
degree  of  symmetry  around  the  central  latitude,  which  was 
conserved  until  the  final  stages  of  the  integration. 

Another  feature  was f that  even-wave  numbers  components  of 

the  disturbance  at  the  central  latitude  were  highly  damped. 

Our  experiments  point  out  that  the  numerical  value  for  the 

5  2  -1 

eddy  diffusion  coefficient,  10-'  m  sec  ,  was  somewhat  too 
large.  A  confirmation  of  our  results  can  be  seen  in  A. 
Grammeltvelt’ s  (19^9)  results  for  the  primitive  equations.  He 
shows  that  using  the  Euler  Backward  time-differencing  method 
for  his  model  gives  a  decrease  in  available  energy  which 
corresponds  to  the  use  of  the  leap  frog  time  differencing  scheme 
with  an  eddy  viscosity  coefficient  of  K  =  3.6  X  10^  m2/sec. 

The  introduction  of  a  more  refined  vertical  resolution 
(i.e.,  with  more  vertical  levels)  together  v/ith  more  refined 
smoothing  devices  like  variable  eddy -diffusion  coefficients  - 
or  more  sophisticated  time -differencing  schemes  -  may  promote 
more  extended  and  physically  significant  long-term  numerical 
integrations. 
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Fig.  1.  Amplitudes  of  the  first  harmonics  of  the  stream 
function  at  the  central  latitude,  at  250  mb  and 
750  mb.  In  m/sec. 
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Fig.  2.  Amplitudes  of  the  odd  wave  numbers  harmonics 

of  the  stream  function  at?the  central  latitude, 
at  the  250  mb  level.  In  m  /sec. 
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Fig.  3.  Amplitudes  of  the  even  wave  number  harmonics 

of  the  stream  function  at2the  central  latitude, 
at  the  250  mb  level  in  m/sec. 
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VARIATION  OF  U4  WITH  LATITUDE  AND  TIME 
(EXPERIMENT  I) 
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Fig.  5.  The  zonal ly  averaged  flow  at  1000  mb  as  a  function  of 
latitude  and  time. 
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Fig.  6.  The  latitude  pressure  distribution  of  the  zonally  averaged 
northward  eddy  transport  of  momentum  in  ro/sec^  on  the  10th 
20th,  30th  and  40th  days  of  the  integration. 


PRESSURE  PRESSURE 


Fig.  7.  Area  and  pressure  averaged  energy  transformations  as  a  function 
of  time. 


Fig 


8  The  area  averaged  transformation  of  frictional  losses  of  kinetic 
energy  bv  lateral  eddy  viscosity  as  a  function  of  time  in 
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area  AVERAGED  kinetic  ENERGY  per  unit  mass  (INmJW?) 
(EXPERIMENT  I) 


Fig.  9.  Area  averaged  values  of  the  kinetic  energy  of  the 
zonal  flow,2the  (jddies  and  of  the  total  kinetic 
energy  in  rn  /sec  as  a  function  of  time. 
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Fig.  10.  Zonally  averaged^influx  of.sensible  heat  into  unit 

column,  in  kj  vT*  sec  xlO^  on  the  10th,  20th,  30th  and 
40th  days  of  the  integration. 


Fig.  11.  Frictional  losses  of  kinetic  energy  as  a  function 
of  time. 
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Fig.  12.  Amplitudes  of  the  first  harmonics  of  the  stream 

functions  at  the  central  latitude  at  250  mb  and  750  mb 
2 

In  m  /sec. 


Fig.  13.  Amplitudes  of  the  odd  wave  numbers 

harmonics  of  the  stream  function  a£  the  central 
latitude,  at  the  250  mb  level  in  rn  /sec. 


Fig.  14.  Amplitu  '  s  of  the  odd  wave  numbers  harmonics  of 
the  st  ;am  function  at  the  central  latitude  at 
the  250  mb  level.  In  m  /sec.  Exp.  Ill,  1967. 


Pig.  15.  Amplitudes  of  the  even  wave  numbers  harmonics  of  the 

stream  function  at  the  central  latitude  at  250  mb  level 
In  m  /sec. 
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Fig.  16.  Amplitudes  of^the  even  wave  numbers  harmonics  of  the  stream 
function  in  m  /sec  (Exp.  III.  1967). 
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Fig.  17.  .he  zonallv  averaged  flow  at  1000  mb  as  a  function  of  latitude 
and  time. 
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POLEWARD  MOMENTUM  flux  (IN  rrflitc*) 
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Fig.  19.  The  latitude  pressure  distribution  of  the  zonally 
averaged  poleward  eddy  transport  of  momentum,  in 
m2/sec2  on  the  48th,  54th,  60th  and  65th  days  of 
the  integration. 


Fig.  20.  Zonally  averaged  heat  balance  components  per 

unit  column  on  the  48th,  54th,  60th  and  65th  days  of 
the  integration. 


AREA  AVERAGEO  KINETIC  ENERGV  PER  UNIT  MASS  (IN  m1  set' ) 
(EXPERIMENT  II) 


TIME,  IN  DAYS 


Area  averaged  values  of  the  kinetic  energy  of  the 
zonal  flow,tl}e  edcjies,  and  of  the  total  kinetic 
energy,  in  m  /sec  as  a  function  of  time. 
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Fig.  22.  The  area  averaged  transformation  of  zonal  potential 

energy  into  eddy  potential  energy  and  of  the  frictional 
losses,  as  a  function  of  time 


Fig.  23.  Area  averaged  prccipitable  water  content  per  unit 
column  in  runfl,  as  a  function  of  time. 


Fig.  27.  The  area  averaged  transformation  of  frictional  losses 

of  kinetic  energy  by  lateral  eddy  viscosity  as  a  function 

^  i 

of  time  in  40  kjm  "sec 

j 

i 


m*/sec 
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Fig.  30.  Amplitudes  of  add  wave  numbers  harmonics  of  the 

stream  function  at  the  central  latitude  at  250  mb 
level.  In  m  /sec. 


Fig.  31.  Amplitudes  of  the  even  numbers  harmonics  of  tne  stream 
function  at  the  central  latitude  at  250  mb  level. in 
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VARIATION  OF  Clt  WITH  LATITUDE  AND  TIME 
(EXPERIMENT  III) 


Fig.  34.  The  tonally  averaged  flow  at  1000  mb  as  a 
function  of  latitude  and  time. 


■vfltA  AVERAGEO  KINEHC  ENERGY  PER  UNIT  MASS  (IN 
(ExPERiMENI  III) 


TIME,  IN  OAYS 


Rig.  36.  Area  averaged  values  of  the  kinetic  energy  of 
the  zonal  flow,  the  ejddie^  and  of  the  total 
kinetic  energy,  in  ni  /sec  as  a  function  of  time. 


159 


ZONAU.Y  AVERAGED  OIABATIC  HEATING  PER  UNIT  COLUMN 
(EXPERIMENT  III) 


LATITUDE 


Fig.  37.  Zonal ly  averaged  influx  of-sensible  heat  into  unit 
column  in  kjm  sec  x  10  on  the  10th,  20th,  30th, 
40th,  50th  and  60th  days  of  the  integration. 
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Fig.  42.  The  time  averaged  mean  geostrophic  poleward  flux  of 
sensible  heat  at  500  mb.  Exp.  No.  2. 


Fig.  43. 


The  time  averaged  observed  mean  geostrophic  poleward 
flux  of  sensible  heat  at  500  mb  -according  to  Mintz. 
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Fig.  44.  The  kinetic  energy  losses  due  to  the  use  of  Matsuno's  time 
differencing  schemes. 
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AN  INTEGRATION  OF  A  THREE-LEVEL  QUASI -GEOSTROPHIC 
MODEL  OF  THE  ATMOSPHERE  WITH  CONSTANT  AND 
PRESSURE  DEPENDENT  STATIC  STABILITY 


by 

A.  Huss  and  M.  Segal 


ABSTRACT 


The  evolution  of  general  circulation  patterns 
has  been  studied  by  means  of  a  three- level  quasi-geostro- 
phic  model . 

Two  experiments  were  performed.  In  the  first 
experiment  static  stability  remained  constant  with  height, 
while  in  the  second  it  varied  with  height. 

The  results  were  compared  to  those  obtained  in 
a  two  level  model  experiment,  with  similar  features,  and 
also  with  distributions  obtained  from  observational  data. 

The  distribution  of  the  various  field  variables 
were  found  to  be  similar  to  those  in  the  two-level  model 
experiment,  however  their  magnitudes  were  found  to  be  as 
a  rule  higher  in  the  present  experiment. 

Magnitudes  of  several  variables,  including  wind 
velocities,  were  unrealistically  high  in  the  last  stages 
of  integration  in  the  second  experiment. 
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I.  INTRODUCTION 


General  circulation  patterns  were  studied  by  means  of 
a  3-level  quasi-geostrophic  and  quasi-hydrostatic  model.  The 
features  of  the  model  are  the  same  as  those  adopted  in  the  two- 
level  model  described  by  Huss  (1967).  The  initial  conditions, 
boundary  conditions ,  differencing  scheme  etc.  will  be  described 
in  the  following. 

In  the  two- level  model  the  thermodynamic  energy  equation 
was  only  applied  at  the  S00  mb  level,  while  in  the  present  experi¬ 
ment  the  thermodynamic  energy  equation  was  applied  at  two  levels, 
the  333  mb  and  the  667  mb  levels.  The  application  of  the  thermo¬ 
dynamic  energy  equation  at  two  levels  instead  of  one,  leads  to  a 
better  definition  of  the  vertical  p-velocity  in  the  model- 
atmosphere. 

No  release  of  latent  heat  and  no  lateral  viscosity  were 
introduced  into  the  model. 

The  influx  of  heat  per  unit  column  was  formulated  in  the 
same  manner,  as  adopted  by  Staff  Members  of  Electronic  Computation 
Center,  Japan  Meteorological  Agency  (1965).  According  to  this 
formulation  the  influx  of  heat  per  unit  column  was  made  proportional 
to  the  temperature  difference  at  the  ocean-air  interface  and  to  the 
geostrophic  wind-speed  at  this  level  -  which  was  taken  as  1000  mb. 

Static  stability  is  a  function  of  the  spatial  coordinates  and 
time.  However  in  quasi-geostrophic  models, as  shown  for  instance  by 
Wiin-Nielsen  (1959),  the  static  stability  should  be  assumed  to 
vary  only  in  the  vertical.  Among  the  models  in  which  the  above 
assumption  was  introduced  we  may  mention  those  of  Staff  Members  of 
Electronic  Computation  Center,  Japan  Meteorological  Agency  (1965), 
Murakami (1964)  and  Cressman  (1963). 
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Two  different  experiments  are  described  in  this 
paper.  In  the  first  experiment,  the  static  stabilities  at 
333  mb  and  667  mb  were  assumed  to  be  equal  to  the  500  mb 
static  stability  of  the  standard  atmosphere.  Thus  in  this 
experiment,  the  static  ability  is  lower  at  333  mb  and 
higher  at  667  mb,  than  the  corresponding  values  in  the  standard 
atmosphere.  In  the  second  experiment  the  standard  atmospheric 
stabilities  are  assumed  at  333  mb  and  667  mb# 
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II.  THE  STRUCTURE  OF  THE  MODEL 


The  x,y,  p  system  of  spatial  coordinates  is  adopted. 

The  0,  167,  333,  500,  667,  833  and  1000  mb  levels  are  indicated 
by  the  subscripts  0,1, 2, 3, 4, 5  and  6  respectively.  The  vorticity 
equation,  from  which  the  evolution  of  motion  is  deduced,  was 
applied  at  p^,  p^  and  p^,  while  the  thermodynamic  energy  equation 
was  applied  at  and  p^. 


The  coriolis  parameter  -p  was  taken  as  the  constant 
P  -4  -1 

r  =  10  sec  ,  except  where  its  y  derivative,  s  ,  appeared.  The 

*  1 1  J 

Rossby  parameter,  3  ,  was  ta.ken  as  the  constant  ji  -  1.6  x  10  m" 


The  streamfunctions  at  p^,  p^  and  p^  are  related  to  the 
geopotential ,  0  ,  by  means  of  the  expression 


CD 


Thus  the  quasi-geostrophic  approximation  is  adopted,  with 


Vj  -  IK*  V  f 


C2) 


1 


C3) 


|<;  =  unit  vertical  vector 


L  -1 
sec  . 


V.  =.  V  =  geostrophic  wind. 

a 
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The  thermodynamic  energy  equation  can  be  written: 


with:  UJ  =  ,  the  vertical  p-velocity 

Q  =  the  potential  temperature 
H  =  the  heat  added  to  unit  mass 
in  unit  time. 

Cp  3  the  specific  heat  at  constant 
pressure  =  1004  kjt  ^deg 

The  hydrostatic  equation  is  given  by  means  of  the  streamfunction 
as  : 


‘  VS  -  Sr 

7)p  F  7>P  fp 


(5) 


with  :  R  =  the  gas  constant  =  287  kj.t  Jdeg  1 

The  static  stability,  t>  , which  is  assumed  to  be  a  function  of 
the  pressure  only,  is  given  by: 


)  —  VS  VSnS 

®  ~  df>  OfD  C6) 

The  constants  i=2  and  i  =  4  are  defined  as: 


with 


op  =333  mb. 
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In  terms  of  the  streamfunction  (4)  is  written: 


(8) 


By  applying  equation  (8)  at  p£  and  p^,  and  replacing  the  pressure 

derivatives  by  cent^-d  differences,  the  following  expressions 

for  w,  and  U>,  are  obtained: 

2  4 


The  frictional  stress  is  assumed  to  be  dependent  on  surface  wind  and! 
the  pressure,  according  to  the  expression 


~  QP)  \ 


(10) 


with  Km  *  mean  surface  stress  coefficient 
K  =  geostrophic  surface  wind 
The  pressure  dependence  of  C(p)  is  defined  by: 


C11) 
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The  amount,  of  sensible  heat  entering  unit  column  of  the 
atmosphere  at  1000  mb  from  the  lower  boundary  (assumed  to 
represent  an  ocean  surface)  was  expressed  by: 

%  ■  \  Vs  -  V  <*2) 

with  :  =  the  wind-speed  at  the  surface 

Ts  =  the  temperature  of  the  air,  near  the 
surface. 

=  the  temperature  of  the  underlying 
surface. 

A  =  Constant, 
s 

The  influx  of  sensible  heat  was  assumed  to  become  distributed 
instantaneously  and  homogenously  throughout  the  whole  column  of 
air.  Thus  the  influx  of  sensible  heat  per  cb  is  given  by: 


(13) 


Radiation  losses  were  assumed  to  be  constant  throughout  the  whole 
region.  Thus,  because  of  a  symmetric  chocs e  of  conditions,  the 
structure  of  the  model  equations,  and  by  assuming  that  an 
instantaneous  heat  balance  was  maintained  at  each  instant  of  time 

i 

the  integration,  zero  net  outgoing  radiation  results. 


Applying  the  vorticity  equation  at  the  three  reference 
levels  P^»P3  andpg,  while  using  U)  =  C  as  a  boundary  condition 
for  u>  ,  leads  to: 
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4%\ur-° 


( %  *  %■  1} I K)  +  ^  +3K£%  --C 


Sv+Vj-V)!^*)-*  a*) 


The  wind  velocity  at  and  was  expressed  by  means 
of  the  values  at  the  two  adjacent  information  levels: 

V  -  ^|-mV* 

v  ^  ^irS  (15) 

1  ri¬ 

ll  e  cause  of  the  quasi-geostrophic  assumption  we  obtain: 

V-  -yi] 

C16) 

\V  ^  ii -%) = y.M\ %■%%■%) 

Eliminating  bd,  and  ^  between  (9)  and  (14),  while  using  (16)  leads 
to  the  final  form  of  the  prediction  equation,  in  terms  of  the  poten¬ 
tial  vorticites,  at  the  three  information  levels  P^.Pj  and  p^: 


ur,p>  -< 


H  =  -vr  *U -W  ->A 

(17) 

^  +V,  ^ 

where  the  potential  vorticities,  ,q3  and  q£  are  defined  by  the 
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following  diagnostic  equations: 

h  =u 


(18) 


The  following  averages  are  used  in  various  computations: 

n»xj« 


o 


(19) 


where  x  is  the  length  in  the  lateral  extent  of  the 
rectangular  region  and  y  its  north-south  extent  . 


Let  the  subscript  N  and  S  denote  the  northern  and 
southern  boundaries  respectively.  With  this  notation  the 
boundary  conditions  were  written: 


(20) 


(21) 


In  the  x-direction  cyclical  continuity  of  the  field  variables 
is  assumed. 
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The  boundary  conditions  for  the  potential  vorticity  will  be 
discussed  later. 

For  computing  the  tonally  averaged  meridional  velocity,  the 
following  requirement  is  made: 


where  is  the  tonally  averaged  meridonal  divergent 

component  of  the  wind. 

The  zonally  averaged  meridional  divergent  velocity  is  derived 
from  the  zonally  averaged  equation  of  continuity: 


Sty  _ 

'  -r 

The  hydrostatic  equation  (5)  can  oe  written  as 


(23) 


d 


K  /  /  Jf 
-7--  dCvp 
r 


1 


(24) 


The  application  of  finite  differencing  at  p2  and  p^leads  to: 


Assuming  the  Standard  Atmosphere  lapse  rate  between  667  mb 
and  the  lower  boundary: 


V.  19 


(26) 
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Integrating  (24)  between  p,.  and  p^: 


5L  JL. 
9 '  °'s 


(27) 


This  equation  gives  the  relation  between  the  streamfunction  at  p^ 
and  pg  by  means  of  the  temperature  at  p^. 


The  finite  difference  form  of  t''e  Jacobian  term: 


The  finite  difference  scheme  for  the  advective  terms 
in  the  prognostic  equations  was  the  scheme  suggested  by  A.  Arakawa 
(1966).  By  using  this  scheme  and  appropriate  boundary  conditions, 
not  only  the  area  integral  of  the  Jacobian  vanishes,  but  also  the 
advective  process  described  by  the  finite-difference  of  the 
Jacobian  does  not  leadto  a  fictitious  generation  of  mean  energy  and 
mean  square  potential  vorticity. 

Denoting  by  ^  the  finite  difference  Jacobian, 
and  i  and  j  as  indices  increasing  along  x  and  y  direction,  and  by  J 
the  total  number  of  grid  points  along  the  y  direction,  the  Arakawa 
analogue  of  the  Jacobian  at  inner  points  will  be: 


(28) 
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and  on  the  boundaries  ; 


(29a) 


-(ifi-p-.*  +^i,  -h-, 

+(fc>  -#)&■»] 

i  !f  fV/jJij-  -  ■is?'  +  y*  »j-(  -J  &)  f;„,j 

4  (<l  1  + 

~(K t->  )h^j-D 


where  'f'.  and  are  the  constant  values  of  'f' 
on  the  southern  and  northern  boundaries  respectively,  arid /IS  is 
the  grid  interval. 

The  scheme  for  calculating  the  newvaluesofthepotential  vorticity 
on  the  boundaries  . 


The  Jacobian  term  on  the  boundary  is  given  by  (29a)  and 


(29b). 


For  the  frictional  term,  we  introduce  the  Laplacian  on 
the  boundaries  in  the  following  way: 


i-K0  (a‘Y  Ity'1-1  ~ 'f°) 
"*•  i  r  fe’- I  % 


The  surface  velocities  required  for  the  calculation  of  the 
diabatic  heating  are  defined  on  the  boundaries: 


1VW 
lvkj,6  ~  h 
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The  calculation  of  the  new  values  of  the  potential  vorticity 
on  the  boundaries,  is  done  by  applying  the  above  boundary  condi¬ 
tions  with  equation  (17). 

The  finite  difference  form  of  the  equations 

In  order  to  perform  the  calculations,  equations  (17) 
and  (18)  have  to  be  re-written  in  finite  difference  form. 

Centered  finite  differences  replace  the  time-derivatives 
appearing  in  (17): 


(cty)  ~  u ,  r.-t— 

l  OfcA,j,T  (32) 

The  finite  difference  approximation  for  the  advection 
term  is  given  by  (28). 

The  Laplacian  of  the  streamfunction  will  be  approximated  as: 


vy  %j) 

Thus  the  equation  (17)  will  be  approximated  by: 

kpjpY 'X-.jji.-i  ~  f~  -i  j 


(34) 
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and  the  diagnostic  equations  (18)  by: 


1 ,  z.  —  J  Sj'/TT  “^V  f  K  Zjjj-C’  fc,  vj,  x)  — 

^  S  j  W-  ■  [Aiy  7  fs, l% jjir +A,  t%  >  si,  r  -  :  j',  c) 

The  energy  transformations 


The  notations  and  the  methdds  which  are  used  for 
obtaining  the  energy  transformations  are  the  same  as  descri¬ 
bed  by  Huss  (•j96'iL  for  the  two  level  model  experiment.  However 
in  the  present  experiment  energy  transformations  must  be 
computed  separately  for  the  lower  half  of  the  atmosphere,  and 
the  upper  half  of  the  model  atmosphere,  and  the  corresponding 
values  for  the  two  layers  are  averaged.  Energy  transformations 
are  computed  per  unit  column;  as  follows  : 


■  X 


(36) 
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The  U  equation 


I 


The  diagnostic  equation  for  the  vertical  velocity  tt, 
is  derived  by  eliminating  time- derivatives  from  the  potential - 
vorticity  equations  and  the  thermodynamic  energy  equations: 


Equations  (9)  can  be  written  as  : 

^  Jl'fl .  l/r)  +UV  +  Hj 


(37) 


By  using  (17)  and  (18)  we  obtain  the  following  equations: 


^  -  (fr)J  vfl 


1  82 


Eliminating  ,  2^  ( y;  —  Vi)  ,  VJ(y>  -  vO 

and  'ft)  by  means  of  (37)  we  obtain  the  final  form  of  the 

diagnostic  equations  for  W: 
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III.  EXPERIMENT  I 


The  static  stabilities  at  Pj  and  p^  were  assumed  to  be 
both  equal  to  the  Standard  Atmospheric  value  at  p^.  The  duration 
of  experiment  I,  was  31  days. 

Fig.  1  shows  the  variation  of  the  amplitudes  of  the  first 
harmonics  of  the  streamfunctions  at  the  central  latitude,  at  167,  500 
and  833  mb,  as  a  function  of  time. 

In  the  first  10  days  an  almost  exponential  increase  in 

amplitudes  is  observed.  During  this  period  the  magnitudes  of  the 

amplitudes  at  833  mb  and  at  500  mb  are  almost  equal.  Following  this 

initial  increase,  a  sharp  decrease  in  amplitudes  was  observed  and 

5  2 

on  the  12th  day  the  amplitudes  reached  local  minima  of  57  x  10  m  /sec 

5  2  5  2 

at  833  mb,  95  x  10  m  /sec  at  500  mb  and  157  x  10  m  /sec  at  167  mb. 

After  attaining  these  minima  (the  only  pronounced  minima  of  the 

entire  run)  the  amplitudes  grew  rapidly  until  the  16th  day. 

Thereafter  and  up  to  the  end  of  the  run,  a  general  trend  of  slow 

increase  in  amplitudes  was  observed,  accompanied  by  small  oscillations. 

The  variation  of  the  zonal ly  averaged  zonal  component  of  the 
wind,  at  1000  mb,  as  a  function  of  latitude  and  time  is  shown  in  Fig.  2  . 
Except  for  the  last  days  of  the  run,  the  boundaries  separating 
easterlies  from  westerlies  were  located  in  the  south,  between  j  ■  4  and 
j  =  5,  and  in  the  north,  between  j  =  11  and  j  =  12.  The  maximum  intensi¬ 
ty  of  the  westerlies  was  obtained  always  at  the  central  latitude. 

The  first  maximum,  observed  on  the  12th  day,  had  a  value  of  about 
10  m/sec,  and  was  accompanied  by  low  values  of  the  amplitude  of  the 
disturbances  in  the  streamfunctions.  Other  dominant  maximum-velocities 
were  attained  on  the  18th,  22th  and  30th  days.  The  last  one  had  a 
value  of  15  m/sec  -  the  highest  value  observed  during  the  entire  run. 

The  highest  speed  of  the  easterlies,  -  6  m/sec,  was  observed  on  the 
18th  and  on  the  30th  day  .  Maxima  of  surface  wind  velocities  were 
larger  than  those  observed  in  the  two  levels  model  experiment.  Thus 
frictional  energy  losses,  and  the  influx  of  sensible  heat  from  the 
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surface,  being  dependent  on  surface  wind  velocity,  could  attain 
larger  maximum  values. than  in  the  case  of  the  two  level  model. 

The  zonally  averaged  zonal  flow  as  a  function  of  latitude 
and  pressure  for  3  selected  days  is  shown  in  Fig.  3.  On  the  13th 
day  the  region  of  easterlies  extended  to  j=5  at  the  south  and  to 
j  =  1 1  at  the  north,  and  penetrated  somewhat  beyond  p^.  The  maximum 
values  of  westerlies  were  attained  at  the  central  latitude.  The 
velocity  at  at  the  central  latitude  was  55  m/sec,  while  at  p^ 
and  p5  the  corresponding  velocities  were  26  m/sec  and  13  m/sec  respec¬ 
tively. 


On  the  21st  day  the  distribution  and  intensity  of 
easterlies  remained  almost  the  same  as  on  the  13th  day.  At  the  upper 
levels  a  split"jet"  had  developed  with  two  maxima  north  and  south 
of  the  central  latitude.  Tiie  velocity  at  p^,  at  the  central  latitudes 
was  32  m/sec,  while  the  maximum  velocity  at  this  level,  37  m/sec, 
was  located  between  j=5  and  j  =6  and  between  j =10  and  j=ll.  The 
splitting  of  the  jet  at  p^was  less  pronounced,  while  at  p^  no 
such  feature  was  observed.  Compared  with  the  13th  day,  the  intensity 
of  the  westerlies  had  decreased,  the  decrease  being  most 
pronounced  at  the  central  latitudes. 

On  the  31st  day  the  northern  and  southern  regions  of 
easterlies  were  considerably  smaller  in  lateral  extent  and  height, 
compared  w'th  correpsonding  regions  observed  on  the  21st  day. 
Easterlies  were  also  observed  between  j=7  and  j=9,  penetrating 
somewhat  beyond  p^.  The  speed  of  the  easterlies  was  very  low  and 
did  not  exceed  -1  m/sec.  No  slight  easterlies  of  this  nature  appeared 
on  the  days  preceeding  the  31st  day. 

A  very  pronounced  splitting  of  the  belt  of  the  westerlies, 
with  two  maxima  at  about  j=5  and  j=ll(was  observed  at  all  levels. 

The  velocity  at  p^  at  the  central  latitude  was  22  m/sec,  compared 
with  a  velocity  of  42  m/sec  at  j=5  and  j=ll.  The  maximum  values 
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at  were  24  m/sec,  while  at  the  central  latitude  a  value  of  only 
1  m/sec  was  observed. 

The  mean  poleward  flux  of  momentum  on  the  10th  and  on 
the  20th  days  as  a  function  of  latitude  and  pressure  is  presented 
in  Fig.  4. 

On  the  10th  day  two  well-defined  regions  were  observed 

in  the  northern  and  southern  halves  of  the  model-atmosphere,  with 

positive  flux  in  the  south  and  negative  flux  in  the  north, 

having  identical  magnitudes  at  symmetrical  location  with  respect  to  the 

central  latitude.  At  p.  the  maximum  positive  and  negative  values,  of 
2  2  1 

182  m  /sec  ,were  observed  at  j=5  and  j=ll  respectively,  compared 

2  2  2  2 

with  maximum  values  of  10  m  /sec  at  p^  and  9  m/feec  at  p^.  Thus, 
momentum  was  transported  predominantly  at  the  upper  layers. 

On  the  20th  day  broad  regions  of  momentum  transports  towards  the 
southern  and  northern  boundaries  (the  highest  value  at  p^  being 
163  m^sec2)  were  observed  in  the  south  and  north  respectively.  Posi¬ 
tive  momentum  transports  were  observed  in  a  relatively  narrow  region 
south  side  of  tbe  central  latitude,  while  corresponding  negative 
momentum  transports  were  observed  north  of  it  .  In  these  regions 
transports  were  less  intense.  On  the  20th  day  as  well  as  on  the  10th 
day  most  of  the  momentum  was  transported  in  the  upper  layers . 

Various  energy  transformations  are  shown  in  Fig.  5  . 

The  energy  transformations  and  were  both 

consistently  positive  during  the  entire  run.  The  changes  of  one  of 
these  transformations  were  parallel  to  those  of  the  other  . 

Except  for  the  first  12  days,  the  differences  between  the  values  of 
these  transformations  were  quite  large.  Thus  only  part  of  the  total 
amount  of  zonal  potential  energy  which  was  converted  into  eddy 
potential  energy,  was  transformed  into  kinetic  energies  of  eddies. 

The  highest  values  of  Jpi  p1^  and  amounting 

-5  -2  -1  -1  _c  _  2  _l  -1 

to  43.2  x  10  kjm  sec  cb  and  19.2x10  kjm  sec  cb  respectively 

were  observed  on  the  30th  day.  The  average  values  of  j  jS,  f>J 
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and  ?  I'  *  A  ■'  were  considerably  higher  compared  with  those 

obtained  in  the  two  levels  experiment. 

During  most  stages  of  the  run  the  [1/  transformation 

was  lower  than  the  Iff'1)  and  transformations. 

Negative  values  of  were  observed  several  times 

during  the  run. 

The  conversion  of  heat  into  zonal  potential  energy  was 
consistently  positive,  while  conversion  of  heat  into  eddy  potential 
energy  was  consistently  negative.  The  numerical  values  of  the  first 
transform  were  slightly  larger  than  those  of  the  second  transform.  Thus 
available  potential  energy  appears  to  increase  somewhat  due  to  heating. 

The  frictional  losses  of  zonal  kinetic  energy  were  small,  with 
-5  -2  -1  -1 

a  maximum  of  32x10  kjm  sec  cb  .  The  frictional  losses  of  eddy 
kinetic  energy  increased  with  time  and  reached  a  value  of  11.3x10  ^ 
kj.m  sec  cb  i  at  the  end  of  the  run. 

The  change  with  time  of  the  various  components  of  the  kinetic 
energy  is  shown  in  Fig.  6.  At  each  level  the  kinetic  energy  of  the  eddies 
increased  in  general  with  time.  Normally  a  situation  of  this  nature  was 
accompanied  by  a  decrease  of  the  kinetic  energy  of  the  zonal  flow. 

However,  no  balance  was  observed  anda  steadyincrease  in  the  total  averaged 
kinetic  energy  took  place.  At  the  two  upper  levels  the  kinetic  energy  of 
the  eddies  exceeded  the  zonal  kinetic  energy  from  about  the  18th  day, 
while  at  p,-  this  was  the  case  already  from  the  12th  day  on.  Following 
these  dates  an  increase  in  the  kinetic  energy  of  the  eddies  and  a  decrease 
in  kinetic  energy  of  tue  zonal  flow  was  observed  at  all  levels. 

The  total  averaged  kinetic  energy  had  a  value  of  about  330 

2/  2  2/2 

mysec  during  the  formative  stages  of  the  run,  and  rose  to  1040  m/sec 
at  the  end  of  the  run. 

Fig.  7  represents  the  distribution  of  the  zonally  averaged  diabatic 
heating  per  unit  column,  for  the  13th, 21st  and  31st  days.  Due  to  the 
structure  of  the  equations  and  the  initial  conditions,  the  distribution 
is  antisymmetrical  around,  the  central  latitude.  On  each  of  those  three 
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selected  days,  the  southern  half  of  the  region  gained  heat,  while  the 
northern  half  lost  heat. 

On  the  13th  day  heating  has  a  maximum  value  of 

-2-1 

0.05  kjm  sec  at  j=3,  and  cooling  the  negative  corresponding 
value  at  j=13.  The  magnitude  of  the  heat  influx  was  considerably 
larger  on  the  21st  and  31st  days.  The  maximum  value  on  the  21st 
day  was  0.29  kjm'  sec  ,  and  on  the  31st  day  0.35  kjm  sec'A. 

These  values  are  about  twice  the  maximum  values  obtained  in  the 
two  level  model  experiment. 

The  latitude-pressure  distribution  of  the  poleward 
flux  of  sensible  heat  is  shown  for  3  selected  days  in  Fig.  8. 

Positive  values  predominated.  On  the  13th  day,  the  negative  sensible 
heat  flux  values  were  observed  at  the  upper  level  near  the  southern 
and  northern  boundaries;  however, the  magnitudes  were  very  small. 

On  the  21st  day  the  negative  values  disappeared,  and  a 

very  pronounced  intensification  of  flux  took  place.  A  maximum  value 
2  -1-1 

of  72x10  kjsec  cb  was  observed  at  p.  at  j*8  while  at  p., 

2  -1-1 

two  maxima  of  50x10  kjsec  cb  were  observed  at  j=7  and  j=9. 

Further  intensification  of  sensible  heat  flux  values 

was  observed  at  p,  on  the  31st  day,  and  a  maximum  value  of 
-1  -1  * 

210  kjsec  cb  was  located  at  j=8. 
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The  amplitudes  of  the  first  harmonics  of  the  disturbances 
at  the  central  latitude,  at  the  three  reference  levels,  as  a  function 
of  time,  are  shown  in  Fig.  9  .  Broadly  speaking,  the  variation 
of  the  amplitudes  at  the  three  reference  levels  was  similar. 


During  the  formative  stages,  the  familiar  pattern  observed 
in  Experiment  I  was  repeated  ;  a  steady,  nearly  exponential  growth 
at  the  three  levels,  with  only  a  small  decrease  within  the  first  day. 

The  first  maximum  was  attained  on  the  13th  day,  amounting  to 

c  2  /  5  2  5 

128x10  m/sec  at  833  mb,  193x10  m  /sec  at  500  mb  and  133x10 

2 

m  /sec  at  167  mb.  From  now  on  up  to  the  end  of  the  run  the  amplitudes 
oscillated  .  These  oscillations  at  the  different  pressure  levels  were 
not  in  phase.  However,  the  dates  at  which  extremal  values  appeared 
at  the  different  levels,  did  not  vary  by  more  than  a  day  or  two. 


Low  amplitudes  were  observed  at  the  three  levels  approximately 

on  the  23rd,  31st,  39th, 46th  and  51st  days.  The  lowest  value 

5  2 

at  p5  was  obtained  on  the  45th  day,  amounting  to  52x10  m  /sec;  at  p3 

on  the  46th  day,  amounting  to  70x10^  m^/sec  and  at  p.  on  the  39th  day 

5  2/  1 

amounting  to  77x10  nysec.Dn  the  whole  the  amplitudes  at  the  two  lower 

levels  had  the  same  order  of  magnitude  as  those  observed  in  experiment  I. 

However,  the  average  amplitude  at  the  upper  level  was  smaller 

in  the  present  experiment,  due  to  the  greater  stability  at  this  level. 


Fig.  10  represents  the  mean  zonal  wind  at  1000  mb,  as 
a  function  of  latitude  and  time.  During  the  first  31st  days  the 
distribution  and  intensity  of  surface  wind  were  similar  to  those  observed 
in  experiment  I.  Thereafter,  the  region  of  westerlies  narrowed  and 
disappeared  entirely  at  about  the  36th  day.  Unitl  the  45th  day  easterlies 
prevailed  at  all  latitudes,  having  a  maximum  of  -5  m/sec  at  the  central 
latitude  .  The  familiar  distribution  of  surface  wind  reappeared  at 
about  the  45th  day. 
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The  first  three  maxima  of  westerlies  were  observed  on  the 
16th,  23rd,  and  31st  days,  reaching  values  of  13  m/sec,  8  m/sec  and 
15  m/sec  respectively.  The  intensification  of  the  westerlies  was 
accompanied  by  an  intensification  of  the  easterlies  flanking  them 
and  hy  a  reduction  of  the  amplitude  of  the  disturbance  at  p,_  . 

The  distribution  of  wind  was  similar  to  that  observed  in 
the  two- level  model  experiment  (except  for  the  period  where  easter¬ 
lies  became  dominant  throughout.the  entire  region)  ,  but  wind  speeds 
were  somewhat  larger  in  the  present  experiment. 

The  latitude-pressure  distribution  of  u  for  5  selected 
days,  is  shown  in  fig.  11.  On  the  13th  day  surface  easterlies  were 
observed  in  the  south  and  in  the  north  regions,  penetrating  close 
to  the  lateral  boundary  beyond  p^.  The  highest  velocity  in  the 

easterlies  was  observed  at  j=3  and  j=13,  with  values  of  -5  m/sec.  The 

westerlies  prevailing  over  most  of  the  region,  are  characterized 
by  smaller  vertical  wind  shears  at  the  lower  layers  compared  with 
those  at  the  upper  layers.  The  maximum  velocity  at  Pj  is  47  m/sec, 
abserved  at  the  central  latitude.  A  similar  pattern  was  observed 
on  the  same  day  in  Experiment  X . 

On  the  21st  day  the  distribution  and  intensity  of  the 
easterlies  remained  almost  unchanged.  A  slight  splitting  in  the 
westerlies  was  observed  at  P3,  without  appearing  at  p1  and  p,.. 

The  westerlies  were  almost  equal  to  those  in  the  previously  described 

day.  Compared  with  the  same  day  in  experiment  I,  the  wind  had  a 
similar  distribution  but  somewhat  higher  velocities. 

On  the  31st  day  a  considerable  intensification  and  extension  of 
easterlies  was  observed  ;  laterally  they  were  found  to  prevail 
beyond  j=5  and  j  =  11  and  vertically  they  had  penetrated  somewhat 
beyond  p^.  Maximum  easterly  velocities  were  attained  near  the 
surface  at  j=3  and  j=  13,  with  -8  m/sec. 
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The  westerlies  also  intensified,  reaching  values  of  67  m/sec 
and  30  m/sec  respectively  at  and  at  the  central  latitude. 

On  the  39th  day  easterlies  were  found  to  prevail  throughout 
a  shallow  layer  near  the  surface,  with  an  upward  bulge  penetrating 
vertically  up  to  p4  in  the  central  part  of  the  region. 

The  westerlies  were  split,  with  two'  maxima  of  75  m/sec  and 
14  m/sec  located  at  j=6  and  j  =  ll  at  p^  and  p^  respectively. 

The  vertical  wind  shear  in  the  upper  layers  was  extremely 
large  compared  with  the  shear  observed  in  the  previous  days. 

On  the  53rd  day  no  easterlies  were  any  longer  observed  at 
the  central  latitudes.  They  recovered  their  normal  spatial  distri¬ 
bution,  with  quite  a  small  intensity,  not  exceeding  -3  m/sec. 

Compared  with  the  39th  day  considerable  intensification 
of  westerlies  was  observed,  the  velocity  at  the  central  latitude  at  Pj 
rising  beyond  95  m/sec. 

On  the  whole  the  wind  velocities  at  the  upper  layers  were 
larger  than  those  observed  in  the  two-level  experiment. 

A  three-level  model,  by  its  very  nature,  enables  a  vertical 
resolution  of  the  vertical  shear.  In  the  present  experiment  a  strong 
shear  was  generated  in  the  upper  layer  and  a  weak  shear  in  the  lower 
layer. 


The  latitude-  pressure  distribution  of  the  poleward  eddy 

momentum  transport  is  shown  in  fig.  12  for  a  number  of  selected 

days.  On  the  13th  day  the  northward  flux  was  contained  within  the 

southern  half  of  the  region,  and  the  southward  flux^with  symmetrical 

magnitudes  around  the  central  latitude  in  the  northern  half.  The  most 

intense  momentum  fluxes  of  +_  79  m  /sec  were  observed  at  p^,  j=5 

and  j=ll.  Highest  values  of  momentum  transport  at  p,  were  observed 

2  2  1 

at  j=4  and  j=12,  amounting  to  72  m  /sec  . 
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The  pattern  on  the  31st  day  remained  almost  the 

same  as  that  on  the  13th  day,  but  considerable  amplification 

took  place.  The  largest  transports  positive  and  negative,  were 

located  at  j=5  and  j=ll  at  all  the  levels  with  magnitudes  of 
2  2  2  2  2/  2 

135  myscc  ,  144  m  /sec  and  57  m/sec  at  p^,  p^  and  p^  respecti¬ 
vely. 


On  tbe  53rd  day  the  pattern  became  quite  complex,  with 
strips  of  negative  values  alternating  with  strips  of  positive 
values.  The  maximum  intensities  were  located  in  the  upper  levels. 
Somewhat  similar  patterns  ,with  about  the  same  intensities,  were 
also  obtained  in  the  case  of  the  two  level  model  at  some  stages 
of  the  integration. 

The  energy  transformations  are  presented  in  Fig.  13. 

Compared  with  experiment  I  the  values  of  [  p",  9’)  and  |p’, 

were  generally  lower.  The  magnitudes  of  these  two  transformations 
during  most  stages  of  the  integrations  differed 
while  the  pattern  of  their  variation  with  time  was  similar. 

Low  values  of  the  transformations  were  accompanied  by  low 
values  of  stream  function  disturbances. 

During  the  first  23  days  the  values  of  ^  ?,  f 'f  and  |  f 
are  quite  low  .  During  ;,the  following  days  a  steady  increase  both 
in  ^p,p'|  and  ^p1,  was  observed  .  On  the  27th  day  the 

maxima  of  20. 4x  10~^  kjm  ^sec  cf and  15.2  x  10'^kjm  ^sec’^di- 
respectively  were  observed ,  followed  by  a  decrease  to  quite  low 
values.  However,  toward  the  end  if  the  run  the  transformations 
increased  again  attaining  their  highest  values  during  the  entire 
run  on  the  58th  day:  30  x  10  ^kjm  ?--"and  21x10  ^kjm  “500 
respectively.  Negative  values  of  both  transforms  were  observed 
between  the  35th  and  37th  days. 
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The  I  transformation  as  i-n  Experiment  I,  was  occasionally 

negative.  Towards  the  end  of  the  run  large-amplitude  fluctuations  of 
this  transformation  were  observed.  Generally  speaking,  the  magnitudes 
of  the  transformations  were  quite  similar  to  the  corresponding  values 
in  the  two  level  model. 

The  conversion  of  heat  into  zonal  potential  energy  was 
consistently  positive,  while  the  conversion  of  heat  into  eddy  poten¬ 
tial  energy  was  almost  always  negative,  except  that  during  some 
stages  of  the  run  this  transformation  vanished. 

Except  for  the  first  20  days,  the  values  of  the  first 
transformation  were  considerably  larger  than  those  obtained  for  the 
second  one.  In  the  two-level  model,  and  in  Experiment  I,  both  trans¬ 
formations  had  almost  the  same  absolute  magnitudes. 

The  maximum  value  of  { <$,  ff  was  observed  on  the 
29th  day,  amounting  to  19.2  x  10'  kjm  secVfe  On  the  same  day  p'J 
had  its  minimum  value^  amounting  to  -3.8x10"^  kjm"^  sec"*£^'! 

It  will  be  noticed  that,  as  has  been  observed  also  in  the  two 
previous  experiments,  maximum  values  of  f  were  generally 

accompanied  with  minimum  values  of  j 

The  values  were  quite  small;  moreover,  during 

some  stages,  this  transformation  became  negative.  Such  negative 
values  occurred  during  the  days  in  which  surface  easterlies  were 
observed  at  all  latitudes. 

The  frictional  losses  of  eddy  kinetic  energy  did  not  exceed 
-5  -2  -i 

7.2  x  10  kjm  sec  .Compared  with  the  corresponding  period  in 
Experiment  I  their  values  were  found  to  be  lower,  but  still  larger  than 
those  observed  with  the  two  level  model. 

The  march  of  the  various  kinetic  energy  components  with  time  is 
shown  in  Fig.  14.  K  values  decreased  slightly 
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during  the  first  21st  days.  Thereafter  a  noticeable  growth 

was  observed,  up  to  2453  in^sec2  on  the  57th  day.  K'167  values 

were  consistently  lower  than  values,  as  was  also 

observed  in  Experiment  1.  Their  magnitudes^ were  lower  than  during 

the  corresponding  period  in  Experiment  I.  On  the  other  hand  K'.,- 

167 

values  were  higher  in  the  present  experiment. 

At  p^  the  kinetic  energy  of  the  eddies  exceeded  during  some 
stages  the  zonal  kinetic  energy  t  however  toward  the  end  of  the  run 
R 500  again  exceeded  .The  maximum  of 

amount  to  160  m%ec2  (the  initial  value  of  was  146  m2/sec2. 

The  maximum  of  k^-  observed  on  the  30th  day,  amounted 

to  219  m/^ec  .  Compared  with  experiment  I  the  average  values 

obtained  in  the  present  experiment,  were  about  the  same,  while  k'toC> 
values  were  considerably  lower. 

Except  for  the  first  19  days,  the  kinetic  energy  of  the 

eddies  at  p5  exceeded  the  zonal  kinetic  energy.  Its  highest  value  of 
2/  2 

120  m/Sec  ,  was  observed  on  the  29th  day,  while  the  maximum  value  of 

2  2 

the  kinetic  energy  zonal  flow  amounting  to  45  m  /sec  ,  was  observed 
on  the  31st  day.  The  average  magnitude  of  both  quantities  was  similar 
to  that  observed  in  Experiment  I. 

The  area  averaged  kinetic  energy,  including  all  components, 
increased  from  the  initial  value  of  328  m^sec2  to  1175  m2/sec2  at  the 
end  of  the  run.  In  the  present  experiment  the  rate  of  increase  in  total 
kinetic  energy  was  larger  than  during  the  corresponding  period  in 
experiment  I.  On  the  31st  day  total  averaged  kinetic  energy  was  almost 
twice  the  value  observed  in  Experiment  I. 

The  zonally  averaged  diabatic.  heating,  .  as  a  function  of  la¬ 
titude,  for  6  selected  days,  is  shown  in  fig.  15.  Except  for  the  31st 

day  the  heating  values  were  quite  moderate,  and  their  maxima  were 
-2  -1  -2  -1 

between  0.07  kjm  sec  to  0.16  kjm  sec  .  Similar  maxima  are  obser¬ 
ved  in  the  two-level  model  experiment.  On  the  31st  day  the  maximum 
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was  0.35  kjm  see  .  Generally,  the  magnitudes  of  both  heating  and 
cooling  were  seen  to  Le  of  die  same  order  of  magnitude  at;  those 
obtained  by  Clapp  (1961)  for  the  winter  period. 

Fig.  In  represents  the  poleward  flux  of  sensible  heat  as  a 

function  of  latitude  and  pressure  for  5  selected  days.  On  the  13th 

day  sensible  heat  flux  had  negative  values  near  the  northern  and 

southern  boundaries;  however  the  absolute  magnitudes  were  small 

compared  with  the  positive  values  around  the  central  latitude,  where 

2  -1  -1 

a  maximum  value  of  42*10  kjm  sec  cb  was  obtained  at  p^ t 
Intensification  and  spreading  of  negative  values  was  observed  on  the 
21st  day,  positive  values  were  dominant  at  the  lower  levels,  while 
negative  values  dominanted  at  the  upper  levels. 

A  similar  distribution  was  observed  on  the  31st  day.  However, 
the  region  with  negative  values  contracted. A  maximum  value  of 

2  i  i 

44x10  kj-sec  cb  was  obtained  at  p.,  at  j=8  .compared  with 

2-1-1  4 
24^10  kj  sec  cb  ,  observed  at  p2  at  j=3.5  and  j=12.5. 

On  the  39th  day  positive  values  prevailed  almost 
thx-oughout  the  whole  regions,  but  the  magnitude  of  the  flux  at 
was  considerably  smaller  than  that  at  p4<  On  the  53th  day  most  values 
were  positive.  Maximum  values  were  obtained  at  p2> 

Generally,  the  variation  of  the  sensible  heat  flux  at  p4  was 
small,  compared  with  pronounced  changes  observed  at  p2. 

Fig.  17  shows  the  time-averaged  zonal  flow  (for  the  38th-60th 
days), Except  for  a  slight  splitting  at  p^  the  westerlies  had  one 
maximum,  at  the  central  latitude.  The  most  intensive  velocity  of 
the  westerlies  was  observed  at  pj.  At  the  central  latitude  the  velocity 
reached  90  m/sec,  compared  with  a  17  m/sec  velocity  at  p^.  Easterlies 
were  observed  at  the  edge  of  the  southern  and  northern  regions,  their 
intensity  being  low,  not  exceeding  -2.2  m/sec. 

The  vertical  shear  of  the  wind  velocity  was  very  pronounced 
in  the  upper  half  of  the  model  atmosphere,  as  compared  with  the  lower 
ha  If. 
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In  the  representation  of  the  zonal  wind  distribution 
in  the  northern  hemisphere,  by  Mintz  (1954), based  on  real  data, 
the  jet  is  situated  around  30°N  in  ^he  winter,  and  its  highest 
intensity  of  40  m/sec  is  found  at  200  mb.  Likewise,  easterlies 
are  observed  at  northern  latitudes  with  maximum  intesity  of 
-5m/sec.  Thus  the  jet  intensity  in  our  experiment  is  about 
twice  that  of  the  real  atmosphere,  while  the  easterlies  are  less 
pronounced. 

Fig.  18  represents  a  time-averaged  eddy  momentum 

transport.  The  values  were  positive  over  the  southern  half 

region,  and  negative  over  the  northern  half.  The  flux  had  extre- 

2  2 

mal  values  of  +  48  m  /sec  at  j  *  5  and  j  *  11,  at  p  .  The 
maximum  values  at  p^mounted  to  8m2/sec2,  and  were  oblerved  at 
the' same  latitudes.  The  momentum  transport  at  p5  was  very 
small  . 

The  poleward  momentum  flux  in  winter  as  estimated 

2/  2 

by  Holopainen  (1967)  has  its  maximum  of  66  m  /sec  at  250  mb 
and  is  located  at  30°N. 

Fig.  19  shows  a  time  average  of  the  poleward  flux  of 
sensible  heat,  as  a  function  of  latitude  and  pressure.  Positive 
values  prevailed  almost  throughout  the  whole  region.  The 
values  at  p2  and  p4  had  similar  magnitudes.  The  sensible  heat 
flux  maxima  were  observed  at  the  central  latitude,  with 
1£  kj-sec  *cb  *  and  15  kj  sec  *  cb  *  at  p2  and  p4  respectively. 
The  values  of  sensible  heat  flux  decreased  gradually  toward 
the  northern  and  southern  boundaries. 
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V  .  CONCLUSIONS 


Two  numerical  experiments  dealing  with  the  evolution 
of  atmospheric  general  circulation  patterns  were  performed  by 
means  of  a  3- level  quasi -geos trophic  and  quasi  hydrostatic  model. 

In  the  first  experiment  static  stability  was  assumed  to  be 
constant  with  height,  while  in  the  second  experiment  it  varied  in  the 
vertical.  These  experiments  were  compared  to  earlier  experiments 
with  a  two  level  model. 

The  resulting  patterns  of  the  various  field  variables 
showed  a  certain  resemblance  to  those  of  the  earlier  experiments; 
however  magnitudes  of  field  variables  were  generally  larger  . 

The  amplitude  of  the  disturbances  in  the  streamfunctions  at  the 
central  latitude  on  the  upper  level  were  smaller  in  Experiment  II, 
than  the  corresponding  amplitudes  in  Experiment  I,  due  to  the 
greater  stability  at  this  level  in  the  second  experiment. 

The  distribution  of  the  zonal  surface  wind  was  generally 
similar  in  both  experiments  to  that  observed  in  the  two  level 
model  experiment.  However,  wind  speeds  were  somewhat  larger  in  the 
present  two  experiments. 

The  time  averaged  zonal  speeds  in  the  westerlies  in  the 
upper  layers  were  considerably  higher  in  Experiment  II  than  values 
observed  in  real  atmosphere. 

The  distribution  of  eddy  fluxes  of  momentum  was  occasionally 
quite  unrealistic,  and  varied  considerably  from  day  to  day.  A  time 
averaged  distribution  for  23  days  ,  in  the  case  of  experiment  II, 
showed  pronounced  negative  eddy  momentum  flux  in  the  northern 
region. 

The  transformations  of  zonal  potential  energy  into  eddy 
potential  energy,  and  from  eddy  potential  energy  to  eddy  kinetic 
energy,  had  a  similar  pattern  of  variation  with  time,  in  both 
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experiments.  Low  values  of  these  transformations  were  accompanied  by  low 
values  of  the  streamfunction  disturbance. 

Due  to  more  intensive  surface  wind  speeds,  transformations 
of  heat  into  potential  energy,  and  frictional  losses  were  higher 
than  those  observed  in  the  earlier  experiment. 

In  both  experiments  a  steady  increase  in  the  total  averaged 
kinetic  energy  was  observed,  due  to  the  f3Ct  that  there  was  no 
balance  between  kinetic  energy  of  zonal  flow,  and  kinetic  energy  of  the 
eddies . 


The  magnitudes  of  the  zonally  averaged  diabatic  heating  were 
generally  higher  in  Kxperiment  I  than  those  observed  in  Experiment  II. 
The  magnitude  of  both  heating  and  cooling  were  seer,  to  be  of  the  same 
order  of  magnitude  as  those  in  the  real  atmosphere  during  winter. 
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AMPLITUDE  OF  FIRST  HARMONIC  OF  STREAMFUNCTIONS 


1.  Amplitudes  of  the  first  harmonics  of  the  stream-functions 

at  the  central  latitude  at  167,  500  and  833  mb. 

2 

In  m  /sec.  Experiment  I. 


VARIATION  OF  Ue  WITH  LATITUDE  AND  TIME 
(EXPERIMENT  I)' 


TIME, IN  DAYS 


The  zonal ly  averaged  flow  ar  1000  mb  as  a  function  of 
Uiiiude  and  time.  Experiment  I 


latituoe 


The  latitude-pressure  distribution  of  the  tonally  averaged 
flow,  UL,  on  the  13th,  21st  and  31st  days  of  the  integration 
Experiment  1 . 
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Fig.  4 


POLEWARD  MOMENTUM  FLUX  (IN  vn2/sec2) 
(EXPERIMENT  1) 


20th  Day 


LATITUDE 


The  latitude-pressure  distribution  of  the  zonally  averaged 
northward  eddy  transport  of  momentum,  in  m%ec2,  on  the  10th 
and  20th  days  of  the  integration.  Experiment  I. 
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Fig.  5.  Area  and  pressure  averaged  energy  transformations,  including 
frictional  losses,  as  a  function  of  time.  Experiment  I. 


Fig.  6.  Area-averaged  values  of  the  kinetic  energy  of  the  zonal^  ^ 
flow,  the  eddies,  and  of  the  total  kinetic  energy,  in  m  /sec  , 
as  a  function  of  time.  Experiment  I. 
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POLEWARD  FLUX  OF  SENSI8LE  HEAT 
(IN  H|  »»c  '  cb  '  xicn  (EXPERIMENT  !) 


Fig.  8.  The  latitude-pressure  distribution  of  the  zonally 

averaged  poleward  flux  of  sensible  heat  in  kj  sec-1cb"^fTO  , 
on  the  13th,  21st,  and  31st  days  of  the  integration.  Expe¬ 
riment  I. 
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Fie.  9.  Amplitudes  of  the  first  harmonics  of  the  stream- functions 
at  the  central  latitude  at  167,  500  and  85^  mb  in  tn  /stv  • 
Experiment  II. 
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Fig.  10.  The  tonally  averaged  flow  at  1000  mb  as  a  function  of 
latitude  and  time.  Experiment  II. 


M£AN  ’OMl  fLoW  iJH  ) 

(EaPLUiMEN'T  Di 

i3U'0a 

:tii  rrn-T  i  i  tit 


LATITUDE 

LATITUDE 

3iih0ay 

39^  Do 

5 _ 

latitude 


0 

U'..7  6 


Fig.  }1.  The  latitude-pressure  distribution  of  the  zonally 

averaged  flow,  ,  on  the  13th,  21st,  31st,  £9th  and 
53rd  days  of  the  integration.  Experiment  II. 
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POltWARO  MOMENTUM  FUJ*  {IN 
{feXPER'MEWTO) 
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Fig.  12.  The  latitude-pressure  distribution  of  the  tonally 

averaged  northward  eddy  transport  of  momentum  in 
2  2 

m  /sec  ,  on  the  13th,  31st  and  53rd  days  of  the  integration. 
Experiment  II 
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Fig.  13.  Area  and  pressure  averaged  energy  transformations, 
including  frictional  losses,  as  a  function  of  time. 
F.xperiment  II. 
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I  arla  AVERAGED  KJNETIC  ENERGV  PER  UNIT 
I  MADS  AT  THE  83Lmb  LEVEL  (IN 

isoL  (Experiment  d) 


Fig.  14. 


Area-averaged  values  of  the  Kinetic  energy  of  the  zonal 

flow,  the  eddies,  and  6f  the  total  kinetic  energy,  in 
2  2 

m  /sec  ,  as  a  function  of  time.  Experiment  II. 
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poleward  elux  or  sensible  heat  (in  x|  »•<'*  cb  xia'i 

(EXPERIMENT  0| 


Fig.  16.  The  latitude-pressure  distribution  of  the  zonally 
averaged  poleward  flux  of  sensible  heat,  in 
kjsec"1cb’1xlO"'‘',  on  the  13th,  21st,  31st,  39th 
and  Sjrd  days  of  the  integration.  Experiment  II. 
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Fig.  17.  Time  averaged  (38th-60th  days)  latitude-pressure 

distribution  of  the  zonally  averaged  flow.  Experiment  II. 

Time  AVERAGED  (Dth-Wthdays)  POLEWARD 


LAT’TUDE 

Fig.  18,  Time  averaged  (38th-60thdays)  latitude-pressure  distribution 

of  the  zonally  averaged  northward  eddy  transport  of  momentum,  in 

m'/sec'-.  Experiment  II. 
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TIME  AVERAGED  (38  th -60  th  DAYS) 
POLEWARD  FLUX  OF  SENSIBLE  HEAT 
(IN  k|  sec  ’  cb-’xlO1).  (EXPERIMENT  n) 


Fig.  19.  Time  averaged  (38th- 60th  days)  latitude-pressure  distribution 

of  the  zonally  averaged  poleward  flux  of  sensible  heat, 

-1-1  -2 

in  kisec  cb  xlO  .  Expenemnt  II. 
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Abstract 

A  three-level  quasi-geostrophic  model  of  the  atmosphere  has  been 
used  to  obtains  three  dimensional  time  dependent  radon  distribution. 

Several  experiments  were  performed  in  which  radon  was  introduced 
at  the  lowest  level.  The  sources  were  distributed  uniformly  -  either  over  the 
whole  region,  or  over  one  half  of  the  total  region.  The  distribution  of  radon 
was  observed  to  respond  quite  realistically  to  the  large  scale  atmospheric 
flow,  and  particularly  to  the  convergence-divergence  patterns,  i.  e.  to  the 
distribution  of  the  vertical  components  of  the  motion. 
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Introduction 

Studies  of  trace  elements  in  the  atmosphere  are  concerned  as 
a  rule  with  two  domains,  that  of  the  large  scale  processes  and  that  of  small 
scale  phenomena,  i.  e.  tracers  have  been  widely  applied  to  the  study  of  prob¬ 
lems  of  stratospheric  and  tropospheric  mixing  (Hunt  and  Manabe,  1968, 

Machta  1966  and  others),  and  to  the  study  of  the  mean  residence  time  of 
aerosols  (for  example  several  papers  reviewed  by  Junge,  1963).  In  these 
studies  synoptic- scale  disturbances  were  treated  as  large-scale  eddies  and 
their  detailed  structure  was  not  taken  into  account.  On  the  other  hand  the 
scavenging  of  tracers  and  their  transport  in  the  boundary  layer  were  treated 
by  several  investigators  (Greenfield,  1957,  Gat  and  Assaf,  1968  and  others), 
under  the  assumption  that  synoptic  features  are  large-scale  phenomena.  Only 
in  a  few  studies  dealing  with  trace  elements  (Danielsen,  1968)  was  the  syn¬ 
optic  scale,  including  its  variations  in  time  and  space, considered  as  the  scale 
of  interest. 

A  scale  analysis  (Assaf,  1969)  assuming  steadiness  and  uniform¬ 
ity,  indicated  that  the  radon  distribution  is  sensitive  to  the  divergence  regime 
and  might  be  used  as  a  m  easure  of  the  vertical  wind  on  the  synoptic  scale.  How¬ 
ever,  these  are  very  strong  assumptions,  if  we  bear  in  mind  that  the  mean 
life  of  radon  is  of  the  same  order  of  magnitude  as  the  period  of  synoptic  oscil¬ 
lations.  Moreover,  changes  in  the  source  distribution  may  become  quite  im¬ 
portant,  particularly  once  we  are  dealing  with  the  effects  due  to  ocean-continent 
contrasts.  The  aim  of  this  paper  is  to  substantiate  this  approach  by  studying 
the  response  of  the  radon  distribution  to  different  source  distributions  and  to 
different  patterns  of  atmospheric  motions. 

It  was  realized  that  it  would  be  easier,  in  the  first  step  in  this 
direction,  to  make  use  of  a  rather  simple  and  crude  numerical  model,  rather 
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than  to  advect  the  radon  by  an  observed  real  distribution  of  winds.  As  the 
simplest  possible  model  adequate  for  our  purpose,  a  three-level  quasi¬ 
hydrostatic  and  quasi-geostrophic  model  was  chosen.  This  model  generated 
a  streamline  field  at  3  levels,  and  vertical  velocities  at  2  levels.  Needless 
to  say,  the  model  does  not  provide  a  very  accurate  analogue  of  the  general 
circulation  of  the  atmosphere,  nor  was  it  intended  to  do  so.  Several  import¬ 
ant  dynamical  processes  were  disregarded,  such  as  small-scale  convection 
and  diffusion.  The  aims  of  the  experiment  were  not  to  obtain  at  this  stage 
realistic  three-dimensional  distributions  of  radon,  but  to  observe  how  radon, 
introduced  at  the  lower  boundary  of  the  model,  would  react  to  a  time-dependent 
circulation  -  and  in  particular,  what  its  schematic  distribution  would  look 
like  in  regions  where  horizontal  divergence  or  convergence  respectively  dom¬ 
inate. 

The  Model 

The  calculation  of  the  velocities  required  to  transport  the  radon 
were  obtained  by  means  of  a  three  level  quasi-hydrostatic  and  quasi-geostro¬ 
phic  general  circulation  model  incorporating  diabatic  heat  sources  and  surface 
friction, whose  details  will  be  described  elsewhere  (see  third  article  in  this 
report)  and  which  resembles  in  many  respects  the  model  constructed  by 
Phillips  (1956).  The  prognostic  equations  were  applied  at  three  pressure 
levels  167,  500,  833  mb,  to  be  denoted  by  indices  1,  3  and  5.  The  predicted 
values  were  the  potential  vorticities  from  which  the  geopotentials  at  the  three 
reference  levels  (proportional  to  the  streamfunctions )  were  derived  by  means 
of  an  inversion  procedure.  From  the  streamfunctions  the  rotational  component 
of  the  horizontal  wind  was  obtained.  The  potential  vorticity  in  the  model  is 
advected  horizontally  only  by  the  rotational  component  of  the  wind,  while  no 
advection  in  the  vertical  is  assumed. 

The  Initial  conditions  correspond  to  a  baroclinically  unstable 
zonal  westerly  flow,  constant  at  each  level,  with  a  small  perturbation  superimposed 
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upon  it.  At  the  lower  boundary  the  velocity  is  assumed  to  vanish.  As  tim  i 
went  on  a  jet  stream  developed  in  the  upper  level  of  the  model  atmosphere, 
and  easterlies  were  generated  in  the  lower  layers  in  the  north  and  in  the 
south.  Lows  and  highs  evolved  at  the  surface,  and  became  displaced  to  the 
north  and  to  the  south  respectively.  Thus,  the  model  was  capable  of  genera¬ 
ting  some  of  the  basic  features  of  the  general  circulation  of  the  stmosphere. 

It  was  considered  essential  to  introduce  vertical  advection  of 
radon.  The  model  allows  a  diagnostic  computation  of  the  vertical  p-velocity, 
^  ,  at  333  mb  and  666  mb,  denoted  as  levels  2  and  4.  These  velocities 
were  used  to  advect  radon  vertically,  uo  was  assumed  to  vanish  at  p  =  0 
and  at  p  =  1000  mb.  Thus  it  is  not  required  to  specify  the  radon  concentra¬ 
tion  at  the  top  or  at  the  bottom  of  the  model  atnosphere.  It  is  quite  obvious 
that  the  introduction  of  vertical  advection  required  that  the  horizontal  advec¬ 
tion  should  be  not  only  due  to  the  rotational  component  of  the  wind,  but  also 
due  to  the  divergent  component,  in  order  that  the  total  mass  of  the  material 
advected  should  be  conserved  (except  for  decay  processes).  The  divergent 
components  were  obtained  by  relaxation  from  the  uj  -  distribution. 

The  prognostic  equation  for  concentration  of  radon  was: 

Of  ■  '  Oyl  )  AC  s  )  (1) 

in  which:  c  =  concentration  m  atoms  per  ton, 
x,  y  -  the  horizontal  coordinates 
u,  v  -  the  horizontal  components  of  wind, 
p  -  pressure, 
cu  -  dp/dt, 

\  -  the  decay  constant, 

S  -  the  source  constant. 
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(5) 
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where  ■£  ,  i  and  j  identify  the  discrete  values  in  time,  x  and  y,  and 


) 


) 

(6) 

Cyclical  continuity  was  assumed  in  the  x  direction.  The  vertical  advectio.n 
and  the  divergence  term  were  neglected  at  the  northern  boundaries.  These 
boundaries  were  treated  from  the  point  of  view  of  the  dynamics  involved,  as 
"solid"  walls.  Time  increments  of  1  hour,  space  increments  of  444  km, 
and  a  grid  size  of  15  x  15  points  were  used. 

The  source  at  level  5,  introduced  into  the  atmoaphere  1100  atoms 

2 

per  ton  per  sec.  ,  which  is  equivalent  to  CCuu  atoms  per  m  per  sec.  Two 
different  source  distributions  were  used:  (a)  A  uniform  distribution  over 
the  whole  region,  (b)  A  uniform  distribution  over  the  western  half  of  the 
region  defined  as  the  "continent",  while  the  eastern  half,  without  sources, 
was  defined  as  an  "ocean".  A  third  numerical  integration  (c)  was  carried 
out  with  the  same  source  distribution  as  (b)  but  with  different  initial  condi¬ 
tions. 

In  (a)  and  (b  )  radon  was  introduced  on  the  10th  day  of  the 
dynamic  experiment,  while  in  (c)  it  was  introduced  on  the  16th  day. 
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Results  and  Discussion 

The  radon  distribution  was  allowed  to  evolve  between  the  10th 
day  and  the  23rd  day  and  between  the  16th  day  and  the  23rd  day  of  the 
circulation  experiment  for  cases  (a),  (b)  and  case  (c)  respectively. 

3  days  were  selected  for  a  detailed  analysis:  the  14th,  the  18th  and  the 
22nd  days.  These  days  are  characterized  by  different  distributions  of  the 
vertical  velocity.  Fig.  1  to  6  show  the  radon  distribution  at  levels  5  and 
3,  superimposed  on  the  vertical  wind  field  at  level  4  on  the  corresponding 
days. 

On  the  14th  day  the  circulation  is  characterized  by  well  defined 
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centers  of  the  convergence  and  divergence.  The  area  with[u>|  7  10  cb/sec 
is  quite  extensive.  The  developing  circulation  up  to  this  stage,  acting  on 
an  initially  uniform  concentration  introduced  at  day  10,  results  in  a  rather 
regular  distribution  of  radon.  The  iso-concentration  lines  of  radon  are 
roughly  parallel  to  the  lines  of  constant  vertical  wind,  and  the  maxima  of 
the  radon  concentrations  at  levels  5  and  3  coincide  with  the  maxima  of  up¬ 
ward  vertical  wind.  This  regular  pattern  is  most  pronounced  in  experiment 
(a)  where  the  source  is  uniformly  distributed  over  the  total  region,  but  even 
in  experiment  (b),  where  the  center  of  the  convergence  is  located  over  the 
"ocean",  i.  e.  over  the  region  which  does  not  introduce  radon  into  the  air, 
tongues  of  radon  are  directed  towards  this  center. 

Day  18  has  a  different  distribution  of  vertical  velocities,  with 

tv/o  separate  cells  of  relatively  small  magnitudes  of  the  vertical  wind.  The 

-4 

area  over  which  |u>|  >  10  cb/sec  is  very  small,  and  in  one  of  the  cells 

-4 

does  not  exceed  0.  3  x  10  cb/sec.  As  a  result  horizontal  transports  become 
dominant,  as  can  be  seen  from  the  elongation  of  the  iso-concentration  curves 
of  radon  due  to  the  zonal  flow. 
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Day  22  is  again  characterized  by  relatively  high  magnitudes  of 
the  vertical  wind.  Here  again  the  correspondence  between  the  vertical  wind 
and  the  radon  concentration  is  quite  remarkable.  Fig.  7  presents  the  dis¬ 
tribution  of  ration  in  experiment  (c)  superimposed  on  the  radon  distribution 
in  experiment  (b).  The  correspondence  between  the  two  is  indicative  of  the 
short  "memory"  of  radon;  its  distribution  is  almost  unaffected  by  the  pre¬ 
vailing  flow  pattern  at  the  moment  at  which  radon  is  introduced  first  into  the 
model  atmosphere. 

The  meridional  distribution  of  the  zonally  averaged  concentration 
of  radon  is  shown  in  Figs.  9  -  14.  A  pronounced  correlation  with  the  vertical 
velocity  was  observed  on  the  days  with  strong  ascent  and  descent  of  air  Thus 
on  days  14  and  22  the  zonally  averaged  vertical  wind  at  level  4  has  a  strong 
effect  on  the  zonally  averaged  distribution  of  radon.  The  distribution  observed 
on  day  18  on  the  other  hand,  was  quite  erratic. 

Fig.  15  gives  the  ratio  between  the  total  radon  content  of  the  west¬ 
ern  half  of  the  region  and  the  total  content  of  the  eastern  half.  In  experiment 
(a),  where  the  source  is  uniformly  distributed,  the  fluctuations  reflect  the 
changes  in  the  convergence-divergence  patterns  generated  by  the  model.  In 
experiment  (b),  with  the  sources  located  in  the  western  half  of  the  region,  the 
trends  are  somewhat  similar  to  those  in  experiment  'a),  but  the  ratio  is  ob¬ 
viously  lower.  The  relative  amounts  of  radon  at  different  levels  at  the  end  of 
the  integration  period  (23rd  day)  are  presented  in  table  1,  in  the  form  of  per¬ 
centages. 
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Sum  m  ary 

Thu  pL'i'lim iiuiry  numerical  experiments,  described  above^  generated 
lime  dependent  radon  distributions,  which  show  that  the  radon  distribution  res¬ 
ponds  to  large  scale  patterns  of  the  vertical  wind,  as  well  as  to  those  of  hori¬ 
zontal  flow.  This  correspondence  was  quite  pronounced  in  the  case  of  a  well- 

defined  cellular  distribution  of  oo  .  In  our  model  the  maximum  vertical  p- 

-4 

velocity  generated  did  not  exceed  t  2  x  10  cb/sec.  Danielsen's  analysis  (1966) 

-4 

resulted  in  values  up  to  t  5  x  10  cb/sec.  It  seems  reasonable  to  suppose  that 
velocities  of  this  magnitude  will  have  an  even  more  marked  effect  on  the  radon 
distribution  Advection  by  small-scale  motion  will  modify  the  distributions 
generated  by  synoptic-scale  flow.  However,  it  is  quite  difficult  to  incorporate 
such  effects  into  present-day  numerical  models.  The  results  obtained  so  far 
are  schematic,  and  it  is  planned  to  repeat  the  calculations  with  more  realistic 
general  circulation  models.  Small  scale  effects  could  be  introduced  parame;r- 
ically  In  any  case,  the  results,  presented  above,  raise  the  possibility,  that  ob¬ 
served  distributions  of  radon  in  the  free  atmosphere  may  help  us  in  the  identi¬ 
fication  of  the  synoptic-scale  patterns  of  vertical  motions. 
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Figs  1  -  6  :  The  distribution  of  radon  and  of  the  vertical  velocity. 

-  -  -  -  lines  of  constant  concentration  of  radon  at  833  mb. 

-  .  -  -  lines  of  constant  concentration  of  radon  at  500  mb. 

- lines  of  constant  vertical  velocity. 

Hatched  areas  are  regions  of  rising  air.  Cross-hatched  areas  are 

-4 

regions  with  oo  >10  cb/sec.  M  and  m  refer  to  maximum 
upward  and  downward  flow. 

The  heavy  vertical  line  in  Figs.  2,  4  and  6  indicates  the  boundary  between 
the  region  containing  sources  of  radon  (on  the  left)  and  the  region  without 
sources. 

Fig.  1:  Day  14,  Experiment  a. 

Fig.  2:  Day  14,  Experiment  b. 

Fig.  3:  Day  18,  Experiment  a. 

Fig.  4:  Day  18,  Experiment  b. 

Fig.  5:  Day  22,  Experiment  a. 

Fig.  6:  Day  22,  Experiment  b. 

Fig.  7:  The  zonally  averaged  zonal  component  of  the  velocity,  as  a 

function  of  latitude  and  altitude,  in  m  sec  .  The  center  of 
the  diagram  models  the  45°  latitude.  Hatched  areas  correspond 
to  easterly  flow. 

Fig.  8:  The  distribution  of  radon  on  the  22nd  day  in  Experiments  (b) 

and  (c), 

- .  lines  of  constant  concentration  of  radon  in  Exp.  (b)  at  833  mb. 

-  lines  of  constant  concentration  of  radon  in  Exp.  (b)  at  500  mb. 

-  -  -  -  lines  of  constant  concentration  of  radon  in  Exp.  (c)  at  833  mb. 

-  -  -  -  lines  of  constant  concentration  of  radon  in  Exp.  (c)  at  500  mb. 

Fig.  9-14:  The  zonally  averaged  vertical  velocity  and  the  zonally  averaged 

radon  concentration  as  a  function  of  latitude. 

-  vertical  velocity.  - radon  concentration  in  Exp.  (a) 

-  radon  concentration  in  Exp.  (b). 

Fig.  15:  The  ratio  between  the  total  radon  content  of  the  western  half 

of  the  region  to  that  of  the  eastern  half  of  the  region,  for  ex¬ 
periments  (a)  and  (b),  as  a  function  of  time. 


Fig.  ? 


latitude 


U) 

a 

D 

cn 

to 

LJ 


CL 


Fig.  8 


DAY  13  500  mb. 


Fig.  12 


DAY  22  833  mb. 


?■ 


Cb/sec]  l  IO"4Cb/saci 


rr.~ 


h 


DOCUMENT  CONTROL  DATA  -  R  &  D  I 

--  -  ‘  at  *hf»,  t;ul\  abxtrmct  and  indewin^  annotation  must  bo  entered  whan  the  overall  reoort  la  ch  ..allied) 

. !ebi>-i.  University  of  .iemsalem 
■  a h vi,  Israel 

2a.  REPORT  SECURITY  CLASSIFICATION  1] 

UNCLASSIFIED 

26.  GROUP 

a 

-  **om  •  f  "  .  *  ■  1  '  Sj 

Numerical  Studies  ot'  Planetary  Circulations  in  a  Model  ' ~mosphere  | 

i  SC  *•  i  f*  *  •  v  *■  no  Tf.i  r  . . i  report  and  inclusive  dates)  f 

Scientific.  Pinal.  \ 

•  fhO'cv  tif sf  aj.  ir,  ■  .J.lio  mlttal,  last  name)  |( 

vo  ra  it  am  lluss  | 

"OM,;A"  March  1970 

7«.  TOTAL.  NO.  OF  PACES  76.  NO.  OF  RETS  ‘1 

233  66  £ 

.»  U  au  T  a 40 

AF  Fa\R  67-40 

“'WO  .  t  v  •  NO 

8004-04 

61445014 

t>S  lain 

Oa.  ORIGINATOR***  REPORT  NUM BER|S) 

9h.  other  report  nois)  (Any  other  number a  that  may  be  aaaiftned  jj 

thla  report)  j 

MSTInbUtlON  ST  *  T  LMKN  T 

This  document  has  been  approved  for  public  release  and  sale;  j 

its  distribution  is  unlimited. 

12.  SPONSORING  MILITARY  ACTIVITY  j 

Air  Force  Cambridge  R;  search  Laboratories  j! 

CCRHD)  L  G  Hanscom  Field,  Bedford,  Mass  ] 

01730  .! 

- i! 

fixtures  of  the  planetary  circulation  over  a  limited  region.  In  the  first 
experiment  the  results  obtained  from  a  primitive  equation  2  level  model  are 
compared  with  the  results  of  a  geostrophic  2  level  model.  The  equations  were 
integrated  over  a  -plane  with  two  lateral  walls  in  the  north  and  in  the  south. 
In  the  second  experiment  the  results  of  a  3  level  geostrophic  model  are 
compared  with  those  of  a  2  level  geostrophic  model.  In  the  third  experiment 
the  smoothing  effects  of  a  lateral  viscosity  term  are  investigated  and  compared 
with  the  smoothing  effect  of  the  Matsuno  approximation,  of  the  time  derivative. 
In  the  4th  experiment  the  spread  of  radon  in  a  model  atmosphere  is  investigated 
schematically,  the  source  of  the  radon  being  located  in  the  lower  boundary. 


DD  ‘”•.1473 


UNCLASSIFIRD 


Security  Classification 


General  Circulation 
Primitive  Equatioi  Model 
|.:i  oral  Viscosity 
Matsuno  Scheme 

Harmonic  Analysis  of  the  Geopotential  Field 
Release  of  Latent  Heat 
Transport,  of  Radon 


